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PREFACE 


; This NASA conference publication contains the proceedings of the Third 

International Symposium on the Science and Technology of Low Speed and Motor- 
less Flight held at the NASA Langley Research Center, Hampton, Virginia, 

March 29-30, 1979. The symposium was cosponsored by the Langley Research 
Center (LaRC) and the Soaring Society of America (SSA) . Oran Nicks, Deputy 
Director of the Langley Research Center, and James Nash-Webber, Massachusetts 
Institute of Technology and past chairman of the SSA Technical Board, were 
general cochairmen. Perry Hanson, NASA LaRC, was the symposium organizer and 
technical program chairman. Hewitt Phillips, NASA LaRC (Retired) ; Joseph Gera, 
NASA LaRC; and Robert Lamson, Chairman of the SSA Technical Board, served as 
chairmen for the technical sessions. 

The purpose of the Symposium was to provide a forum for the interchange of 
information on recent progress in the science and technologies associated with 
low speed and motor less flight. Twenty-eight papers were presented in the areas 
of low speed aerodynamics, new materials applications and structural concepts, 
advanced flight instrumentation, sailplane optimal flight techniques, and self- 
launching and ultralight glider technology. This NASA conference publication 
contains these presentations and a paper, which was not presented, on proposed 
' definitions for various categories of sailplanes and gliders. 

The use of trade names or manufacturer's names in this publication does 
not constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by NASA. The included papers are largely as submitted. 

The physical quantities, whether in the International System of Units (SI) or 
■ U.S. Customary Units, are retained as submitted by the authors. 
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LOW-SPEED SINGLE-ELEMENT 


AIRFOIL SYNTHESIS 

John H. McMasters and Michael L. Henderson 
Boeing Commerical Airplane Co. 

SUMMARY 

Large quantities of experimental data exist on the characteristics of 
airfoils operating in the Reynolds number range between one and 
ten million, typical of conventional atmospheric wind tunnel operating 
conditions. Beyond either end of this range, however, good experimental 
data becomes scarce. Designers of model airplanes, hang gliders, 
ultralarge energy efficient transport aircraft, and bio-aerodynamicists 
attempting to evaluate the performance of natural flying devices are hard 
pressed to make the kinds of quality performance/design estimates taken 
for granted by sailplane and general aviation aerodynamicists. Even 
within the usual range of wind tunnel Reynolds numbers, much of the data 
is for "smooth" models which give little indication of how a section will 
perform on a wing of practical construction. 

The purpose of this paper is to demonstrate the use of recently developed 
airfoil analysis/design computational tools to clarify, enrich and extend 
the existing experimental data base on low-speed, single element 
airfoils, and then proceed to a discussion of the problem of tailoring an 
airfoil for a specific application at its appropriate Reynolds number. 
This latter problem is approached by use of inverse (or "synthesis") 
techniques, wherein a desirable set of boundary layer characteristics, 
performance objectives, and constraints are specified, which then leads 
to derivation of a corresponding viscous flow pressure distribution. In 
this procedure, the airfoil shape required to produce the desired flow 
characteristics is only extracted towards the end of the design cycle. 
This synthesis process is contrasted with the traditional "analysis" 
(either experimental or computational) approach in which an initial 
profile shape is selected which then yields a pressure distribution and 
boundary layer characterisitics, and finally some performance level. The 
final configuration which provides the required performance is derived by 
cut-and-try adjustments to the shape. 

Examples are presented which demonstrate the synthesis approach, 
following presentation of some historical information and background data 
which motivate the basic synthesis process. 

INTRODUCTION 

Since the dawn of human flight, enormous efforts have been expended on 
the design of efficient wings and their constituent airfoil sections. As 
such development became a race for ever increasing speed, the problems of 


low-speed flight frequently became relegated to the status of 
"off-design" conditions, with performance requirements met by fitting 
"high speed" cruise airfoils with increasingly complex and sophisticated 
high-lift devices. During the past forty years, relatively little 
attention has been given to the development of "optimized" low-speed 
airfoils by other than academicians and "cut-and-try" experimenters. 

While frequently outside the mainstream of modern commerical interest, 
the range of low-speed flying devices (characterized by generally low 
values of the scale parameters Reynolds and Mach number) covers an 
enormous portion of the feasible flight spectrum. To place the 
subsequent discussion in a proper global context. Figure 1 has been 
prepared to demonstrate quantitatively the relationships between 
low-speed flight vehicle size and performance and the sometimes arcane 
parameter, Reynolds number. While "low-speed" generally implies low 
Reynolds and Mach numbers, it is worth noting that recent interest in 
ultralarge transport aircraft has now expanded the low-Mach number flight 
Reynolds number range from that typical of small insects (10<Rn<104) 
through devices like huge wing-in-ground effect aircraft (ref. 1) which 
may have chord Reynolds numbers approaching one billion at flight speeds 
on the order of 100 m/s (M'^0.3). Even a "small" monster like the Boeing 
747 (average wing chord approximately 10 m) becomes a low-speed aircraft 
during approach, with typical average Reynolds numbers for the wing of 40 
million at M*^0.2. 

To discuss the full range of problems associated with wing/airfoil design 
for the range of vehicles shown in Figure 1 would require several books. 
The present paper is limited to a discussion of two aspects of the 
overall problem: 

1. A brief survey of historical trends in low-speed, single-element 
airfoil development, culminating in a review of the present 
state of the art in analytic design methodology. 

2. A demonstration of the value of modern computational 
capabilities to, first, clarify the performance characteristics 
of several existing low-speed airfoil sections for which 
experimental data exist; and then show how one may proceed to 
"synthesize" a suitable section for a specific application from 
a desired specification of boundary 1 ayer/pressure distribution 
characteristics. 


NOTATION 


c 


c 


AR 


b 


Aspect ratio = b/c = b^/S 
Wing span (m) 
chord (m) 

Average chord = S/b (m) 
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Cj Section drag coefficient 

Skin friction coefficient 
Wing lift coefficient = lift/qS 
Section lift coefficient 
Cn Pressure coefficient = (p-p_)/q 

Section pitching moment coefficient 
H Boundary layer form parameter =6/0 

M Mach number 

p Static pressure (N/m2) 

q Dynamic pressure = (N/m) 

Rn Reynolds number = Vc/ v 

5 Wing area (m2) 

t Airfoil thickness (m) 

V Velocity (m/s) 

V Local velocity (m/s) 

W Weight (N) 

X Chordwise coordinate 

z Coordinate normal to chord 

Greek symbols: 

a Angle of attack (degrees) 

(5 Boundary layer displacement thickness = f (l - ^)dz 

Voo 

e Section lift-drag ratio = Cq /C , 

X. Q ^ 

6 Boundary layer momentum thickness = 

V Kinematic viscosity (1.46 x 10"5 m^/s standard sea level) 

p Air mass density (1.225 kg/m3 standard sea level) 

Superscript: 

( )* indicates "design condition" 
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Subscript: 


( )r recovery point or region 

( )tr transition point or "trip" location 
( )fp fair point (see Fig. 9) 

( )TE trailing edge 

( ) 00 free-stream condition 

( )u airfoil upper surface value 

HISTORICAL PERSPECTIVE 

To clarify the present status of low-speed airfoil development, it is of 
interest to briefly review the history of how we got from there to here. 

A map of the route is shown in Figure 2 . It is important to note that 
well into the present century airfoil "design" was a largely empirical 
process, drawing its main inspiration from natural models (i.e., birds), 
and only partially clarified and systematized by recourse to potential 
flow theory (e.g., Joukowski airfoils). Elaborate testing programs at 
Gottingen and by the NACA, among others, guided by intuition, experience, 
and inviscid theory eventually lead to the accumulation of masses of data 
and subsequent publication of airfoil section catalogs to aid designers. 

It was not until the mid-1930's that the influence of viscous "scale 
effects" was appreciated, and boundary layer theory well enough 
developed, to allow the qualitative incorporation of viscous flow 
concepts into the design of "low-drag" sections. The main upshot of 
these new considerations was the famous NACA 6-series "laminar flow" 
airfoils. The accumulated results of fifty years of empiricism 
culminating in the matrix of 6-series sections are covered extensively in 
the classic catalogs by Abbott and von Doenhoff (ref. 2), Riegels (ref. 

3) and reports such as those by Jacobs and Sherman (ref. 4). 

The preeminence of the 6-series sections (slightly altered on occasion to 
the taste of the individual designer) lasted for nearly twenty years, and 
these sections have only been overshadowed since the late 1950 's by the 
emergence of the revolution ushered in by the computer. While the 
equations of advanced potential flow methods and viscous flow theory can 
be concisely written, it is quite another matter to routinely solve 
analytically the complex flow fields around even "simple" airfoils in a 
real fluid. Thus, until the advent of large computers, theory could only 
guide what remained a largely experimental development effort. 

The wind tunnel is a marvelous tool for describing what happens, but 
seldom provides much guidance on why a particular event (e.g., boundary 
layer separation) occurs. To go beyond the level of "design by testing," 
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practical quantitative solutions to the equations of viscous flow were 
required to supplement empirical experience. 

The remarkable success of computer based methods in improving airfoil 
performance beyond the NACA 6-series level is well demonstrated in the 
catalog of Wortmann FX-series sections (ref. 5) and the reports and 
papers listed in refs. 6 and 7. Despite this new progress, designers 
without access to a computer of sufficient size, or those lacking a 
sophisticated background in theoretical aerodynamics and mathematics are 
still forced to rely on catalog data and outmoded "simplified" theory. 
With very few exceptions (notably ref. 8), available good catalog data i 
for "ideal" surface quality wind tunnel models operating in the range 
7 X io5 < Rn < 107. As a summary of the preceding historical discussion 
Figure 3 shows some representati ve airfoils sections used, or 
specifically designed for, various categories of low-speed aircraft 
during the last eighty years. The variety of shapes even within a given 
category is sometimes bewildering. 

LOW-SPEED AIRFOIL DESIGN 

The general principles of low-speed, single-element airfoil design in 
light of modern theory have been discussed in detail by several authors, 
notably Wortmann (ref. 9-11), Mi ley (ref. 12) and Liebeck (ref. 13). A 
brief review is presented here in Appendix A. 

Whether one is designing a new airfoil section or attempting to select 
one from a catalog, it is important that all the relevant criteria are 
kept clearly in mind. The author's list is as follows: 

Basic Airfoil Selection/Design Criteria 

1. Basic Operating Conditions (superscript * indicates design point): 

a. Lift Coefficient Range (o < < C^ < 

drag 

b. Reynold Number Range Rn f^i^max) 

c. Mach Number Range ^0 < M*< 

2. Airfoil Characteristics Desired (Priorities to be established for 
each specific application): 

a. Low Drag (e.g., absolute minimum drag at C^ , "low" drag over 

operating C^^^range). ^ I' 

b. High Lift (e.g., absolute C , moderate C. with "gentle" 

stall). ^max ^max 


c. 


Pitching Moment (e.g., positive moment for flying wing 
applications, low negative moment to minimize horizontal tail 
trim loads or aeroelastic effects on wing). 


3. Practical Constraints: 

a. Required thickness-chord ratio and/or required local structural 
thickness. 

b. Anticipated surface quality (e.g., skin joints or slat/airfoil 
junctions which might force boundary layer transition). 

High-Lift/Low Drag Design 

From the preceeding list it can be seen that the airfoil selection/design 
process is complex and this partially accounts for the wide variety of 
section shapes shown in Figure 3, each intended to strike some 
particularly beneficial compromise between often conflicting 
requirements. It is seldom possible to state categorically that a 
particular section is the "best" one even for a given type of aircraft. 

Within the overall low-speed performance spectrum, however, one is 
generally forced to bias the selection/design toward achievement of 
either: (a) low-drag, or (b) high-lift. No general rules can be given 

for how much "high-lift" one can achieve with a "low-drag" section or 
vice versa, although clues are beginning to emerge from modern viscous 
flow theory. General guidelines for good design can be formulated, and 
these are briefly reviewed in Appendix A. 

It should be noted that the NACA 6-series airfoils are basically "low 
drag" sections. Their long reign is due more to the fortuitous fact that 
they scaled well with Mach number, rather than providing the long runs of 
laminar flow which was the original design objective. Only in the 
special case of applications to sailplane wings was the original 
objective met, practical construction and operational problems (bugs, 
paint, rivets, dimples, etc.) tending to abort the "laminar flow" 
behavior in other applications. None of the 6-series sections can be 
categorized as "high-lift" airfoils. 

Empirical Data 

With the preceeding list of airfoil selection/design criteria in mind, 
one can consult the various catalogs to see if a suitable section 
exists. Data from these standard sources (e.g.. Refs. 2-5, 7, 8) is 
summarized in global terms in Figure 4. 

Within the range of Reynolds number for which large quantities of data 
exist, a diligent searcher can find some apparently curious anomalies - 
specifically the "spectacular" Liebeck sections (ref. 13). That the 
Liebeck sections achieve the high-lift performance shown is no longer in 
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serious question, nor are the reasons such performance is achieved. What 
remains unclear is the nature of the trade-offs in section 
characteristics which are available between the "feasible upper bound" 
represented by the Liebeck sections and the "top-of-the-line" 
conventional sections within the shaded bands shown in Figure 4. 


As a prerequisite to discussion of systematic methods for evaluation of 
these trade-offs, some appreciation of the parameters of boundary layer 
theory as they relate to airfoil performance is required. Figures 5 
through 8 show some examples of the boundary layer characteristics of 
several familiar sections and the relationships between this data; and 
the more traditional display of global performance data, section geometry 
and pressure distributions is discussed in detail in Appendix B. 

AIRFOIL SYNTHESIS 

To advance beyond an empirically based approach to airfoil selection, or 
to consider the prospect of tailoring airfoil sections to a specific 
application, it is necessary to understand the difference between a 
design approach based on "analysis" as contrasted with one based on 
"synthesis." The synthesis (inverse) approach to airfoil design begins 
with the boundary layer characteristics as they effect the pressure 
distribution and ultimately define and limit the performance of a section 
in every way. The airfoil shape is derived last in this process, and is 
that physically realizable contour which provides the desired flow 
characteristics. Synthesis is almost the direct opposite to the 
traditional "empirical" (analysis) approach wherein one begins with a 
shape which yields a pressure distribution and a set of boundary layer 
characteristics, and thus initial values of lift, drag and moment. 
Performance requirements are finally met by trial and error modification 
of the shape. Whether these modifications are made to a wind tunnel or 
computer model, the basic process is one of iterative cut-and-try until 
the solution "converges." 

AN INVERSE AIRFOIL DESIGN TECHNIQUE 

While the possibility of synthesizing an airfoil has been recognized for 
many years, it has only been possible to implement satisfactory inverse 
methods (based on modern boundary layer theory) since the advent of the 
computer. Synthesis approaches have been employed by Wortmann (ref. 9) 
and more recently by Liebeck (ref. 13). A very general technique for 
airfoil synthesis (applicable to both single- and multi-element section 
components) has recently been developed by Henderson (ref. 14), based on 
proven integral boundary layer techniques described largely in 
Schlichting (ref. 15). While the specific techniques used in the overall 
program may seem almost old fashioned, the program has proven to be very 
satisfactory in practice and is quite a powerful tool for both single and 
multi-element airfoil synthesis (particularly when coupled with the 


methods described in ref. 16). Details of the method are described in 
reference 14, and only the basic elements are listed here for reference. 

Elements of an Inverse Boundary Layer Analysis and Design Technique 


Component 

Laminar Boundary Layer 
Laminar Separation 
Laminar Separation Bubble 
Transition 

Turbulent Boundary Layer 


Turbulent Separation 
Compressibility Corrections 
Profile Drag 


Theory (ref. 15 except *) 

Pol hausen 
Polhausen 

Henderson (empirical)* 

Granville 

Momentum integral 

Power law velocity profile 

Garner's eqn. for form 

parameter 

Ludwieg-Til Iman eqn. for wall 
shear stress 

H > 3.0 

Karman-Tsi en* 

Squire and Young 


Utilizing the methodology outlined above, it becomes possible to 
implement the airfoil design process shown in Figure 9. Once an 
"optimized" viscous flow pressure distribution and linear theory airfoil 
shape have been determined, the powerful methods described by Henderson 
in reference 16 (which also account for separated flows) are applied to 
arrive at the final airfoil geometry which yields that pressure 
distribution, and final analytic performance predictions are made. 


Several points in this synthesis process need to be clarified. For 
example, any "airfoil" shape will produce a unique pressure 
distribution. The converse is not generally true. In order to assure 
that an initial "designed" pressure distribution will result in a closed, 
non-reentrant airfoil shape, an upper surface pressure distribution is 
designed free of geometrical constraints, and a lower surface pressure 
distribution is defined as that which will result in a section with an 
NACA OOXX thickness form. This yields a total pressure distribution 
vyhich will result in a realizable airfoil of desired thickness. This 
initial lower surface pressure distribution and its corresponding 
boundary layer characteristics are usually poor. In the initial stage, 
however, it is the upper surface which is being optimized, and it is a 
simple matter to subsequently reconfigure the lower surface (guided by 
the preliminary result) to a more desirable form as indicated in Figure 9. 
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The program allows a rather arbitrary specification of upper surface 
recovery region form parameter (H) variation as a primary input. Thus 
one can systematically study the effect of this important parameter 
easily and in some detail before proceeding to more detailed design 
calculations. This feature will be demonstrated shortly. The 
significance of various form parameter variations is discussed in 
Appendix B. 

The most difficult parameter to specify correctly at the outset is the 
trailing edge pressure coefficient. This parameter has a very powerful 
effect on the design lift level a theoretical section will achieve, and 
to date the determination of its final "correct" value has generally 
required an iterative approach. The problem is discussed at some length 
by Liebeck (ref. 13). 

Probably the weakest part of the theoretical performance estimation 
procedure is calculation of profile drag. In principle, at the final 
stage in the design cycle one can integrate the total pressure and skin 
friction drag components and arrive at a total profile drag coefficient. 
Experience to date with viscous flow programs which accurately predict 
pressure distributions and hence lift and pitching moments gives 
generally less accurate drag estimates. This is due primarily to the 
fact that drag is usually two orders of magnitude lower than lift, and 
whereas errors in lift computations are small with a good pressure 
! distribution predictor, errors in pressure integration (particularly in 
the leading edge region) tend to be on the same order as pressure drag 
values. Thus for simplicity, the present state of the art is to rely on 
the method of Squire and Young (ref. 15) for total drag prediction and, 
in the present case, a supplementary calculation of skin friction drag to 
provide a clarification of the magnitude of this component within the 
total drag value. This procedure has been found to be reasonably 
adequate, at least for purposes of comparing the drags of single-element 
sections. While absolute values of Squire and Young drag may sometimes 
be questionable, anyone experienced with the pecularities of 
two-dimensional wind tunnel testing (particularly at high-lift values) 
must realize the magnitude of the error band in "good" experimental drag 
data. 


SOME RESULTS 

To indicate the use of the above methodology, two examples have been 
, chosen to demonstrate several aspects of the influence of Reynolds number 
I on airfoil characteristics. Figure 10 demonstrates the results 

I obtainable from a parametric study of the influence of variations of 

■ recovery point location and Reynolds number on a family of sections with 
simple roof-top pressure distributions (cf. Fig. 9), and a common 
specified exponential form factor variation in the recovery region. The 
principal observations to be made in this example are the significant 
difference in "optimum" recovery point between sections designed (for 


high lift-drag ratios) at two million and thirty million Reynolds number, 
and the ultimate desirability of designing to full-scale Reynolds number 
conditions (i.e., 30 x 106 in this case) to achieve maximum performance, 
despite the fact that such results may appear inferior to those obtained 
from a design optimized at wind tunnel conditions when both are tested at 
low Reynold numbers. 

Figure 11 shows the effect of a systematic variation of recovery region 
form parameter on the shape and characteristics of three airfoils 
designed to the same lift coefficient level at a Reynolds number of 
five-hundred thousand. The performance characteristics of these sections 
are summarized in Figure 12, and clearly show the trades available in 
lift, drag, pitching moment and stall break from different specifications 
of recovery region characteristics. 

The results shown in Figure 12 are generally nonobvious and are of some 
interest in view of the discussion in Appendix B and the fact that 
relatively little modern experimental data exists for sections designed 
specifically for this low value of Reynolds number. The stall behavior 
of the three sections can be understood on the basis of the discussion in 
Appendix B regarding the correlation between boundary layer form 
parameter (H) variation and upper surface separation progression. 

A more subtle and remarkable aspect of the results shown in Figure 12 is 
that the net Squire-Young drag of all three sections at the design point 
lift coefficient is nearly the same. The rate at which the drag rises 
between the design point and maximum lift coefficients will be different, 
however, reflecting the way in which flow separation progresses on the 
three sections as stall is approached. The example calculations also 
show the relative values of upper surface recovery region (turbulent) 
skin friction coefficient relative to the total upper surface profile 
drag coefficient. Although the highly concave recovery pressure 
distribution of Airfoil C (which approaches a Stratford type recovery, 
c.f. Appendix B) shown in Figure 11 has the lowest skin friction 
coefficients, it also has the highest rate of growth (and final trailing 
edge value) of boundary layer momentum thickness. Thus while Airfoil C 
has the lowest skin friction drag it has the highest pressure drag and in 
the overall balance, all three sections exhibit similiar net profile drag 
values. This effect is not limited to the low Reynolds number case 
shown. As Reynolds number increases, the pressure drag becomes the 
increasingly dominant drag term, and minimization of the recovery region 
turbulent skin friction coefficient by employing a Stratford type 
recovery becomes increasingly less satisfactory. 

CONCLUDING COMMENTS 

A review of the history and present state of the art of low-speed 
single-element airfoil design has been presented, leading to a 
description of a powerful new inverse boundary layer scheme which can be 
used to synthesize an airfoil section tailored to the requirements of a 
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specific aircraft. The basic intent of this paper has been to provide 
background and motivation for this alternative approach to airfoil 
design, as contrasted with the more traditional "design by 
experiment/analysis" approach to the problem. Along the way (Appendix B) 
it has been possible to clarify the performance characteristics of 
sections of quite different geometry and design objectives, and indicate 
the influence of Reynolds number on both "low-drag " and "high-lift" 
sections. Several examples of parametric analyses using the "synthesis" 
methodology have been presented which only hint at the potential of these 
new techniques. 

It has been shown that airfoil design (even when limited to very low Mach 
numbers and single-element sections) is a hugely complex problem to which 
no single "best" solution exists even for a single specialized category 
of aircraft type. On the other hand, it is clearly possible to derive a 
section biased and optimized to the taste of an individual aerodynamic! st 
with a great deal more intelligence than was possible less than a decade 
ago. Much work still needs to be done, however, to finally free the hang 
glider designer from reliance on his present very slender catalog of 
airfoil candidates. 
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APPENDIX A: BASIC AIRFOIL DESIGN 

The purpose of this appendix is to provide a brief tutorial review of 

some of the principles of airfoil design. The discussion follows that of 

Wortmann (ref. 11), Miley (ref. 12) and Liebeck (ref. 13). 


All practical airfoils will carry some lift loading (whether high, low, 
or moderate) at some desired operating condition, and this will be 
characterized by generation of some peak level of negative pressure 
coefficient on the upper surface of the section, followed by recovery to 
near free-stream conditions at the trailing edge. The pressure loading 
on the lower surface will depend on factors like required maximum section 
thickness, establishment of favorable pressure gradients for low-drag at 
the section design lift level, and the requirements of satisfactory 
"off-design" performance at low section lift coefficients. At some point 
on both surfaces of the contour, the initial run of laminar boundary 
layer flow will transition to turbulent flow, the particular transition 
points being strongly dependent on the Reynold number, the form of the 
pressure distribution (or the profile shape which generates it), the 
surface quality of the section, and the free-stream turbulence level. 

All other factors being equal, the natural transition point will move 
forward on the profile as Reynolds number increases. 

At this point there is a parting of the ways as one seeks either 
high-lift, or low-drag performance at low-to-moderate lift coefficients. 

To achieve low-drag, the longest possible runs of laminar flow are 
desired on both surfaces of the section followed by an orderly transition 
to thin turbulent boundary layer flow as the pressure recovers to 
trailing edge conditions; and separation is to be avoided like the plague. 

In the high-lift case, attention mainly focuses on the upper surface. As 
in the low-drag case laminar flow is sought, together with high negative 
pressures over the forward portion of the section. The problem in the 
high-lift case is not necessarily to delay the onset of turbulent flow, 
but rather to cause an orderly transition at some optimum point to a 
healthy thin turbulent boundary layer over the pressure recovery region 
to allow the flow to decelerate from the high peak values reached on the 
forward portion without significant separation. The "optimum" high-lift 
upper surface pressure distribution will thus be constructed to produce 
the highest possible loading on the forward portion of the profile, 
consistent with the recovery capability of the turbulent boundary layer 
beginning at an "optimum" transition point. At low Reynolds numbers, 
getting rid of laminar flow at the recovery point and avoidance of large 
scale laminar separation become a major consideration. 

A major constraint on the high-lift section is the character of the stall 
break; all things being equal, a gradual stall progressing from the 
trailing edge is desired. It should also be noted that the bulk of 
' "good" high- lift sections achieve their maximum lift coefficients after 
upper surface (trailing edge) separation has begun. Controlled laminar 
separation bubbles may even be tolerated vf they lead to orderly 
transition to turbulent flow in the pressure recovery region and do not 
burst before trailing edge separation is well developed. 

In the high-lift case, the lower surface pressure distribution will be 
tailored in much the same fashion as in the low-drag case, although the 


lower surface pressure distribution can be made to produce a significant 
portion of the net lift and/or alter the pitching moment 
characteristics. This factor and the influence of various forms of upper 
surface distribution on section pitching moment coefficients are 
indicated in Figures 9 through 12 and in Appendix B. 

APPENDIX B: SOME RELATIONSHIPS BETWEEN AIRFOIL 

PERFORMANCE AND BOUNDARY LAYER CHARACTERISTICS 

While most aerodynami cists have some appreciation of the section geometric 
parameters (e.g. thickness, camber, leading edge radius, trailing edge 
angle) which may influence performance, relatively few have a 
corresponding "feeling" for the fundamental parameters of boundary layer 
theory (e.g. form parameter, momentum thickness), and how these 
parameters are influenced by scale effects. The purpose of this appendix 
is to provide a brief evaluation of the boundary layer characteristics of 
several representative airfoils, and a description of how these 
parameters relate to the more familiar presentations of pressure 
distributions and global performance characteristics. An understanding 
of the connection between boundary layer behavior, pressure distribution, 
and section geometry as they influence performance is essential to 
success in the synthesis approach to design. 

The performance characteristics of four familiar sections are shown in 
Figure 5. Two of these sections (the NACA 633-OI8 and Wortmann FX 
61-184) have been designed primarily for low-drag, and the other two (the 
FX 74-CL6-140 and Liebeck L1003) for high-lift. These sections actually 
represent something of a continuum in that the NACA section is a classic 
"minimum drag" shape while the Liebeck is a pure "high-lift" section. 

The Wortmann FX 61-184 (ref. 5, 11) is a classic 1960 vintage sailplane 
section designed for "low-drag" over a "wide" range of lift coefficients, 
with a compromise struck between absolute low drag, thickness, and a very 
benevolent stall behavior at a moderate maximum lift coefficient. 

The FX 74-CL6-140 (ref. 18) on the other hand, represents an attempt to 
design a section with the same level of maximum lift coefficient as the 
Liebeck, but with a biased compromise again being struck between 
thickness, maximum lift, wide "drag bucket" and satisfactory stall 
characteristics. All four sections are quite different in shape, and in 
the absence of detailed information on the types of pressure distribution 
and boundary layer characteristics (including an evaluation of the 
post-separated flow region) one is provided only superficial clues to why 
each of these sections exhibits such different performance characteristics. 

As an aside, the influence of flow separation on the performance of a 
section and the importance of accurately modeling this effect in a 
theoretical design exercise have been graphically demonstrated by 
Henderson (ref. 16). Figure 6 shows an experimental lift curve for the 
NASA GA(W)-1 section (ref. 17) in comparison with theoretical 
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calculations made with increasingly sophisticated analytical techniques. 

For this particular section. Figure 6 shows that modeling the attached 
boundary layer flow remains inadequate in predicting the variation in 
lift with angle of attack beyond 75% of the final maximum lift 
coefficient value. The full theory developed by Henderson (ref. 16), 
which models both the boundary layer and separation, provides excellent 
predictions however. This improved methodology (which extends to 
multielement sections) represents a major, and so far unique, advance in 
computational capability. 

To better understand the differences in performance and shape between the 
sections shown in Figure 5, it is necessary to evaluate in detail the 
pressure distributions and boundary layer parameter (specifically the 
form parameter, H) variations for each section. Example data for the 
NACA 633-OI8 (ref. 19) at 2° angle of attack (within the drag bucket of 
the section) are shown in Figure 7 for three widely different Reynolds 
numbers. The classic 6-series aft-end shape corresponds to a roughly 
linear rise in the recovery region pressure distribution, and consequent 
form parameter (H) variation shown. The influence of Reynolds number on 
the location of the point of natural transition is indicated, and clearly 
shows the difficulty of achieving long runs of laminar flow as Reynolds 
number increases. 

As shown in Figure 11, the shape and magnitude of the form parameter (the 
ratio of boundary layer displacement thickness to momentum thickness) 
variation in the pressure recovery region of the airfoil correlate in 
general with the shape of the pressure distribution in this region. The 
specification of recovery region form parameter variation is one of the 
central inputs in the Henderson inverse method described previously. As 
discussed in Schlichting (ref. 15), laminar separation occurs when H 
reaches 3.5 and turbulent separation begins when H exceeds about 3.0. 

The influence of the H-factor variation on airfoil stall behavior will be 
discussed presently. 

Wortmann (refs. 9-11) has argued that there are advantages to a "concave" 
recovery pressure distribution (with near constant value of recovery 
region form parameter) for drag reduction, compared to the linear or 
convex pressure distributions associated with earlier profiles, including 
many of the Gotti ngen/Joukowski airfoils (c.f. Figure 3). The basic 
principles of the design of Wortmann's sailplane and related sections 
(including the FX 61-184) with concave pressure rises have been 
thoroughly discussed in references 9 through 11, and by Miley (ref. 12). 
These references also discuss the importance of properly contouring both 
the upper and lower surfaces of low-drag profiles. 

Turning attention to the high-lift airfoils cases, it is interesting to 
compare the pressure distributions and boundary layer characteristics of 
the Wortmann FX 74-CL4-140 (ref. 18) and Liebeck L1003 (ref. 13) shown in 
Figure 8, and contrast this data with that for the NACA 633-OI8 in Figure 7. 
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The Liebeck sections are of great theoretical interest for several 
reasons. Members of the family apparently approach the upper limit of 
lift coefficient achievable with a single-element section without 
mechanical boundary layer control. The sections also exhibit commendably 
low drag coefficients in the region of the design lift coefficient and 
low pitching moments. In exchange for these desirable characteristics, 
the stall behavior is wretched and the undersurface separates at rather 
high (positive) lift coefficients, thus limiting "high-speed" 
performance. This latter factor can be partially ameliorated by use of a 
camber changing trailing edge flap; however, the abrupt stall behavior is 
a fundamental characteristic of the basic family. 

The Liebeck sections have been theoretical ly designed by the previously 
described synthesis process, in this case by use of a Stratford recovery 
region pressure distribution (ref. 20) to establish the maximum level of 
negative pressure on the upper surface "roof top" region of the section. 
The Stratford recovery region pressure distribution is that which, for a 
turbulent flow, results in a boundary layer which is everywhere equally 
close to separation. Thus, to within the accuracy of the Stratford 
formulation, the recovery region boundary layer is either completely 
attached or completely separated - there is no (theoretical) middle 
ground. This factor accounts for the very abrupt stall behavior of the 
sections. Thus, by reliance on the Stratford distribution, Liebeck 
generated the single class of high lift sections which can be "optimized" 
and analyzed without recourse to explicit partially separated flow 
calculations. Herein lies the success Liebeck had in designing to very 
much higher lift coefficients and section lift-drag ratios than had once 
been thought possible for a single-element section. The resulting shapes 
and pressure distributions for Liebeck sections are entirely non-obvious 
and the prospects of happening on them by "cut-and-try" were remote. 

This example provides a strong motivation for use of inverse methods. 

The experimental verification of the predicted performance of the Liebeck 
sections, and by extension the validation of the Stratford theory, 
apparently opens a whole new prospect in high-lift airfoil design. 
However, the inability of Liebeck 's methodology to account for partially 
separated flows, and the resulting formal reliance on the Stratford 
distribution, severely circumscribe the range of sections which can be 
designed. The possible trade-offs in performance between the Liebeck 
sections and the range of conventional sections shown in Figure 4 
remain obscure. 

The result of a highly sophisticated attempt to design such an 
"intermediate" airfoil, which trades some drag and thickness for a better 
stall behavior, while acheiving the same high-lift level, is represented 
by the Wortmann FX 74-CL(X)-140 pair discussed in ref. 18. Referring to 
Figure 8, one sees that the Liebeck and Wortmann pressure distributions 
are quite different, although both have "concave" distributions in the 
recovery region. Where Liebeck uses a well defined "instability" region 
as described by Miley (ref. 12) to achieve orderly transition to 
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turbulent flow in the recovery region, Wortmann forces the formation of a 
"well-behaved" thin laminar separation bubble which acts as a passive 
boundary layer trip. 

Reviewing the performance curves for the Wortmann and Liebeck high-lift 
sections shown in Figure 5, one sees the consequences of the two 
approaches to the design problem. Looking at the resulting airfoil 
shapes and pressure distributions in Figure 8, one sees little in common 
between the two sections however. To see how "equally" high-lift 
coefficients are generated by two such dissimilar sections, one must 
refer to the details of the boundary layer characteristics for the two 
airfoils. 

For both the Liebeck and Wortmann sections, recovery begins at about 40% 
of the chord aft of the leading edge. Prior to this, the "laminar H" for 
the Liebeck section is nearly constant through the instability region, 
falling abruptly to an initial "turbulent" value as the flow 
transitions. By contrast, on the Wortmann section the laminar H rises 
abruptly prior to transition until a value of H for laminar separation is 
reached, following which a "short bubble" is formed leading to transition 
and turbulent reattachment at the beginning of the recovery region. 

Once into the recovery region, the turbulent form parameters on the 
Liebeck section rise rapidly to an initially high value and then begin a 
further very gradual linear rise to a point just short of the trailing 
edge. This recovery region form parameter variation is characteristic of 
a Stratford imposed pressure distribution. 

On the Wortmann section, the turbulent form parameter does not jump 
initially, but rises instead from its starting value behind the laminar 
bubble at a nearly identical rate to that of the Liebeck/Stratford, until 
it hooks upward at the end. The result is again a generally concave 
pressure distribution on the recovery portion of the Wortmann section. 

Comparison of these form parameter variations for two very different 
"looking" sections clarifies much of the difference in stall behavior 
between the sections. On the Liebeck section, as angle of attack is 
increased beyond the "design" value (design lift coefficient equal to 
1.8), the recovery region form parameter level is shifted progressively 
upward until a value of approximately 3.0 is reached, at which point 
turbulent separation begins. With the Liebeck/Stratford recovery 
pressure distribution, the form parameter level is almost constant across 
the bulk of the recovery region. Thus, if nothing else (a laminar short 
bubble for example) interferes, the whole recovery region becomes 
"critical" with respect to separation at nearly the same time, and an 
abrupt stall subsequently occurs. By contrast, the recovery region form 
parameter on the Wortmann section does not reach so uniform a critical 
level as angle of attack is increased towards stall. This is reflected 
in the more gradual stall break for the Wortmann section. The existence 
of the short bubble ahead of the recovery point on the Wortmann section 


throughout this approach to stall clouds the issue of how the stall 
progresses, and the critic will note that the stall behavior is not that 
much better than the Liebeck. That the stall progresses 
non-catastrophically (at least initially) from the trailing edge is 
indicated (c.f. Fig. 5) by the creeping drag rise as stall is approached 
and entered. 

The preceeding examples are intended to be illustrative of a few well 
known sections and demonstrate some specific trends. The results shown 
are not necessarily typical of wide classes of sections and the possible 
ranges of form parameter variation and pressure distribution are 
enormous. These limited examples do, however, demonstrate the level of 
detailed analysis which modern theory can provide, and the necessity of 
delving this deeply into detail in order to understand differences and 
similarities between airfoils with different shapes and global 
performance characteristics, and finally to design an optimized profile 
for a given application. Obviously, much more could and should be said 
on these topics. In addition, much needs to be said regarding the 
problems of "optimizing" both upper and lower surface contours, and the 
influence on drag of form parameter variation, boundary layer momentum 
thickness, transition point, etc. All of these investigations require a 
technique by which the important variables of the problem can be varied 
in an orderly and systematic fashion, particularly as a function of 
Reynolds number. Such a technique has been described in this paper. 
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AN EXPLORATORY INVESTIGATION OF THE EFFECT OF A PLASTIC COATING 

ON THE PROFILE DRAG OF A PRACTICAL-METAL-CONSTRUCTION 

SAILPLANE AIRFOIL 

Dan M. Somers 
Langley Research Center 

SUMMARY 


An exploratory investigation was performed in the Langley 1 ow-turbul ence 
pressure tunnel to determine the effect of a plastic coating on the profile 
drag of a practical -metal -construct! on sailplane airfoil. The model was tested 
with three surface configurations: (1) filled, painted, and sanded smooth; 

(2) rough bare metal; and (3) plastic-coated. The investigation was conducted 
at Reynolds numbers based on airfoil chord of 1.1 x 10^, 2.2 x 10^, and 
3.3 X 10^ at a constant Mach number of 0.10. 

The results indicate that, at all three Reynolds numbers, the order of 
the drag values of the three surface configurations, starting with the highest 
drag, was: filled, painted, and sanded smooth; rough bare metal; and 

plastic-coated. 


INTRODUCTION 


Research on advanced technology airfoils has received considerable 
attention over the past several years at the Langley Research Center. As part 
of this overall research program, the present investigation was conducted to 
determine the effect of a plastic coating on the profile drag of a practical- 
metal -construction sailplane airfoil. Accordingly, a two-dimensional wind- 
tunnel model was constructed by an American sailplane manufacturer employing 
the same sheet-metal fabrication techniques used in constructing the corre- 
sponding production wing. Three surface configurations were investigated: 

(1) as received (filled and painted); (2) bare metal; and (3) plastic-coated. 

The plastic-coating procedure is described in detail in reference 1. The air- 
foil, which corresponds to the FX 67-K-170/17 airfoil designed by F. X. Wortmann, 
is representative of state-of-the-art laminar airfoils having variable geometry 
(in this case, a plain flap). The experimental section characteristics of the 
FX 67-K-170/17 airfoil are reported in reference 2. 

The investigation was performed in the Langley 1 ow-turbul ence pressure 
tunnel (ref. 3). The profile-drag coefficients of the three configurations 
were obtained at Reynolds numbers based on airfoil chord of 1.1 x 10^, 

2.2 X 10^, and 3.3 x 10^ at a constant Mach number of 0.10. The geometric 
angle of attack varied from -5° to 10°. The results have been compared with 
data from reference 2. 
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SYMBOLS 



d 

h 

M 

R 

s 

X 

2 

a 


pressure coefficient 

airfoil chord, cm (in.) 

section profile-drag coefficient. 




point drag coefficient (ref. 4) 
section lift coefficient 

section pitching-moment coefficient about quarter-chord point 
surface-waviness-gage reading, cm (in.) 
vertical distance in wake profile, cm (in.) 
free-stream Mach number 

Reynolds number based on free-stream conditions and airfoil chord 

arc length from leading edge, cm (in.) 

airfoil abscissa, cm (in.) 

airfoil ordinate, cm (in.) 

angle of attack, deg 


MODEL, APPARATUS, AND PROCEDURE 


Model 

The constant-chord wind-tunnel model was constructed by an American 
sailplane manufacturer employing the same sheet-metal fabrication techniques 
used in constructing the corresponding tapered production wing. The structure 
consisted of a spar and four stringers to which a 0.81 mm (0.032 in.) skin was 
flush-riveted. In addition, four ribs were flush-riveted to the skin at 
30.48-cm (12.00-in.) intervals spanwise. The model had a chord 66.47 cm 
(26.17 in.) and a span of 91.44 cm (36.00 in.). A plain lower-surface-hinged 
flap having a chord of 0.17c was fixed at 0° deflection (fig. 1). The flap gap 
was sealed with tape along the lower surface. No orifices were installed in 
the model . 

Three surface configurations were investigated (fig. 2). Configuration 1 
(as received) (fig. 2(a)) had a factory finish - a painted epoxy primer (filler) 
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which had been sanded to insure an aerodynamically smooth surface. Configura- 
tion 2 (bare metal) was obtained by chemically removing the paint and primer. 
(See fig. 2(a).) The surface of configuration 2 (bare metal) was very rough 
because it had been mechanically roughened at the factory to provide a good 
bonding surface for the epoxy primer (fig. 3(a)). A plastic film was then 
bonded to the metal of configuration 2 (bare metal) to obtain configuration 3 
(plastic-coated) (figs. 2(b) and 3(b)). It should be noted that the rough sur- 
face of configuration 2 (bare metal) can be seen through the plastic film and 
adhesive of configuration 3 (plastic-coated) (fig. 3(b)). The thickness of 
the plastic film was approximately 0.1 mm (0.005 in.) whereas the adhesive 
averaged about 0.25 mm (0.010 in.) in depth. The thickness of the plastic film 
and the adhesive together was nearly equal to that of the paint and filler as 
illustrated in figure 2(c). Configuration 1 (as received) and the 
FX 67-K-170/17 airfoil are compared in figure 2(d). The coordinates of the 
three configurations together with those of the FX 67-K-170/17 airfoil are 
listed in table I. 

A relative waviness survey was made at the midspan of configuration 3 
(plastic-coated). (See fig. 4.) A surface-waviness gage as described in 
reference 5 was used. The distance between the feet of the gage was 
approximately 6.4 cm (2.5 in.). 


Wind Tunnel 

The Langley low-turbulence pressure tunnel (ref. 3) is a closed-throat, 
single-return tunnel which can be operated at stagnation pressures from 10.13 
to 1013 kPa (0.1 to 10 atm) with maximum tunnel-empty test-section Mach numbers 
of 0.46 and 0.23, respectively. The minimum unit Reynolds number is approxi- 
mately 0.66 X 10° per meter (0.20 x 10^ per foot) at a Mach number of about 
0.10, whereas the maximum unit Reynolds number is approximately 49 x 10^ per 
meter (15 x 10^ per foot) at a Mach number of 0.23. The test section is 
91.44 cm (3.000 ft ) wide by 228.6 cm (7.500 ft) high. 

Hydraulically actuated circular plates provide positioning and attachment 
for the two-dimensional model. The plates, 101.6 cm (40.00 in.) in diameter, 
are flush with the tunnel sidewalls and rotate with the model. The model ends 
were mounted to rectangular model -attachment plates as shown in figure 5. 

Wake-Survey Rake 

A fixed, wake-survey rake (fig. 6) was cantilevered from the tunnel 
sidewall at the model midspan and approximately 0.9 chords downstream from the 
trailing edge of the model. The wake rake employed 91 total -pressure tubes, 
0.152 cm (0.060 in.) in diameter, and 5 static-pressure tubes, 0.318 cm 
(0.125 in.) in diameter. The total-pressure tubes were flattened to 0.102 cm 
(0.040 in.) for a length of 0.61 cm (0.24 in.) from the tips of the tubes. 

Each static-pressure tube had four flush orifices located 90° apart, 8 tube 
diameters from the tip of the tube in the measurement plane of the total - 
pressure tubes. 
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Instrumentation 


Measurements of the wake-rake pressures were made by an automatic 
pressure-scanning system. Basic tunnel pressures as well as the wake-rake 
pressures were measured with variable-capacitance precision transducers. Geo- 
metric angle of attack was measured by a calibrated digital shaft encoder driven 
by a pinion gear and rack attached to the circular plates. Data were obtained 
by a high-speed data-acquisi tion system and were recorded on magnetic tape. 

Tests and Methods 

The airfoil was tested at Reynolds numbers based on the airfoil chord of 

1.1 X 10^, 2.2 X 10®, and 3.3 x 10® at a Mach number of 0.10 over an angle-of- 
attack range from -5° to 10°. For several test runs, the upper surface of 
configuration 3 (plastic-coated) was coated with oil to determine the location 
as well as the nature of the boundary-layer transition from laminar to turbulent. 

Section lift coefficients and pitching-moment coefficients about the 
quarter-chord were determined with the viscous-flow airfoil method of refer- 
ence 6 because no orifices were installed in the model. Section profile-drag 
coefficients were computed from the wake-rake total and the wake-rake static 
pressures by the method of reference 4. 

Standard low-speed wind-tunnel boundary corrections (ref. 7), approximately 
2 percent of the measured coefficients, have been applied to the drag data. 

DISCUSSION 

Pressure Distributions 

The theoretical chordwise pressure distributions at the approximate limits 
of the laminar low-drag range are shown in figure 7. At an angle of attack of 
00 (c^ = 0.5) which corresponds to the lower limit of the laminar low-drag 
range, a favorable pressure gradient was predicted on the upper surface to about 
x/c = 0.40 whereas a zero pressure gradient was predicted on the forward portion 
of the lower surface. As angle of attack was increased, the calculated pres- 
sure gradient on the lower surface became more favorable whereas that on the 
upper surface became less favorable. At an angle of attack of 6° (c^ = 1.2), 
the upper limit of the low-drag range, a favorable pressure gradient was pre- 
dicted to about x/c = 0.60 on the lower surface whereas a zero pressure gradient 
was predicted on the forward portion of the upper surface. Between the lower 
and upper limits of the laminar low-drag range, favorable pressure gradients 
were predicted on the forward portions of both surfaces. 

Section Characteristics 

The section characteristics of the three configurations are shown in fig- 
ure 8 and tabulated in table II. The lift and drag coefficients of the 
FX 67-K-170/17 airfoil are shown for comparison, having been interpolated from 
the data of reference 2, which were obtained at Reynolds numbers of 1.0 x 10^, 
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1.5 X 10^, 2.0 X 10^, and 2.5 x 10 ^, As previously mentioned, both the lift and 
pitching-moment coefficients of the three configurations were generated by the 
theoretical method of reference 6, which appears to give excellent agreement 
with experiment where no trail ing-edge separation is present (ref. 8). Accord- 
ingly, the c^ - versus - a and the c^ - versus -c^ portions of figure 8 are 
entirely theoretical whereas the c■^ - versus -c^j portion consists of the theo- 
retical lift coefficient plotted against the experimental drag coefficient. No 
quantitative measure of maximum lift coefficient is possible because of a lack 
of separation modelling in the theory of reference 6. 

The mechanism of boundary-layer transition from laminar to turbulent on this 
airfoil at these Reynolds numbers is a laminar separation bubble as shown in fig- 
ure 9 and illustrated in the sketch below. 



The bubble was caused by a slight adverse pressure gradient immediately 
downstream of the minimum pressure on the upper surface. (See fig. 7.) This 
slight adverse gradient was a design feature of the airfoil, as discussed in 
reference 9. 

The section characteristics at a Reynolds number of 1.1 x 10^ are shown 
in figure 8(a). The drag of configuration 1 (as received) was the highest, 
the drag of configuration 2 (bare metal) lower, and the drag of configuration 3 
(plastic-coated) the lowest. The drag coefficients interpolated from the data 
of reference 2 for the FX 67-K-170/17 airfoil fell between those for 
configurations 1 (as received) and 2 (bare metal). 

One possible explanation for the above order, based upon an understanding 
of laminar separation bubbles and the data presented in references 9-11, 
follows. The lower drag coefficients of configurations 2 (bare metal) and 3 
(plastic-coated) have been attributed to reductions in the size of the laminar 
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separation bubble on the upper surface of the airfoil. These reductions were 
probably caused by two different mechanisms. For configuration 2 (bare metal), 
introduction of additional disturbances into the laminar boundary layer by the 
roughness of the surface (fig. 3(a)) apparently did not cause premature transi- 
tion because they were too small. Once the laminar boundary layer had separated, 
however, the disturbances would grow rapidly, resulting in transition and, 
finally, turbulent reattachment. These additional disturbances, therefore, 
probably reduced the distance between laminar separation and transition (i.e., 
a shorter bubble). For configuration 3 (plastic-coated), introduction of dis- 
turbances into the laminar boundary layer by the waviness of the surface 
apparently affected the length of the laminar separation bubble as did the 
disturbances caused by the roughness of configuration 2 (bare metal) with an 
even shorter bubble for configuration 3 (plastic-coated). The waviness of the 
configuration 3 (plastic-coated) surface (fig. 4) was probably caused by hand 
application of the plastic film on very thin sheet metal. 

6 6 

The section characteristics at Reynolds numbers of 2.2 x 10 and 3.3 x 10 
are shown in figures 8(b) and 8(c), respectively. The drag of configuration 1 
(as received) was again the highest, the drag of configuration 2 (bare metal) 
was lower, and the drag of configuration 3 (plastic-coated) was again the lowest. 
The drag coefficients interpolated from the data of reference 2 for the 
FX 67-K-170/17 airfoil were higher than those for configuration 1 (as received) 
at a Reynolds number of 2.2 x 106 . The explanation for these results is probably 
the same as that for a Reynolds number of 1.1 x 10^. 

Results similar to those described above have been reported by other 
investigators. A substantial drag reduction was obtained by using a trip wire 
to eliminate the laminar separation bubble on the upper surface of an airfoil 
(ref. 9). Reductions in the sizes of the laminar separation bubbles on two 
different airfoils through the introduction of disturbances by roughness and 
trip wires were reported in references 10 and 11, respectively. 


SUMMARY OF RESULTS 


An exploratory investigation was performed in the Langley low-turbulence 
pressure tunnel to determine the effect of a plastic coating on the profile 
drag of a practical -metal -construction sailplane airfoil. The model was tested 
with three surface configurations: (1) filled, painted, and sanded smooth; 

(2) rough bare metal; and (3) plastic-coated. The resulting data have been 
compared with data for the design airfoil (Wortmann FX 67-K-170/17) from 
another low-turbulence wind tunnel. The investigation was conducted at Reynolds 

numbers based on airfoil chord of 1.1 x 10^, 2.2 x 10^ and 3.3 x 10^ at a 
constant Mach number of 0.10. 

At all three Reynolds numbers, the drag of the filled, painted, and 
sanded smooth configuration was the highest, followed by the drag of the rough 
bare metal configuration, and finally the drag of the plastic-coated 
configuration. 
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TABLE I.- AIRFOIL COORDINATES 


(a) Configuration 


(as received) 


[c = 66.4827 cm (26.1743 in.)] 


Upper surface 

Lov/er surface 

x/c 

2/C 

x/c 

z/c 

0.000000 

0.000000 

0.000000 

0.000000 

.000004 

-.000042 

-.000008 

-.000038 

.000497 

.004443 

.000500 

-.002048 

.000993 

.005651 

.001001 

-.003026 

.001490 

.006827 

.001490 

-.003805 

.001987 

.007721 

.002017 

-.004489 

.002980 

.009452 

.002988 

-.005391 

.003973 

.011133 

.003992 

-.006185 

.004982 

.012722 

.004986 

-.006839 

.006980 

.015611 

.006980 

-.007943 

.009968 

.019389 

.009968 

-.009360 

.014957 

.024635 

.014950 

-.011274 

.019939 

.029055 

.019932 

-.013127 

.029907 

.036658 

.029907 

-.015611 

.039879 

.043390 

.039879 

-.017712 

.049843 

.049549 

.049850 

-.019515 

.059814 

.055214 

.059814 

-.021063 

.069782 

.060498 

.069786 

-.022427 

.079754 

.065381 

.079754 

-.023680 

.089725 

.070080 

.089725 

-.024677 

.099693 

.074363 

.099689 

-.025735 

.119633 

.082123 

.119636 

-.027550 

.149540 

.092274 

.149540 

-.029773 

.199386 

.106295 

.199386 

-.032016 

.249225 

.116867 

.249237 

-.033212 

.299087 

.124687 

.299080 

-.034221 

.348922 

.130342 

.348922 

-.035031 

.398773 

.133142 

.398769 

-.035500 

.448627 

.132928 

.448623 

-.034924 

.498466 

.129925 

.498470 

-.033602 

.548301 

.124095 

.548313 

-.031749 

.598163 

.114696 

.598159 

-.028734 

.648017 

.101921 

.648006 

-.024612 

.697845 

.086738 

.697852 

-.020153 

.747699 

.070241 

.747703 

-.015389 

.797538 

.053824 

.797546 

-.010468 

.847392 

.037896 

.847392 

--.006823 

.897247 

.026320 

.897231 

-.001941 

.947089 

.012975 

.947085 

.001242 

.967025 

.008638 

.967029 

.000455 

.976997 

.006629 

.976993 

-.000317 

.986964 

.004569 

.986964 

-.001035 

1.000000 

.001486 

.999828 

-.001284 


TABLE I.- AIRFOIL COORDINATES - Continued 
(b) Configuration 2 (bare metal) 

Ic = 66.4670 cm (26.1681 in.)] 


Upper surface 

Lower surface 

x/c 

z/c 

x/c 

z/c 

0.000000 

0.000000 

0.000000 

0.000000 

.000050 

..000000 

.000046 

.000000 

.000520 

.003604 

.000501 

-.002109 

.000994 

.004632 

.001376 

-.004330 

.001490 

.005568 

.001494 

-.004548 

.001987 

.006370 

.001987 

-.005381 

.002992 

.008117 

.002992 

-.006496 

.003986 

.009749 

.003990 

-.007184 

.004987 

.011342 

.004979 

-.007819 

.006955 

.014227 

.006993 

-.008931 

.009966 

.018075 

.009966 

-.010433 

.014957 

.023238 

.014949 

-.012508 

.019944 

.027568 

.019936 

-.014132 

.029914 

.035069 

.029918 

-.016730 

.039884 

.041807 

.039892 

-.018832 

.049855 

.047856 

.049858 

-.020552 

.059832 

.053542 

.059836 

-.022061 

.069799 

.058827 

.069810 

-.023364 

.079777 

.063746 

.079784 

-.024564 

.089747 

.068255 

.089747 

-.025627 

.099721 

.072543 

.099724 

-.026723 

.119661 

.080262 

.119665 

-.028489 

.149579 

.090385 

.149583 

-.030694 

.199434 

.104536 

.170547 

-.042151 

.249292 

.115190 

.249288 

-.033166 

.299162 

.123020 

.299154 

-.034683 

.349005 

.128653 

.348997 

-.035532 

.398871 

.131511 

. 398886 

-.036090 

.448726 

.131064 

.448745 

-.035731 

.498592 

. 1 28087 

.498603 

-.034404 

.548443 

.122191 

.548443 

-.032326 

.598297 

.112702 

.598312 

-.029295 

.648171 

.099874 

.648167 

-.024966 

.698014 

.084599 

.698010 

-.020380 

.747876 

.067773 

.747865 

-.015947 

.797735 

.051677 

.797735 

-.011487 

.847589 

.036457 

.847593 

-.007341 

.897455 

.024645 

.897436 

-.003424 

.947314 

.012198 

.947314 

-.000910 

.967254 

.007417 

.967262 

-.000657 

.977243 

.005583 

.977232 

-.000703 

.987290 

.002939 

.987271 

-.000734 

.999889 

.000046 

1.000000 

-.000378 


TABLE I.- AIRFOIL COORDINATES - Continued 
(c) Configuration 3 (plastic-coated) 
Ic = 66.4860 cm (26.1756 in.)] 


Upper surface 

Lower 

surface | 

x/c 

z/c 

x/c 

z/c 

0.000000 

0.000000 

0.000000 

0.000000 

-.000023 

.000061 

.000997 

-.002124 

.000508 

.005108 

.001494 

-.003266 

.001016 

.006449 

.002006 

-.003973 

.001494 

.007664 

.002995 

-.005058 

.001998 

.008607 

.003985 

-.005891 

.002991 

.010017 

.004989 

-.006582 

.003988 

.011522 

.006984 

-.007683 

.004978 

.012978 

.009963 

-.009077 

.006972 

.015935 

.014964 

-.011136 

.009960 

.019904 

.019935 

-.012817 

.014960 

.025111 

.029898 

-.015369 

1 .019938 

.029356 

.039873 

-.017470 

.029898 

.036943 

.049844 

-.019213 

.039858 

.043674 

.059815 

-.020710 

.049844 

.049772 

.069786 

-.022085 

.059811 

.055414 

.079746 

-.023235 

.069783 

.060671 

.089713 

-.024339 

.079746 

.065634 

.099696 

-.025417 

.089725 

.070222 

.119634 

-.026991 

.099688 

.074432 

.149544 

-.029302 

.119634 

.082168 

.199369 

-.031399 

.149540 

.092517 

.249228 

-.032599 

.199377 

.106561 

.299069 

-.033432 

.249221 

.117101 

.348917 

-.034211 

.299069 

.125105 

.398742 

-.034830 

.348905 

- .130603 

.448593 

-.034547 

.398749 

.133323 

.498441 

-.033069 

.448593 

.132983 

.548282 

-.031124 

.498445 

.130102 

.598122 

-.028175 

. 548285 

.124181 

.647985 

-.024259 

.598130 

.114786 

.697806 

-.019388 

.647981 

.102202 

.747666 

-.014789 

.697818 

.087146 

.797510 

-.010208 

.747654 

.070237 

.847350 

-.006284 

.797518 

.053909 

.897179 

-.002124 

.847346 

.038440 

.947042 

.000004 

.897194 

.027927 

.967015 

.000531 

.947046 

.014960 

.976944 

.000604 

.966985 

.010823 

.986923 

.000714 

.976956 

.008496 

1.000000 

.001108 

.988837 

.005998 



.999924 

.003687 

1 



TABLE I.- AIRFOIL COORDINATES - Concluded 
(d) FX 67-K-170/I7 airfoil 


Upper 

surface 

Lower surface 

x/c 

z/c 

x/c 

z/c 

0.00000 

0.00000 

0.00000 

0.00000 

.00107 

.00653 

.00107 

-.00217 

.00428 

.01292 

.00428 

-.00514 

.00961 

.02012 

.00961 

-.00815 

.01704 

.02765 

.01704 

-.01057 

.02653 

.03487 

.02653 

-.01321 

.03806 

.04309 

.03806 

-.01580 

.05156 

.05158 

.05156 

-.01827 

.06699 

.06011 

.06699 

-.02062 

.08427 

.06856 

.08427 

-.02282 

.10332 

.07685 

.10332 

-.02490 

.12408 

.08490 

.12408 

-.02682 

.14645 

.09263 

.14645 

-.02856 

.17033 

.09994 

.17033 

-.03011 

.19562 

.10677 

.19562 

-.03146 

.22221 

.11305 

.22221 

-.03261 

.25000 

.11870 

.25000 

-.03354 

.27866 

.12365 

.27866 

-.03425 

.30866 

.12783 

.30866 

-.03474 

.33928 

.13119 

.33928 

-.03499 

.37059 

.13370 

.37059 

-.03501 

.40245 

.13526 

.40245 

-.03480 

.43474 

.13571 

.43474 

-.03435 

.46730 

.13490 

.46730 

-.03365 

.50000 

.13274 

.50000 

-.03272 

.53270 

.12919 

.53270 

-.03155 

.56526 

.12429 

.56526 

-.03012 

.59755 

.11808 

.59755 

-.02844 

.62941 

.11063 

.62941 

-.02654 

.66072 

.10208 

.66072 

-.02437 

.69134 

.09263 

.69134 

-.02187 

.72114 

.08259 

.72114 

-.01896 

.75000 

.07233 

.75000 

-.01572 

.77779 

.06229 

.77779 

-.01236 

.80438 

.05287 

.80438 

-.00913 

.82967 

.04437 

.82967 

-.00625 

.85355 

.03689 

.85355 

-.00386 

.87592 

.03040 

.87592 

-.00197 

.91573 

.01991 

.91573 

-.00037 

.94844 

.01201 

.94844 

-.00124 

.97347 

.00631 

.97347 

-.00105 

.99039 

.00243 

.99039 

-.00044 

.99893 

.00027 

.99893 

-.00005 

1.00000 

0.00000 

1.00000 

0.00000 


TABLE II.- SECTION CHARACTERISTICS 
(a) R ^1.1 X TO®, M = 0.10 


Cl 

a, deg 

"d 

"m 

Configuration 

(as received) 

n?o 

-5.05 

.0130 

-.0978 

.ORO 

-<*. 1 -» 

.0117 

-.0986 

. IPR 

-3.0 7 

.0109 

-.1043 

. 

-?.09 

.0101 

-.1055 

.41 ? 

-1.05 

. OOR? 

-.1 0R7 

.530 

-.00 

.0086 

-.11?0 

.f>5n 

1 .07 

.OORR 

-.115? 

.75R 

?. 0? 

.OOR4 

-.1178 

.R67 

3.07 

.00P4 

-. 1?03 

.974 

4.1? 

.0080 

-.!??<» 

1 .07? 

5.1? 

.008? 

-. 1?31 

1.190 

R.lR 

.0079 

-. 1?70 

1 .?97 

7.13 

.0075 

-.1310 

1.410 

H.?? 

.0078 

-.1345 

1 .51? 

9.?7 

.0079 

-.1366 

1 .RO? 

1 0.?4 

.0090 

-.1390 

Configuration 2 (bare metal) 

-. ORO 

-5.1? 

.01?R 

-.0906 

.035 

-4.?1 

.0117 

-.09?5 

. IR? 

-3.04 

.0110 

-.0950 

.?R0 

-?.ll 

.0101 

-.0965 

.3R3 

-1.03 

.0073 

-.0980 

.477 

.01 

.00R9 

-.1003 

.5R5 

.90 

.007? 

-. 10?3 

,R94 

1 .90 

.0074 

-.1045 

.R03 

3.03 

.007? 

-.1067 

.915 

4.07 

.0076 

-.1091 

1 .0?R 

5.13 

.0078 

-.11?? 

1.139 

R. 14 

.OOR? 

-.1157 

1 .?55 

7.?1 

.0077 

-.1195 

1 .357 

R. 19 

.0071 

-. 1??4 

1.441 

9.19 

.0076 

-. 1?3? 

1 .943 

10. ?7 


-.1?47 

Configuration 3 (plastic-coated) 

-.090 

-5.10 

.01?3 

-.0848 

.0?? 

-4 . OP 

.0110 

-.0870 

.130 

-3.10 

.0100 

-. 0H94 

.?30 

-2.05 

.0090 

-.0913 

.345 

-1.04 

.007? 

-.0940 


.0? 

.0066 

-.0967 

.575 

1.03 

• OORR 

-.099? 

.RP5 

?.04 

.0070 

-.1015 

.795 

3. OR 

.0068 

-.1035 


4 .05 

.0068 

-.105? 

1.010 

5.10 

.0071 

-.108? 

1.1?0 

R. 1 1 

.007? 

-.1118 

1 .?33 

7.1? 

.0069 

-.1153 

1 .340 

8. IR 

.0069 

-.117? 

1 .4?5 

9. IR 

.0074 

-. 1 1H7 

1.5?? 

10.17 

.0115 

-. 1?10 


TABLE II.- Continued 


Cb) R« 2.2 X 10®, M -= 0.10 


Cl 

a, deg 


'■m 

Configuration 1 (as received) 

-.o?o 

- 5.05 

.0096 

-.0979 

.090 

- 4.05 

.0095 

-.1006 

.204 

- 3.02 

.0090 

-.1037 

.307 

- 2.06 

.0083 

-.1062 

.420 

- 1.04 

.0059 

-.1093 

.538 

.01 

.0057 

-.1127 

.855 

1.03 

.0057 

-.1159 

.770 

2.06 

.0056 

-.1193 

.885 

3.09 

.0060 

-.1223 

.995 

4.09 

.0062 

-.1250 

1.105 

5.13 

.0064 

-.1275 

1 .220 

6.19 

.0064 

-.1292 

1.305 

7.14 

.0065 

-.1303 

1.400 

8.20 

.0067 

-.1345 

1.505 

9.32 

.0093 

-.1365 

1.600 

10.29 

.0224 

-.1392 

Configuration 2 (bare metal) 

-.067 

- 5.14 

.0103 

-.0905 

.037 

- 4.20 

.0096 

-.0927 

.170 

- 3.00 

.0094 

-.0955 

.281 

- 1.96 

.0085 

-.0977 

.378 

- 1.01 

.0063 

-.0996 

.487 

-.02 

.0056 

-.1020 

.613 

1.04 

.0056 

-.1054 

.727 

2 . 0 b 

.0054 

-.1078 

,839 

3.05 

.0056 

-.1096 

.936 

4.07 

.0060 

-.1113 

1.039 

5.10 

.0063 

-.1132 

1.154 

6.15 

.0066 

-.1164 

1.260 

7.17 

.0062 

-.1190 

1.352 

8.26 

.0084 

-.1220 

1.442 

9.30 


-.1232 

1.527 

10.18 


-.1250 

Configuration 3 (plastic-coated) 

-.089 

- 5.08 

.0098 

-.0848 

.025 

- 4.06 

.0092 

-.0874 

.140 

- 3.03 

.0088 

-.0899 

.246 

- 2.03 

.0082 

-.0922 

.360 

- 1.00 

.0061 

-.0950 

.473 

.01 

.0053 

-.0977 

.590 

1.05 

.0053 

-.1008 

.708 

2.08 

.0053 

-.1038 

.817 

3.06 

.0055 

-.1065 

.930 

4.10 

.0057 

-.1087 

1.035 

5.15 

.0057 

-.1100 

1 .143 

6.18 

.0058 

-.1120 

1.247 

7.15 

.0058 

-.1152 

1.333 

8.17 

.0061 

-.1175 

1.417 

9.16 

.0107 

-.1193 

1.518 

10.18 


-.1215 


TABLE II.- Concluded 


Cc) R 3.3 X 10®, M = 0.10 



Cl , deg 

«=d 

'=m 


Configuration 

(as received] 


-.027 

- 5.09 

.0084 

-.0977 

.090 

- 4.03 

.0081 

-.1007 

.542 

-. 0 ? 

.0050 

-.1130 

.657 

1.00 

.0050 

-.1164 

.775 

2.02 

.0050 

-.1197 

.893 

3.08 

.0053 

-.1230 

1.005 

4.08 

.0055 

-.1260 

I. Ill 

5.09 

.0060 

-.1288 

1 .207 

6.11 

.0063 

-.1307 

1.293 

T .13 

.0060 

-.1295 

1.389 

8.14 

.0076 

-.1330 

1 .487 

9 . IS 

.0093 

-.1365 

1.595 

10.28 


-.1395 

Configuration 2 (bare metal ) 

-.067 

- 5.14 

.0093 

-.0906 

.053 

- 4.06 

.0085 

-.0932 

.158 

- 3 . 1 ? 

.0084 

-.0955 

.275 

- 2.09 

.0078 

-.0980 

.394 

-.99 

.0066 

-.1008 

.498 

.01 

.0050 

-.1030 

.610 

1.00 

.0049 

-.1058 

.737 

2.07 

.0047 

-.1094 

.849 

3.08 

.0049 

-.1116 

.950 

4.09 

.0052 

-.1127 

1 .060 

5.17 

.0057 

-.1150 

1.160 

6.12 

.0058 

-.1170 

1 .255 

7.19 

.0056 

-.1185 

1.340 

8.18 

.0099 

-.1210 

1 .437 

9.23 


-.1237 

1.533 

10.23 


-.1255 

Configuration 3 (plastic-coat 

ed] 

-.099 

- 5.14 

.0089 

-.0847 

.025 

- 4.06 

.0081 

-.0875 

.138 

- 3.08 

.0079 

-.0910 

.248 

- 2.08 

.0075 

-.0926 

.367 

- 1.00 

.0061 

-.0955 

.477 

. 0 ? 

.0046 

-.0983 

.596 

1.04 

.0048 

-.1015 

.710 

2.03 

.0047 

-.1045 

.827 

3.06 

.0047 

-.1074 

.942 

4.10 

.0050 

-.1102 

1 . 05 ? 

5.10 

.0053 

-.1130 

1.150 

6.14 

.0058 

-.1130 

1.245 

7.16 

.0057 

-.1149 

1.328 

8.18 

.0073 

-.1170 

1.420 

9.19 

.0308 

-.1195 

1.523 

10.19 

.0799 

-.1219 
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Configuration I (os received) 




Figure 2.- Comparisons of configurations 1 (as received), 2 (bare metal), 3 (plastic-coated), 

and FX 67-K-170/17 ordinates. 



(c) Configurations 1 (as received) and 3 (plastic-coated) . 



Figure 2.- Concluded. 


'X) 



(a) Configuration 2 (bare metal) . 

Figure 3.- Surfaces of configurations 2 (bare metal) and 3 (plastic-coated) 



(b) Configuration 3 (plastic-coated) . 
Figure 3.- Concluded. 
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Figure 4.- Surface waviness of configuration 3 (plastic-coated) 



Tunnel side walls 


Diam. = 1.53 c 






Airflow 


A 

1 _. 


Circular plate 

S2//777777A 
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Top view 


A 

1.38 c 




Model attachment 



Figure 5.- Airfoil model mounted in wind tunnel. All dimensions are in 
terms of model chord, c = 66.47 cm (26.17 in.). 
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Total-pressure probe 



(a) Drawing of wake-survey rake. All dimensions are in terms 
of model chord, c = 66.47 cm (26.17 in.). 


Figure 6.- Wake-survey rake 


(b) Photograph of wake-survey rake in the Langley low-turbulence 

pressure tunnel. 
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Figure 6.- Concluded. 



Figure 7.- Theoretical chordwise pressure distributions for configuration 1 (as received) 
at lower and upper limits of laminar low-drag range for R = 2.2 x 1 q 6 and M = 0.10. 



(a) R = 1.1 X 1 q6; m = 0,10. 


Figure 8 


Comparisons of section characteristics of configurations 1 (as received) 
2 (bare metal), 3 (plastic-coated), and FX 67-K-170/17 airfoil. 
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(b) R = 2.2 X 106 ; M = 0.10 
Figure 8 .- Continued. 



(c) R = 3.3 X 1q6; m = 0.10. 


Figure 8.- Concluded. 
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Leading edge 


(a) R = 1.1 X iq6, m = 0.07, and a = 0°. 

Figure 9.- Oil flow photographs of upper surface of configuration 3 

(plastic-coated) . 
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Leading edge 


Laminar 


Turbulent 


Flap joint 


Trailing edge 


Laminar separation 
Turbulent reattachment 


(b) R ~ 1.5 X iq6, m = 0.10, and a = 0°. 
Figure 9.- Continued. 
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Figure 9.- Continued. 




(d) R » 2.5 X 1q6, m = 0.16, and a = 0°. 


Figure 9.- Concluded. 
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OPTIMUM TAIL PLANE DESIGN FOR SAILPLANES 


Kay Mayland 

Technische Hochschule Darmstadt 


SUMMARY 


Classical drag equations in a modern version have been used 
to calculate the influence of tail modifications on the drag of a 
standard class sailplane. The profile drag v/hich depends on the 
Reynolds number is included in the calculations. Minimum drag is 
compared with real drag for two lift coefficients. 

Some results have no clear tendency but low tail area and 
relatively low tail aspect ratio give some advantages. Optimum 
and real lift ratios between wing and tail plane are compared for 
the original sailplane. 


INTRODUCTION 


Since the energy crisis in 1973 there is a lot of interest in 
reducing the trimmed drag of airplanes (Refs. 1-5). One contribu- 
tion to the trimiTied drag is the wing/tail interference drag. This 
interference drag had been interpreted as a component of the t^il 
lift vector due to local downwash angle at the tail position. Sachs 
(Refs. 6, 7, 8, 9) has shown that this interpretation is not cor- 
rect. The exact method is to calculate the interference drag with 
the aid of the downwash angle at downstream infinity. This new 
explanation corresponds to the well-known biplane theory of Prandtl 
and Munk (Ref. 10) which was also used in some new papers (Refs. 4, 
11, 12). This theory in the modern version was used in this paper 
to show the relation between optimum and real load distributions 
between wing and tail. 

Another purpose of this paper is to show the influence of 
tail plane design on total drag. It is important that the 
addition of the Reynolds number dependent profile drag has great 
influence on the optimum design. 

All calculations are performed for a standard class 
sailplane . 
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SYMBOLS 


Values 

are given in SI units. 

AR 

aspect ratio 

b,b^ 

wing span, tail span, m 

★ 

drag coefficient 

drag coefficient (tail profile drag + total induced 
drag) 


lift coefficient 


pitching moment coefficient 

c 

mean aerodynamic chord, m 

D 

drag, N 

g 

2 

acceleration of gravity, m/s 

h,h , /h. 
wb t 

distance in chord lengths from leading edge of wing 
to c.g., wing-body aerodynamic center, tail 

k . , k, 
wb t 

induced drag factor for wing-body, tail 

L 

lift, N 

m 

mass, kg 

q 

dynamic pressure, N/m^ 

Re 

Reynolds number 

s,s^ 

wing and tail area, m^ 

e 

downwash angle 

^00 

downwash angle at dov/nstream infinity 

e* 

downwash factor, e = e* k , C. 

^wb 

M 

span ratio b^/b 

a 

interference factor 
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Subscripts: 


min 

minimum 

o 

zero lift 

opt 

optimum 

t 

tail 

wb 

wing-body 


BASIC RELATIONS 


The following fundamental relations v/ere used, assuming that 
the_aerodynamics are linear and that the dynamic pressure ratio 
is q^/q = 1: 



C 


m 



(h 

^wb 


h V.) 

wb 



0 


( 1 ) 


( 2 ) 


= C_ + C, (h - - C, ^ (h, - h.,^) = 0 


'm mo , ■'L 

wb 


S t wb 


2 t 2 

+ k , + -5- (C^ + k. + e ) 

D Do , wb L , S Do, t L. “ L. 

wb wb t t t 


(2a) 

(3) 


The last term within the round brackets of the drag equation (3) 
is the wing/tail interference drag. The derivation was given by 
Sachs and shall not be repeated here. This interference drag is 
the product of the tail lift and the dov/nwash angle at downstream 
infinity and corresponds to the one given by Prandtl but the new 
expression is much easier to use in calculations. The downwash 
angle at downstream infinity may be expressed as (Ref. 13) 


>^wb U* 0 - 3) (4) 

wb 

The downwash factor e* = 1 corresponds to a rectangular, e* = 2 
to an elliptic and e* = 3 to a parabolic spanwise lift 
distribution of the wing. 
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Influence of Reynolds number 

Wortmann (Ref, 14) has designed and measured a lot of 
excellent profiles for sailplanes. The wind tunnel test results 
are published in the "Stuttgarter Profilkatalog" . Figure 1 shows 
some test results demonstrating the influence of Reynolds number 
on profile drag for several profiles. The solid lines are 
according to the relations used in this paper for calculating 
the influence of Reynolds number: 


'Do 


wb 


'DO, 


0.009 


Re 


0.3 


'wb 
0.007 


Re 


0.3 


Re in millions 


(5) 


These relations are only valid for the above mentioned profiles 
and for a special Reynolds number range. 


Real lift ratio 


( 2 ) 


For balance in equilibrium flight {C-^ = 0) equations (1) and 
can be solved for the lift ratio between tail and wing 




=const 
= 0 



mo 


wb 


(h - 




( 6 ) 


mo 


wb 


" 'h - *'wb> 


- (h^ - 


h k) 
wb 


This lift ratio depends on fixed quantities and on the parameters, 
total lift coefficient and c.g. position. It is possible to elimi- 
nate the c.g. position by considering the stability requirement. 
The static margin may be expressed as 


ac 

Jc 


m 

L 


(h - h,b) 


^ S 

±J 

a 


- 




Solving equation (7) for (h - h^^^) and combining this with 
equation (6) gives the following equation: 


(7) 
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mo 


wb 


■'wb 


9C S 

rri ott t/'uu \/T 

jc; " -57 — ^ ■ 3 ^r> 

L 


Ct =const 

c^=o 

m 


mo 


wb 


3C 


L S 

^ ^ -CT^ -# (h-W (l-lf)-(Vhwb) 

Lj Li 

a 


( 8 ) 


Regarding equation ( 8 ) it is possible to say that the real load 
distribution between tail and wing depends on several fixed values 
and on the parameters, stability margin and total lift coefficient 


Optimum load Distribution 

Prandtl has published the optimum tail/wing lift ratio in 
his biplane theory (Ref. 10) : 


wb opt 



with y 



(9) 


The interference factor a has to be taken out of diagrams (Ref. 6 ) 
Therefore it is easier to use the following relations. Comparing 
Prandtl 's equation for the interference drag 


- 2 ^ ^t ^wb 
% 5 b,. b 


( 10 ) 


with equation (3) it is possible to rewrite equation (9) as 


wb opt 



(11) 


The optimum tail lift may be either positive or negative 
depending on the downwash factor e*. Only an elliptic spanwise 
lift distribution over the wing requires zero tail lift. The 
combination of equations ( 6 ) and ( 11 ) gives the optimum c.g. 
position 


mo 


(h - 


h u) 4 . 

wb opt 


wb 


1 - — 
2 


(h. - 


(— ) 


^wb^ 


( 12 ) 
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and the optimum stability margin 


9C, 


m 


(h - h ,v.) 


^ S 
at t 



9C 


L opt 


wb^opt ~ S ^ t 


a 


The minimum induced drag was also given by Prandtl and may be 
rewritten in the following equation 


The term in the brackets demonstrates the decreasing induced drag 
of the wing/tail combination compared with the wing alone. 


All calculations are performed with the data set of a typical 
standard class sailplane; see table 1. 

The pilot of this airplane wants to know whether he can 
reach good performance by choosing the correct c.g. position.. 

In figure 2 the equations (8) and (11) are evaluated. It is easy 
to see that the optimum load distribution depends on the 
downwash factor e* but not on the lift coefficient while the 
real lift ratio has inversed dependencies. The downwash factor 
is normally not known exactly; therefore it is not easy to reach 
exact conclusions from figure 2. A more accurate way is to eval- 
uate the equation (13); see figure 3. It is easy to see that 
the optimum stability margin for lift coefficients > 0.4 is 
obtained by a normal (stable) c.g. position. This statement is 
valid for downwash factors between 1 and 2. Good performance for 
a wide range of lift coefficients will therefore be obtained by 
choosing a medium or forward c.g. position. Only high speed flight 
requires aft c.g. positions. 


While the pilot of a sailplane is interested in the optimum 
c.g. position, the sailplane designer is interested in the tail 
design to meet stability requirements and to achieve good per- 
formance. For a fixed wing geometry it is possible to vary two 
main parameters: tail area and tail span. It is important not to 



2 


(14) 


NUMERICAL RESULTS 


Optimum c.g. position 


Tail modifications 
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evaluate the parameter variations only by regarding the total 
induced drag because the tail profile drag will also change. 
Therefore the criterion is the sum of both the total induced drag 
and the tail profile drag (equation 5) 


* 




(15) 


The induced drag is calculated as in equation (3) with the lift 
coefficients obtained from the real load distribution (equation 8). 
The stability margin is assumed to be “ -0.15, and the 

downwash factor is assumed to be e* = 1. (The calculations were 
also performed for e* = 2; the tendencies correspond to e* = 1.) 

The real drag coefficient is compared with the theoretical 
minimum drag coefficient using equation (14) for the minimum 
induced drag. 

In figure 4 this drag coefficient is plotted versus the tail 
span with constant tail m.a.c. for two lift coefficients. Param- 
eter in this diagram is the tail profile drag at Re = lO^. The 
minimum induced drag (C^q. = 0) is decreasing with increasing tail 
span but the higher profile drag due to increasing tail area is 
predominant. It is suitable to design the tail with the minimum 
possible area to satisfy stability requirements. 

Assuming a minimum tail area of 1 m^, another question is, 
what span or what AR^ is optimum? There are two effects: 

- With increasing tail span the minimum induced drag of the 
complete sailplane is decreasing 

- Increasing tail span means decreasing Reynolds number resulting 
in increasing tail profile drag. 

In figure 5 it is shown that the superposition of these two 
effects results in no clear tendency. With increasing profile drag 
the Reynolds number effect becomes predominant. Assuming a tail 
profile drag coefficient of Cdo^- = 0.01 at Re^ = 10® a reduction 
of tail span from 2.4 m (original value) to perhaps 2.0 m will 
give some little advantages. 

The diff erences betv/een minimum and real drag coefficients 
are generally small; only low total lift coefficients (C^ = 0.2) 
require high tail downloads (see figure 2) resulting in greater 
differences . 


Wing and tail modifications 

The wing geometry is included in the variations. The only 
restrictions are now a wing span of 15 m and a total area (wing + 
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tail) of 
wing and 
suitable 
than the 


11 m2. To simplify the calculations it is assumed that 
tail have the same profile drag at Re = 10^. It is 
to regard the drag or the drag areas directly rather 
drag coefficient: 


o 

q 


Sc S 




Re 


0.3 


Re 


0 . 3 ^ 


(16) 


D . 
1 


min 


k , C? S 
wb L 


(EL_2.) i 

q b 



(17) 


Figure 6 shows the total profile drag area (equation 16) plotted 
versus the ratio of tail area to total area. For normal tail 
aspect ratios between 3 and 6 the total profile drag will 
increase with increasing tail area due to decreasing medium 
Reynolds number of the total area. 

The combination of equations (16) and (17) gives the minimum 
total drag which is plotted versus the ratio of tail area 

to total area in figure 7. The trends are clear: the lowest 
possible drag is obtained with low ratios of tail area to total 
area . 


CONCLUSIONS 


The influences of c.g. position and of tail plane design on 
the performance of a standard class sailplane have been shown. 

One important result is that the optimum c.g-. position is for a 
wide range of lift coefficients within the normal c.g. range. The 
calculations for the tail plane design have shown that the 
reduction of induced drag due to higher tail span is less 
important than the influence of profile drag. Low tail area and 
relatively low tail aspect ratio will give some advantages. It 
is remarkable that the best standard class sailplanes of today 
have a tail area of 1 m^ and a tail aspect ratio of AR^ = 5 

while older sailplanes have for example = 1.5 m^ and AR^ = 6. 
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TABLE 1.- DATA SET OF A TYPICAL SAILPLANE 
[Reference lo] 

Mass, kg 300 

Wing : 

Span, m 15 

2 

Area, m 10 

Mean aerodynamic chord, m 0.67 

Aspect ratio 22.5 

Pitching moment coefficient, zero lift -0.1 

Horizontal tail: 

Span, m 2.4 

2 

Area, m 1.0 

Aspect ratio 5.76 

Distance between wing and tail aerodynamic centers, m 3.85 

Stability margin -0.15 
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DRAG COEFFICIENT 



Fig. 1 : Influence of Reynolds number on profile drag 


0.06 



Fig. 2; Lift ratio between wing and tail 
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Fig, 4 : Variation of tail span and area 
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Fig. 5: Variation of tail span, tail area = const 
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Fig. 6: Profile drag area versus tail area 



Fig. 7 : Minimum drag area versus tail area 
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THE EFFECT OF DISTURBANCES ON A WING 


RICHARD EPPLER 
UNIVERSITAT STUTTGART 
STUTTGART, WEST GERMANY 


SUMMARY 


Disturbances such as flap and aileron hinges and poorly faired spoilers 
were simulated in a computer wind tunnel. The total drag of a single 
roughness element does not depend only on the size of that element. Its 
position on the wing has a surprisingly strong effect. In particular, a 
roughness element on the convex side of a deflected flap or aileron causes 
a very siabstantial increase in drag. Ve 2 ry few experimental data are 
available for comparison. Good agreement with experiment can be achieved, 
however, by adapting a fictive "step size." The correlation between the 
real roughness-element size and the drag increase remains to be determined. 
Simple, fundamental experiments are suggested which will allow a 
theoretical estimation of the drag increase due to roughness elements. 


INTRODUCTION 


Disturbances on a wing due to flap and aileron hinges , variable chord 
arrangements, poorly faired spoilers, etc. become more significant as 
airplanes become more efficient. In other words, as the profile drag 
decreases, so must the parasitic drag. Performance differences between 
airplanes of the same type and performance differences between similar 
types have been measured several times. These differences indicate that 
some airplanes have parasitic drag due to seemingly insignificant details. 
It is necessary to investigate such details in the sense of Bruce 
Carmichael's study "What Price Performance?" (ref. 1). As long as we 
spend lots of money on variable chord concepts, we should at least be 
sure to take every opportunity to realize less expensive performance 
improvements. One such improvement could be the reduction of disturbances 
connected with flap and aileron hinges and with spoiler gaps and steps. 
These two-dimensional disturbances usually occur at wing positions where 
the boundary layer is already turbulent. Because there is very little 
experimental information on such disturbances, a theoretical disturbance 
model has been developed. It yields relative effects and indicates which 
simple, fundamental experiments are necessary to obtain a method for the 
estimation of the absolute amount of these effects. 
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THE DISTURBANCE MODEL 


The boundary-layer flow in the region surrounding a two-dimensional 
disturbance of height h perpendicular to the wall is shown in figure 1 . 
The velocity u{y) in the boundary layer at y = h is called Uj^. It is 
plausible that the influence of the disturbance will depend mainly on h 
and uj^. This influence will be evident several step heights downstream 
(fig. 1). As long as h is not too large, the velocity u(y) depends 
only on the wall shear stress Tq. This was shown by Ludwieg and Tillman 
(ref. 2) and reconfirmed by Kader and Yaglom (ref. 3) . A good 
approximation for u(y) is 


u(y) 


u 


T 


yu^ 

a log + b 

^ V 


( 1 ) 


whe re v 


is the kinematic viscosity. 



is the wall shear-stress 


velocity, and p is the density; a (~5) and b (~6.5) are constants. For 
use in a boundary-layer computation method, it is better to transform 
equation ( 1 ) by means of the local skin-friction coefficient 
T 

into 

pu 


^ [2.17 In ^ ^ 


( 2 ) 


where U is the local potential-flow velocity, R = — ^ is the overall 

Reynolds nxamber of the flow, is the free-stream flow velocity, and L 

is the reference length which, for a wing, is the chord length c. 

As long as this approximation is valid for u^, the influence of the 
disturbance will depend only on the local disturbance Reynolds number 

V 

V 


It is even plausible that the displacement thickness 6 ^ will change 
linearly with Rj^. This means that the distance between the undisturbed 
and disturbed velocity distributions is proportional to uj^h. The same is 
approximately true for the momentum thickness 62 and the energy thickness 
(S3. Therefore, the model to be used assumes that 62 is increased by a 
value. 
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(3) 


A62 = ku^h 


as a result of a disturbance of height h. It can be shown that the 
additional assumption 


A6 


3 



(4) 


is also reasonable provided u^^ is not too large. The most difficult 
problem is determining the value of the proportionality constant k. The 
value of this constant will depend on the precise shape of the disturbance. 
It should not be too difficult to obtain accurate values for k from 
simple experiments. Such experiments have been planned by F. X. Wortmann 
and D. Althaus of the Institut fiir Aerodynamik und Gasdynamik at the 
Universitat Stuttgart. 

If the boundary layer is laminar at the position of the roughness 
element, the computation predicts transition at that position. 


DISTURBANCES CW A WING 


A value of k = 0.15 was used in equation (3) as a rough approximation 
for a simple roughness element like a trip wire. A computer program was 
used to evaluate the effects of disturbances on a wing. Given the airfoil 
shape , the program computes the velocity distributions corresponding to 
the various input angles of attack. For all velocity distributions, 
boundary- layer computations are performed for the different input Reynolds 
numbers. Disturbances can be specified at up to two different positions 
on each surface of the airfoil. 

Several examples illustrate the capabilities of this disturbance model. 
For the first example , one disturbance of height h = 1 mm was introduced 
at various positions along the upper surface of airfoil E603. The velocity 
distributions for various lift coefficients are shown in figure 2. 

The roughness element was introduced at three different chordwise 
positions; x/c = 0.4, x/c = 0.6, and x/c = 0.81. The theoretical polars 
at R = 1 X 10^, which corresponds roughly to low-speed flight in a 
sailplane, and R = 3 x 10^, which corresponds to high-speed flight in a 
sailplane, are presented in figure 3. The results clearly indicate that 
the potential- flow velocity at the position of the roughness element has 
a strong influence on the drag. The additional drag nearly always 
increases with lift coefficient as does the local velocity. The polar 
for the most forward disturbance, however, has a different character. 

At C| « 1.1, the most upstream roughness element has less influence than 
the more downstream ones. By looking into the details, it was determined 
that not only the potential-flow velocity, but also the skin-friction 
coefficient at the position of the roughness element has a strong influence. 
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If transition occurs in an adverse pressure gradient and the roughness 
element is shortly behind this transition, the turbulent boundary layer 
will not be fully developed at the position of the roughness element 
and the effect on the drag will, therefore, be quite small. Of course, 
the effect on the drag will be quite large if the roughness element 
shifts the transition point toward the leading edge- 

These results are entirely theoretical. It is, of course, desirable 
to obtain a correlation between theory and experiment. It would be easy 
to perform wind-tunnel experiments which correspond to this example. 

As previously mentioned, these experiments are planned but have not yet 
been performed. 

Experimental data applicable to this problem are rare. There have 
been many experiments concerning the influence of roughness elements on 
transition, but few on the effects of roughness elements on a boiindary 
layer which is already turbulent. One such experiment was performed in 
1971 by D. Althaus in a low-turbulence wind tunnel at the Universitat 
Stuttgart. In that test, the polars of an airfoil (FX 62-K-153/20) 
with a conventional, center-hinged flap (gap sealed) were measured first. 
Then, the polars were determined for the same airfoil with a so-called 
"Elastic Flap" (ref. 4) . These experiments are valuable for evaluating 
the theory because the two models differed only by the radius of the arc 
between the forward portion of the airfoil and the flap and by the single 
step (roughness element) which is a part of every flap hinged in the 
conventional manner. 

The envelopes of the polars for both configurations are shown in 
figure 4. Each envelope was obtained by plotting the lowest drag 
coefficients for the various flap deflections at a given lift coefficient. 
This means that the drag coefficients for zero or negative (up) flap 
deflections aria used for low lift coefficients and the drag coefficients 
for positive (down) flap deflections are used for high lift coefficients. 
Thus, the envelope is defined by the data for the optimal flap deflections. 
The differences between the two curves is quite small for the lower lift 
coefficients and surprisingly large for the higher lift coefficients. 

For some time, no explanation could be found for this apparent 
anomaly. After introducing the disturbance model into the computer 
program, however, it was not difficult to analyze these two configurations 
theoretically. The velocity distributions for the FX 62-K-153/20 airfoil 
(ref. 5) using the original coordinates are shown in figure 5. Not unlike 
many Wortmann airfoils, the coordinates are not smoo'oth. The velocity 
distributions show irregularities, the worst one occurring at the leading 
edge on the lower surface. In the practical use of this airfoil and for 
the wind-tunnel model, these irregularities have probably been smoothed 
out. Therefore, it was reasonable to smooth the coordinates before 
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proceeding with further computations. The boundary-layer method is very 
sensitive to such irregularities, especially with regard to the prediction 
of transition. The velocity distributions for the smoothed airfoil with 
0° and 10° flap deflection are shown in figure 6. The differences between 
the two flap configurations (plain and elastic) are evident only on the 
upper surface in the region around the hinge. The elastic flap causes a 
much lower suction peak at the hinge than the normal, plain flap. Moreover, 
the plain flap introduces a certain, single roughness element at that 
point. It was not possible to specify the height h of the roughness 
element in the disturbance model directly from the step height at the hinge. 
Instead, several different values for h were tried. The theoretical 
results for h = 0.6 mm at R = 1 x 10^ and 3 x lo^ are shown in 
figure 7. Interestingly, the theory shows exactly the same phenomenon 
as the experiment. At 0° flap deflection, the differences between the two 
flap configurations are small, and at 10° deflection, the same roughness 
element causes a considerable drag penalty for the plain flap. For positive 
flap deflection, the roughness element is located precisely at the 
position of the suction peak, which means that it is in a region of high 
potential- flow velocity. Also, the increased favorable pressure gradient 
for the plain-flap configuration ahead of the element causes an increase in 
wall shear stress which further amplifies the drag penalty. For 0° flap 
deflection, the potential-flow velocity and the wall shear stress are much 
lower at the position of the roughness element, which explains the small 
difference for this case. 

Some of the experimental data from figure 4 are included in figure 7. 
The drag penalties predicted by the theory agree well with the experiment. 

It must be emphasized that the absolute value of the drag penalty is 
not the significant result. This value was achieved by selecting the 
right value for h. The ratio between the drag differences with and 
without flap deflection, however, must be pointed out as a fundamental 
result which agrees well with experiment. This result, of course, has 
practical applications and can eventually explain some of the performance 
differences between similar airplanes. The order of magnitude of these 
drag differences should not be neglected in performance calculations. 

The maximum lift coefficient and its decrease due to the roughness 
element were not predicted as well as the drag penalty. It should be 
noted, however, that the wind-tunnel results for maximum lift coefficient 
must be suspect. The wind tunnel used has a closed, rectangular test 
section with the model spanning the tunnel from wall to wall. The lift 
coefficient is determined by measuring the pressures along the top and 
bottom tunnel walls. Thus, these measurements yield the average 
over the span of the model. Note also that the span is roughly equal 
to the chord. Accordingly, any separation at the jxancture of the tunnel 
wall and the model influences the measured lift coefficient considerably. 


85 


CONCLUSIONS 


Single roughness elements have been theoretically modeled. The data 
from a previously run experiment on a flapped airfoil with and without a 
disturbance at the hinge was used for comparison. The drag penalty 
predicted by the theory and measured in the experiment was large enough to 
account for performance differences between similar airplanes. It must 
be concluded, therefore, that more effort should be spent on dealing with 
this phenomenon. Some simple experiments should be performed to support 
the theory. More attention should be paid to the roughness elements on 
airplanes which originate near spoilers and near flap and aileron hinges. 
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Figure 1.- Influence of disturbance on flow in boundary layer. 



Figure 2.- Velocity distributions of airfoil E603. 


Profile 603 with a 



FX 62-K-I53 



Figure 4.- Experimental envelope for FX 62-K-153/20 airfoil. 



Figure 5.- Velocity distributions for FX 62-K-153/20 airfoil; 

original coordinates. 



Figure 6.- Velocity distributions for FX 62-K-153/20 airfoil; 

smoothed coordinates. 
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GENERATION AND BREAKDOWN OF AERODYNAMIC LIFT: 


PHYSICAL MECHANISM 


Wolfgang Liebe 
Techn. University Berlin 


INTRODUCTION 

At high angles of attack the condition of attached flow is dependent upon 
the balance between forces on the leading edge and on the trailing edge of the 
wing. In the case of low-speed flight the wing operates at the upper limit 
imposed by this balance. In order to calculate and possibly influence this 
limit, a detailed understanding of the physical mechanism is required. We do 
know how to generate lift forces and we are able to calculate their magnitude 
as well as their distribution along the wing span. We do not know, however, 
the real physical mechanism of lift generation. 


THOMSON’S THEOREM 

The lift force is a result of circulation, i.e., a net flow around the air- 
foil. A differential pressure results with a corresponding force perpendicular 
to the main flow; this lift force is proportional to the airspeed and to the 
intensity of the circulatory flow. In order to discuss the problem of lift 
generation, attention has to be focussed on the circulation itself. 

We learn from standard literature (refs. 1 and 2) the way circulation is 
produced: A layer of separation arising at the trailing edge coils up to a 

starting vortex. According to Thomson's theorem (ref. 3) the circulation along 
a closed flow path situated in a homogeneous inviscid fluid remains constant 
with time. Consequently the formation of the starting vortex requires the 
generation of an opposite circulatory flow. Since the total circulation is 
zero to begin with, the magnitude of the opposite circulatory flow is such as 
to compensate the starting vortex. Hence the starting vortex gives rise to a 
superimposed additional velocity of the fluid particles in the vicinity of the 
wing surface. 

At this point one may have difficulties realizing the logical sequence of 
events. Of course there are no doubts about the validity of Thomson's theorem. 
The application of this theorem, however, does not explain the physical origin 
of additional forces acting on the fluid particles. Presently we need some 
kind of electrodynamic "far field effect" to explain this fluid dynamics 
problem. 
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FLOW AROUND THE TRAILING EDGE 


The fluid particles pass the trailing edge with extremely high and local- 
ized velocities, especially during the early phase of the motion. Figure 1 
illustrates the corresponding flow pattern and the resulting low pressure zone. 
Static pressure differences are produced in the vicinity of the trailing edge. 

As a result more and more of the flow close to the surface moves towards the 
low pressure region even against the main flow. The "depression zone" is filled 
up by spirally moving particles forming a vortex sink. With increasing vortex 
diameter the flow velocity around the trailing edge decreases. This in turn 
reduces the suction and the transport of material towards the sink will decay 
as shown in figure 2. Finally the vortex reaches a critical size which is 
characterized by zero suction and vanishing flow around the trailing edge. At 
this moment the vortex is free; it separates from the wing surface (fig. 3). 



Figure 1.- Suction head at the trailing edge. 



Figure 2.- Accumulation in a 
vortex sink. 



As soon as the starting vortex drifts away, the condition of attached flow 
is attained. Smooth flow is established all the way down to the trailing edge. 
A combined local suction and material transport mechanism has been able to 
initiate a circulation. "Far field effects" are no longer required. 


VORTEX ROLL 

Filling up the suction zone is characterized by the formation of a typical 
flow pattern which generates what may be called a vortex roll. The intermittent 
phenomena taking place at the wing surface can be simulated by a continuous 
source-sink mechanism as shown in figure 4. 
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source 



Figure 4.- Source-sink mechanism. 


The mass flow discharging from the source equals the one Q 2 P 

entering the sink at any instant of time. Continuity of the incompressible 
flow requires (dots denote partial differentiation with respect to the time) 

dQ]^ + dQ2 = Q-j^ dt + Q 2 dt = 0 (1) 

In addition to this the stimulated circulation F is directly proportional to 
the volume flow according to 


dQ = -S dr (2) 

The quantity S represents the active span of the wing. 

Combining equations (1) and (2) yields 

dr^ + dV2 = dt + r2 dt = 0 (3) 

which indicates that flow continuity and the proportionality (2) are reflecting 
the relevant physical phenomena at the wing; in particular, equation (3) 
satisfies Thomson’s theorem at any instant of time. 


Finally at a time T the formation of the vortex roll is terminated which 
is characterized by 



(4) 


Consequently a circulation around the wing has been built up, having the same 
magnitude but different orientation from the final vortex roll drifting away. 


LIFT 


The lift force at the wing can be directly calculated from the properties 
of the vortex roll. According to Kutta's theorem the fluid velocity V pro- 
duces a lift force F at the active wing span S of 

F = r p V S (5) 
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r represents the steady state circulation around the airfoil, which according 
to equation (4) is replaced by -T 2 . Following figure 5 the vortex-roll 
circulation To is 


r, = -TT D W (6) 

where the quantity D represents the final diameter of the vortex roll and W 
is the local fluid velocity passing the trailing edge. 

Therefore the lift force is 

F = TTDWpVS 



Mechanical similarity requires a simple proportionality between geometry 
and velocity ratios as shown in figure 5: 

(D/L) = m(U/W) = m tan a (7) 

Herein L is the wing chord and U is the transverse velocity component at the 
trailing edge according to the angle of attack a. The factor m is a dimen- 
sionless coefficient. 

With the relations, 

D = m L tan a 

W = cos a V 

it is possible to calculate the lift force 

2 

F = m TT L tan a V cos a p V S = m tt sin a p V S L 
or specializing for small angles of attack a (wing area S L = A) , 

F = 2 TT m a (p/2)V^A 
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one obtains 


For the lift coefficient C 


'A 



= 2 TT m a 


( 8 ) 


(p/2)v2a 


Experiments indicate that the dimensionless coefficient m is close to 1,0 so 
that the following approximation is justified; 


The condition of attached flow dominates just after separation of the 
starting vortex. This condition, however, is not stable since dissipation and 
other effects are disturbing the flow. Consequently, the fluid again passes 
the trailing edge causing suction, which is able to correct for the disturbance. 
This mechanism at the trailing edge continuously and effectively maintains the 
condition of attached flow. The sharper the edge is, the more effective the 
mechanism is. 

At high angles of attack, however, a counteracting effect is initiated at 
the leading edge. 


Usually the leading edge is rounded; nevertheless similar processes take 
place as at the trailing edge. High fluid velocities are accompanied by strong 
suction. Local backflow is initiated in the boundary layer, but no vortex roll 
can be formed in the front. This is not due to rounding of the leading edge but 
to the fact that the front depression zone has no direct contact with a region 
of significantly higher static pressure which is able to fill a vortex roll. 

With increasing angle of attack there is an Increasing static pressure gradient 
from the front to the back of the upper wing surface. But the boundary layer 
is able to resist a major backflow. So in spite of this pressure gradient, the 
front low pressure zone remains isolated. 

As far as this isolation is concerned, there is a fundamental difference 
between the processes at the leading edge and those at the trailing edge. In the 
case of the back depression this zone is being rapidly filled from a very close 
reservoir, the high pressure stagnation zone. In the case of the front region 
the low pressure is continuously maintained since there is no high static 
pressure reservoir available which could form and fill a vortex roll. 

At high angles of attack the pressure gradient along the upper surface 
rises considerably; this changes the situation drastically: The thickness of 

the boundary layer increases, giving rise to a backflow at the wing surface. 

This reverse wedge flow expands from the back and reaches the front depression 
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THE TRAILING-EDGE MECHANISM 


BREAKDOWN 
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zone. At the moment of contact a rather violent inflow takes place towards 
the centre of the zone. A flow unbalance results since the depression zone 
takes in more material than the main flow is able to deliver. As a conse- 
quence the low pressure region is filled very fast, which in turn leads to 
a rapidly growing vortex roll. Finally the flow separates as indicated in 
figure 6. 


These phenomena now correspond directly to what has previously been 
described as happening at the trailing edge, except for one important differ- 
ence: There is a net flow around the airfoil which reduces the circulation and 

causes lift breakdown. 

It has been shown that flow separation is not only initiated by static 
pressure rise and friction; a third condition has to be satisfied: The condi- 

tion of contact between the low pressure zone with regimes at a higher static 
pressure. 


The range of steady lift generation could be extended if one were able to 
prevent contact of flow from the trailing edge with the front depression. One 
device for this purpose is shown in figure 7. 


o 


{> 



Figure 6.- Above: Expanding reverse flow. 

Below: Contact - the flow separates. 


COUNTERMEASURES 



Figure 7.- Pockets at the upper surface. 
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Such "pockets" at the upper wing surface are well known from the wings of 
birds. At high angles of attack the backflow in the boundary layer causes these 
elastic elements to raise. This way the expansion of the backflow is effec- 
tively stopped; contact is prevented and lift breakdown is delayed. 

The outlined principle is not yet in use in aeronautics, but it has proven 
successful in biotechnics (fig. 8) . So far only one case is known where such a 
pocket-type device has been tested on a stalling airplane — with good results. 



Figure 8.- Heron during approach for landing. 

(G. Rueppell, Vogelflug, Kindler 1975) 

Up to now only two-dimensional flow conditions were considered. In real 
situations lift breakdown mostly starts locally somewhere along the wing span. 
Localized lift breakdown, however, being limited to a short part of the span 
may lead to a three-dimensional flow. Now the low pressure zone in the vicinity 
of the local breakdown fills up from separated regions causing sideways inflow 
of material. At high angles of attack that sideways influx rapidly propagates 
to the wing tip. Like a chain reaction the lift collapses all of a sudden 
along the whole wing as shown in figure 9. 

As a countermeasure an effective device has been suggested: The boundary 

layer fence. A simple shroud is mounted on the wing in order to protect the 
outer part of the wing against infiltration. Thin threads have been fastened 
to the wing surface to make the flow pattern visible during flight (fig. 9) . 

The fence does prevent sideways contact and subsequent lift collapse. 
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Figure 9.- Flow pattern on a stalling airplane. 

Right: Sideways influx rapidly propagating 

to the wing tip. Left: A fence protects the 

outer part against infiltration. 

Quite frequently the boundary layer fence is used in conjunction with 
sweptback wings. Sweepback means shifted airfoils which is accompanied by 
sideways pressure gradients. Thus the low pressure zone is filled sooner from 
the adjacent side, causing not only reduced lift but also unfavorable pitching 
moments. These inherent effects can be prevented by fences. 


WING FLAPPING 

The outlined extremely high velocities around the trailing edge can be 
provoked by a transverse motion of this edge relative to the main flow, for 
example by moving a trailing edge periodically up and down. Corresponding 
vortex rolls are generated and material is absorbed from the boundary layer, 
thus reducing its thickness. The reverse wedge flow decreases and again lift 
breakdovm is delayed. 

The volume flow due to wing flapping can be calculated from the following 
relation: 


Qf 


= 2 TT 


3 ,,2 p2 r^S L^ 


(9) 
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where 


volume flow due to flapping 
k dimensionless coefficient 
n flapping frequency 

Pq angular amplitude of flapping 
r radius of moving edge 

S span of active wing 

L wing chord 

V fluid velocity 

The material taken out from the boundary layer is accelerated backwards by 
flapping, which as a reaction produces wing propulsion. The thrust attained in 
this manner is negligibly small in gases; in high density fluids it is possible 
to produce significant propulsive forces by this technique. 


SUMMARY 

A contribution has been given to an old problem: The explanation of the 

generation of aerodynamic lift. New physical models are described which provide 
a better understanding of the phenomena involved. The suggested viewpoint 
leads to new technological implications. The formation of both a starting vor- 
tex and a circulation can be conceived as the filling of a vortex sink at the 
trailing edge. Fluid is absorbed by the vortex, which causes it to expand to 
a vortex roll. The lift force can be calculated from the properties of the 
vortex roll. 

Once the starting vortex drifts away, the condition of attached flow is 
attained. With increasing angle of attack this condition is disturbed by low 
pressure close to the leading edge. Finally this depression zone fills from 
the back of the wing, which induces a countercirculation and lift breakdown. 

Filling requires the low pressure region to have contact with flow regimes 
of higher static pressure. Flow separation caused by filling of the vortex at 
the leading edge can be influenced by anti-contact devices such as pockets or 
fences. 

A periodic flow can be superimposed around the wing by a forced oscillatory 
motion of the trailing edge. The periodic formation and separation of small 
vortex rolls reduce the drag or even produce propulsion. 
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INTRODUCTION TO THE ARCOPTER ARC WING 


AND THE BERTELSEN EFFECT FOR POSITIVE PITCH STABILITY AND CONTROL 

William D. Bertelsen 
Bertelsen, Inc. 

SUMMARY 


Studies in the realm of low-speed and motorless flight have 
traditionally produced the most creative approaches to the prob- 
lem of flight. The problem is the same today as always, namely, 
the search for higher performance with complete safety. Towards 
that end a brief report is offered on a new wing design, new in 
geometry, construction, and flight characteristics. This report 
includes preliminary wind tunnel data on a three-dimensional mod- 
el as well as some full-scale man-carrying test results. There 
are photos of all phases of the experiments and some figures which 
serve to illustrate the Bertelsen Effect, a unique focus of aero- 
dynamic forces in the arc wing system which allows the attainment 
of high lift coefficients with the maintenance of pitch stability 
and control. 


INTRODUCTION 


The name "Arcopter" comes from a combination of the Latin 
word "arc" for segment of a circle with the Greek "pteron" for 
wing. The name thus embodies the basic geometric configuration of 
the device. In this case the arc refers not to any chordwise air- 
foil curvature but to a regular spanwise curvature describing an 
arc like a rainbow over the lateral pitching axis of the system. 
From antiquity the arch has been an element of structural design 
and it has come to be a symbol of strength and simplicity. This 
paper introduces the arc wing configuration as a novel aeronauti- 
cal device which embodies certain valuable aerodynamic properties 
in a light-weight, self-constituted physical unit of inherent 
strength and simplicity. 


ABBREVIATIONS AND SYMBOLS 


Values are given in both SI and U.S. Customary Units. The 
measurements and calculations were made in U.S. Customary Units. 


103 


VTOL 

m. p . h . 
R.P.M. 


C.L. 

C.G. 


L 

D 

T 

W 

a 

X ' 

z 

c 



vertical takeoff and landing 

miles per hour 

revolutions per minute 

center of lift 

center of gravity 

aerodynamic center 

lift, kg (lbs. ) 

drag, kg (lbs.) 

thrust, kg (lbs.) 

weight, kg (lbs.) 

angle of attack, degrees 

lift moment arm, m (ft.) 

drag moment arm, m (ft.) 

wing chord length, m (ft.) 

velocity, m/sec (ft/sec) 

coefficient of lift 

coefficient of drag 

moment coefficient about a.c. 

moment coefficient about C.G. 


INVENTION OF THE ARC WING 


The VTOL Design Problem 

In the 1950‘s the National Advisory Committee for Aeronautics 
(NACA) recognized the usefulness of the helicopter because of its 
ability to operate from very small bases. The advantages to be 
gained with an airplane that incorporated both the small-field ca- 
pabilities of the helicopter and the high-speed potential of con- 
ventional airplanes became readily apparent (ref. 1 ). One possi- 
ble means of achieving these advantages was seen to be an engine/ 
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propeller combination capable of providing static thrust in excess 
of gross weight. Lift for vertical takeoff could then be obtained 
by deflecting the propeller slipstream downward by means of large- 
chord wing flaps, retractable for high-speed cruising flight. 
Accordingly, an investigation of various wlng/flap configurations 
was conducted in the 7- by 10-foot tunnels at the Langley Aero- 
nautical Laboratory in an effort to develop relatively simple 
arrangements that could deflect propeller slipstreams downward 
for vertical takeoff. References 1, 2, and 3 present the charac- 
teristics of slotted, sliding, and plain flaps, respectively. 

The slotted-flap configuration was effective in achieving a 
slipstream turning angle corresponding to a rotation of the effec- 
tive thrust vector upward about 73 °, with the ratio of resultant 
force to thrust varying from about 1.00 nearest the ground to 
about 0.86 out of the ground effect region. With this configura- 
tion it was concluded Jhat vertical takeoff could be made with an 
initial attitude of 17 ° and at airplane weights up to 90 percent 
of the total propeller thrust. 

Similar results were achieved with the plain flap configura- 
tion, but only after the installation of auxiliary vanes which 
greatly complicated the arrangement. The slotted-flap configura- 
tion, while seen as somewhat simpler, had the disadvantage of ex- 
hibiting rather large diving moments, caused partly by the fact 
that as the flaps extended they moved appreciably rearward and the 
effective axis of the redirected slipstream was relatively far be- 
hind the quarter-chord point of the wing. For the same turning 
angle the diving moments associated with the slotted-flap config- 
urations were found to be approximately twice as large as the div- 
ing moments for the configurations with plain flaps and two auxil- 
iary vanes. However, the process of retracting and storing the 
two auxiliary vanes necessary on the plain flap system was seen to 
present serious mechanical problems, nearly prohibitive to the de- 
sign of a practical, high-speed VTOL aircraft. 

Subsequent investigation of the Ryan VZ-3RY VTOL prototype 
(ref. 4) under the auspices of the new National Aeronautics and 
Space Administration in 1959 underscored the serious limitations 
of the conventional approach to the design of double-slotted flaps 
for VTOL applications. While the aircraft could take off verti- 
cally, longitudinal stability was said to be Impossible to realize 
below 46 km/hr (29 m.p.h.) with the existing center of gravity lo- 
I cation. Pitch control in hover and transition was difficult and 
critical even with a complicated jet-reaction control located in 
the tail. 

In summary, the experiments established (1) that VTOL capa- 
bilities are possible with slipstream deflection by means of con- 
ventional, large-chord, double-slotted flap arrangements and (2) 
that, because of the large diving moments associated with extended 
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double-slotted flaps, a VTOL using such an arrangement may be 
longitudinally trimmed and controlled In either the hovering 
mode or the transitional mode but not In both, at least not by 
simple means. A VTOL aircraft should be stable and controllable 
in hover, transition, and high-speed cruise. Conventional tail 
surfaces are totally ineffective at zero forward speed and almost 
ineffective in transition. The NACA-NASA studies of the 1950' s 
indicated that an unconventional approach would be required to 
meet the VTOL design challenge. 


The Arcopter VTOL 

The Arcopter wing system was the direct result of the ef- 
forts of Dr. William R. Bertelsen to develop a slotted flap con- 
figuration for deflecting a propeller slipstream through the 
large turning angles required for vertical takeoff without the 
deleterious diving moments or complexity which accompanied the 
NACA experiments. The cited NACA technical notes touch on the 
importance of the center of gravity location in analysis of VTOL 
wing and flap pitching characteristics. 

In the Arcopter system it is proposed that if an aircraft 
extends single or multi-element flaps and/or slats for high lift 
or slipstream deflection, then those flaps should rotate, while 
extending, about an area In which the center of gravity lies, so 
that the summation of flap resultant forces converges at all 
times in the vicinity of the center of gravity. Each flap may be 
considered in such regard as an entity with its own force focus 
coinciding with the others near the center of gravity. The cen- 
ter of gravity is preferably below the center of lift of the air- 
foil combination in order to effect stability regardless of the 
attitude of the aircraft with respect to gravity. The concentra- 
tion of wing forces, coupled with flap and wing slot augmentation, 
all converging about the center of gravity of the aircraft, thus 
comprises an engineering principle called the Bertelsen Effect. 

Figure 1 shows how multi-element flaps might be arranged to 
take advantage of this principle in a VTOL of the deflected slip- 
stream type. It can be seen that if flaps B and C retract and 
extend by pivoting on a radius centered at the C.L./C.G. focal 
point there will be little or no diving moment at any flap set- 
ting. The aircraft can therefore make the transition from hover 
to flaps-up cruise smoothly and predictably. Because of the fo- 
cus of flap resultant force through the C.G. area, the system does 
not depend on propeller thrust to achieve longitudinal trim in 
any mode. Forces remain balanced at all power settings including 
power off, affording an extra measure of safety in controlling a 
power-off descent. The full lifting capability of the wing sys- 
tem can be utilized in all modes without the usual loss in effec- 
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tive lift coefficient owing to negative tail loads. The need for 
such negative loads is effectively eliminated in the Arcopter sys- 
tem. 


Figure 1 diagrams the general arrangement of wing and flap 
elements around the center of gravity. It is left to specify the 
most practical physical form to be taken by an aircraft which is 
to employ the Bertelsen Effect. It has been established that mul- 
ti-element slotted flaps of large chord can deflect a propeller 
slipstream through the large turning angles required for vertical 
takeoff (ref. 1). Reference 5 suggests the effectiveness of large 
end plates in augmenting flap efficiency as regards the ratio of 
resultant force to thrust, especially in ground effect. The Arc- 
opter system proposes a synthesis of the wing and end plates into 
a spanwise, semicircular arc as being the most efficient configu- 
ration for confining and deflecting the slipstream of one large- 
diameter propeller or two dual-rotating propellers on a single 
thrust axis. At the same time the necessary rotational motion of 
arc-shape flap elements can be easily achieved, owing to the con- 
venient coincidence of element pivot points on an axis across the 
diameter of the arc. This location and coincidence of wing element 
pivots substantially simplifies the mechanism for flap actuation. 

Figures 2-6 are photos of the Arcopter VTOL flying model 
which was built to demonstrate the Arcopter design principle and 
the Bertelsen Effect. Figure 2 shows the arc wing and flap ele- 
ments fully extended. Such arched structure is inherently strong 
while being light in weight. Because of the great tensional 
strength of the arch structure, there is no longer a design re- 
quirement for thickness in the structure of the main wing. Air- 
foils can be chosen without regard for structural considerations . 
The wing and flaps on the model are constructed of molded Plexi- 
glas sheet. Aluminum tubes attached to the model are for handling 
and serve no aerodynamic function. Figure 3 shows the complete 
VTOL model in a three-quarter front view, flaps fully extended. 
Simple canard control vanes have been included in the slipstream 
to counteract propeller torque and provide positive three-axis 
control at all speeds including zero and reverse. It can be seen 
that any residual diving moment can be dealt with by Increasing 
Incidence on the horizontal canards In such a way as to contribute 
to the overall slipstream-turning and lift effectiveness of the 
system. 

Figure 4 is a direct front view showing shortness of the 
wingspan. If a single thrust axis is to be used on an Arcopter 
VTOL, the wingspan should be somewhat less than the diameter of 
the prop or rotor. Short span saves weight and reduces drag in 
high-speed cruise. The tubular diametric spar visible in figure 4 
is oversized for rough handling. Figure 5 shows a side view of 
the Arcopter VTOL in high-speed cruise conf iguration with flaps 
retracted. Thin flap segments easily nest in the main arc wing 


107 


after simple rotational motion, pivoting about the diametric axis 
through the wing tips at the point where the tightening nut at- 
taches the handle. A glow plug engine drives the propeller via an 
extended shaft to help maintain proper C.G. location. Viewed from 
the top (fig. 6) the arc wing is seen to have an elliptical plan- 
form, and therefore a near-ideal lift distribution. 

The Arcopter VTOL model was demonstrated (In and out of ground 
effect) in hover, in slow flight fore and aft, for control effec- 
tiveness, etc., in the 7- by 10-Foot Tunnel at Langley Aeronautical 
Laboratory on January 23, 1958. As a solution to the VTOL design 
problem, the Arcopter offers a simpler, safer alternative to the 
helicopter through implementation of the Bertelsen Effect. But the 
invention of the arc wing as an element in itself, an arch-tension 
structure with centralized force focus, offers possible solutions 
to a variety of aeronautical design problems, especially those 
where lightness of weight and structural simplicity are prime 
considerations . 


THE ULTRA-LIGHT ARCOPTER WING 


As indicated previously, the arc wing may be considered in 
multi-element combinations , as in the VTOL discussion, or each el- 
ement may be considered separately as an entity with its own cen- 
tralized force focus. Figure 7 represents the Bertelsen Effect as 
it applies to a single-element arc wing. It shows how lift force 
acts in a direction perpendicular to imaginary lines tangent to 
each point along the semicircular arc wing span. The magnitude of 
a local lift force through the point of tangency is proportional 
to the local wing chord length and angle of attack. On an arc 
wing with an elliptical planform, the greatest lift will develop 
near the crown of the arch where the wing chord length is great- 
est. The vector L represents the relative magnitude and direction 
of lift force acting on this point with respect to the lift forces 
which act simultaneously on every other point along the span. Be- 
cause the arc wing is a semicircle as viewed from the front, it 
becomes clear that all lift forces, regardless of magnitude, aim 
through a common point at the geometric center of the arc. This 
point is the true center of lift in the Arcopter system. 

On the right in figure 7 is a side view of the arc wing fo- 
cus. This side view shows how lift and drag forces interact at 
each local section center of pressure to focus a resultant force 
directly through the geometric center of the wing arc. If the 
aircraft C.G. is also located near this point, the vector sum of 
the forces is zero, and there is no pitching moment about the C.G. 
The longer broken lines denote the outline of the arc wing lead- 
ing and trailing edges as seen from the side. Vector R' has the 
same magnitude and direction as the resultant R and acts through 


108 


the same point. R' is simply a restatement of the resultant R for 
convenience in graphically adding R to the weight W and thrust T. 


I Structure of the Ultra-Light Arc Wing 

It was decided to design and build a single-element arc wing 
to analyze its aerodynamic properties, including lift, drag, and 
static pitch stability. The basic simplicity and tension strength 
of the arc geometry implied that an ultra-light structure could be 
devised which could support a very large wing area. Intuitively, 
the semicircular shape is suited for confining high pressure air 
underneath the wing surface by effectively restricting spanwise 
flow. If the wing were properly designed, the lift force created 
by the free stream should stretch a single- or double-surface fab- 
ric membrane into an efficient airfoil curve without the necessity 
for any rib structure whatsoever, at a great saving of weight, 
cost, and complexity. 

The ultra-light, adjustable-camber arc wing evolved during 
numerous experiments with models, kites, wind tunnel tests, full- 
scale force tests, and finally, man-carrying, powered, free-fllght 
tests. The wing Is essentially a fabric tension structure utilizing 
the dynamic force of the air to stretch the wing fabric on the bias,- 
thus maintaining a single-surface airfoil curve. (See figure 8.) 

An aluminum tube forms the basic arch inside the fabric cuff at 
the wing's leading edge. This aluminum (or fiber glass) arch is 
anchored in sockets at opposite ends of a rigid tubular spar. 

This arch and spar assembly forms a "D" shape unit which has prov- 
en extremely rugged and damage resistant on test craft of every 
size. The sail is patterned after the elliptical planform of the 
Arcopter VTOL flap elements, with maximum chord length at the 
crown of the arch. 

Maximum chord length was specified arbitrarily to be one- 
third the length of the wingspan for all size test aircraft, fix- 
ing the aspect ratio at about 3*9 to 1. The wing fabric itself is 
non-porous urethane-coated nylon weighing 88 g/rn (2.6 oz/yd^). 

There is no continuous rigid structure shaping the wing sail ex- 
cept for the arch tube in the leading edge. At zero airspeed the 
fabric droops limply from the arch. The only other members re- 
quired for proper shape in flight are a number of rigid tubes or 
sticks which extend between the leading and trailing edges of the 
wing at various stations on the span. The length of each of these 
chordwise members is adjustable, making it possible to change air- 
foil camber between flights. Shortening the tube increases the 
camber. The tubes are all double hinged at their leading-edge 
point of attachment to allow the sail to hang down at zero forward 
speed. Nylon webbing straps are sewn to the leading and trailing 
edge of each wing tip to transfer flight loads to the spar. The 
trailing edge webbing straps also serve an important pitch control 
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function to be discussed later. 


The preceding description of the ultra-light arc wing struc- 
ture is brief but complete. It is a supremely simple structure 
with few parts, but its arc conf iguration and adjustable-camber 
surface enable it to develop respectable lift coefficients. At 
flying speed all waviness and wrinkles disappear as the fabric 
stretches to its cambered airfoil contour without the use of ribs 
or battens. The natural load distribution of the arc configura- 
tion seems to prevent fluttering of the trailing edge without the 
need for battens. Also contributing to efficiency is the elimina- 
tion of the usual fuselage junction losses which disturb most wing 
mid-sections. The device shown in figure 8 can be built to almost 
any size without complicating the design. The arc wing photo- 
graphed in figures 8-10 has a wingspan of 3 m (10 ft.). This unit 
was used extensively to develop structural design and fabrication 
techniques, as well as to study pitch stability and control in 
tethered flight. 


Pitch Stability and Control 

As apparent in figures 8-10 the Arcopter wing has inherent 
positive static pitch stability in flight without the addition of 
auxiliary stabilizing surfaces which most aircraft require. This 
stability is largely independent of airfoil section characteris- 
tics. Any airfoil section can be employed on an arc wing accord- 
ing to performance requirements. Moreover, the angle of attack at 
which the wing stabilizes can be completely controlled by varying 
the tension in the trailing edge of the wing. This is easily ac- 
complished by tightening or loosening the nylon webbing strap 
which anchors the wing fabric to the spar at the trailing edge. 

Figure 11 shows the full-size Arcopter wing built to carry a 
man. Wingspan of this unit is 7-3 (24 ft.). Clearly visible at 

the wingtip trailing edge is a steel cable attached to the nylon 
anchor strap. When the cable is connected to a trim tab crank or 
control stick, the pilot can control the wing's pitch attitude In 
flight. Pulling on the cable causes the wing to stabilize at a 
higher angle of attack. Releasing tension causes the wing to 
pitch down to a more shallow angle of attack. Recovery from a 
completely luffed condition resulting from negative angles of at- 
tack can be made at once by pulling on the cable. Continuing to 
draw the trailing edge down results in stabilization at extremely 
high angles of attack, upwards of 40°. Experience has shown that 
at high angles of attack the arc wing behaves like a parachute and 
cannot be stalled in the normal sense. Releasing some tension on 
the trailing edge produces immediate wing response, restabilizing 
it at some lower angle of attack. Total cable travel required for 
the whole flight range is only about 15 cm (6 in.). 


no 


Center of Gravity Location 


As in the case of the Arcopter VTOL, the pitch stability and 
control behavior of the ultra-light arc wing is primarily related 
to the location of the center of gravity with respect to the vec- 
tor sum of all aerodynamic forces acting on the wing. Figure 12 
is a representation of wing and low C.G. location which is some- 
what like the arc wing situation. Values can be assigned to the 
lift moment arm x' and the drag moment arm z for each angle of at- 
tack to be considered. The value of x' is taken to be negative 
when the C.G. lies ahead of the aerodynamic center (a.c.). The 
value of z is taken to be positive when the C.G. is below the a.c. 
The formula in figure 12 is developed in reference 6 to express 
the pitching moment about the C.G. when the C.G. location and rel- 
ative forces are known. Positive values of indicate tendency 

to pitch up and negative values Indicate tendency to pitch down. 

Using the formula, a family of curves of versus angle of 

attack can be developed to predict the basic pitch stability char- 
acteristics of a conventional wing with a C.G. located 1.5 chord 
lengths below the section a.c. Figure 13 is such a plot using co- 
efficient values of an NACA 23012 wing of aspect ratio 6. The an- 
gle of attack corresponding to a pitching moment of zero is called 
the "trim point". The slope of the curve at the trim point is an 
indication of the static pitch stability of the system; the more 
negative the slope, the more statically stable the wing. Aft 
movement of the C.G. results in a trend toward increased stability 
at higher angles of attack. 

Figure 13 is a hypothetical case not meant to represent the 
exact behavior of an arc wing, but it does indicate the large in- 
fluence C.G. location has on static pitch stability. Minor ad- 
justments in C.G. location might be made accordingly which would 
enable an arc wing to use any airfoil section and yet retain a ze- 
ro pitching moment about the C.G. at the design lift coefficient. 
When the C.G. is fixed, minor shifts in the focal position of the 
force vectors (from changes in trailing edge tension) give total 
pitch control on the ultra-light arc wing. 


Preliminary Wind Tunnel Tests 

Through the cooperation of the late Dr. H.S. Stillwell, then 
head of the University of Illinois Department of Aeronautical and 
Astronautical Engineering, a brief series of tests were conducted 
on a single-element arc wing of ultra-light construction in the 
university's 1.5-m by 1.5-m (5 ft. by 5 ft.) low-speed wind tunnel 
in 1973 and 1975- The 1975 data is included in figures 14 and I 5 . 
Figure I 6 shows the model installed inverted in the test facility. 
Wingspan was 1.2 m (4 ft.). The model was of the same construc- 
tion in all respects as described previously, including the wing 
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fabric of non-porous urethane nylon, sewn to allow some bias 
stretch . 

Several problems combined to interfere with the accuracy of 
the test results. In the first place, the model was perhaps too 
large for the tunnel. Secondly, the smooth airfoil camber which 
characterizes all the larger arc wings failed to develop on the 
small wind tunnel model. Thirdly, during the course of the test- 
ing, the tunnel screens were discovered to be dirty, thus creating 
extra turbulence. The screens were removed for the tests labeled 
"MAX" and MED" camber in figures 14 and 15 . This raised the tun- 
nel Reynolds number based on mean geometric chord from about 

0.24 X 10° to 0.32 X 10°. Nevertheless data was taken and tabula- 
ted for values of C^, and for three varying degrees of 

airfoil camber, the extremes of which can be seen in figures I 7 
and 18 . Because of inability of the fabric to stretch naturally 
into airfoil camber on the small model, it had to be induced by 
bending the three most central chordwise tubes. 

The lift and drag measurements indicated disappointing per- 1 
formance by the model compared with expectations based on experi- j 
ence with the large arc wings in the field. However the maximum 
value of did show increase with increasing camber as might be | 

expected. One interesting result was that the wing never reached « 
the stalling point in any of the tests. Limitations of tunnel j 

balance apparatus precluded investigation of very high angles of I 
attack, but it can be seen in figure l4 that the wing with maximum j 
camber did not stall even after a 29 ° increase in attack angle, , 

beginning at of about 0.23* 

The pitching moment data (fig. 15) taken about the horizontal 
spar shows a negative slope, indicating a degree of positive stat- 
ic pitch stability, in all three tests. Increasing airfoil camber 
appears to produce greater positive (nose-up) values of the C.G. 
moment coefficient at low angles of attack. Positive moments re- 
main near the maximum even in the vicinity of the zero-lift angle. 
Unfortunately , angles of attack below the zero-lift angle were not 
investigated. The pitching moment data implies that the arc wing 
will retain a measure of positive static stability about the C.G. 
no matter what airfoil curvature is employed. Increased camber 
seems to have a favorable effect on static stability. 


Piloted Tests of the Full-Scale Arc Wing I 

I 

By 1976 the 7*3 m (24 ft.) span Arcopter wing was ready for 
limited flight testing with a pilot aboard. The wing itself, as 

2 2 

shown previously in figure 11, of projected area 13-9 m ( 15 O ft ) ,, 
was fitted to a heavy-duty tricycle landing gear for auto towing. 
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In addition, a unique annular rudder-elevator provided yaw control 
and contributed to pitch control (fig. 19) • Like the wing, the 
"ring tail" is a light-weight tension structure with a circular 
rigid hoop inside the fabric cuff at the leading edge. In the 
same manner as the wing, the tail design provides a maximum of ef- 
fective area with a minimum of structure. 

The auto-tow tests of the Arcopter "sailplane" were very lim- 
ited, intended only to gauge the minimum flying speed at gross 
weight. On one such experiment observed by Dr. Stillwell, the 
wing lifted a total of I 70 kg (375 lbs.) at 10.7 m/sec (24 m.p.h.). 
No drag measurements could be made, but the low-speed lifting po- 
tential of the full-scale wing was substantiated. Some pilotless 
tethered flying was also conducted in moderate wind of 8.9 m/sec 
(20 m.p.h.), as shown In figure 20. Empty weight was 107 kg 
(235 lbs.). These tethered flights indicated the wing to be so 
stable in pitch at moderate to high angles of attack that the ring 
tail could not effect any visible change in pitch attitude. Sub- 
sequently all pitch control was accomplished by regulating trail- 
ing edge tension according to the method described earlier. 

In an attempt to come as close as possible to full-scale 
flight conditions for the purpose of measuring drag on the Arcop- 
ter sailplane, a trailerized mobile force balance was constructed 
in 1977 . The complete airframe, including pilot, was mounted on 
an articulated steel pylon via a ball-and-socket joint at the air- 
craft C.G. (figs. 21, 22 ). The airframe was thus free to pivot 
about three axes, making it possible to check out control systems 
as well as to monitor and record airspeed and drag values from 
calibrated pressure gauges connected to small hydraulic cylinders. 
Cylinder pressure, being a function of the total drag, was contin- 
uously recorded on movie film, as was airspeed from a boom-mounted 
pitot tube. It was also intended to measure lift with the trailer 
apparatus but the lift balance failed to function in the predicted 
manner . 

Good drag data was obtained by towing the rig on smooth 
blacktop. The aircraft was set to stabilize at an angle of attack 
of 13 - 15 ° so a plot could be made of drag versus airspeed. Two 
days of testing produced the data presented in figure 23. On a 
number of test runs the arc wing itself was removed from the rest 
of the airframe in order to measure and compare drag on pilot and 
supporting structure alone. The resulting figures could then be 
subtracted from the total drag to gain a more meaningful idea of 
the drag on the wing as a separate entity. The drag figures ob- 
tained on the Arcopter sailplane were low enough to suggest that a 
very small engine would be sufficient to propel the aircraft and 
pilot in flight, without necessitating an increase in wingspan. 

Besides facilitating drag measurement, the mobile force bal- 
ance made it possible to safely observe the behavior of the 7.3 m 
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arc wing at speed. The ball-and-socket coupling at the C.G. al- 
lowed the entire aircraft adequate freedom to pitch, roll, and 
yaw. Test runs were made with and without the ring tail at vari- 
ous speeds in an attempt to ascertain general handling qualities 
and control responses. The following conclusions were drawn con- 
cerning stability: 1. Without the ring tail, the arc wing has 

only neutral static yaw stability. 2. The arc wing has positive 
static pitch stability over a wide angle of attack range, with or 
without the ring tail. 3- The arc wing is neutral in roll sta- 
bility, but gets increasingly positive as the C.G. is lowered be- 
low the center of lift focus. There is no damping in roll. To 
effect roll control, the arch structure was hinged on wingtip 
"toggles" to enable the pilot to shift the entire wing and center 
of lift to the right and left relative to the C.G., but response 
was sluggish and inconsistent . Pilot weight shift did produce a 
slow but sure response without adverse yaw. 


The Powered Ultra-Light Arcopter B-IA 

The experiments conducted with the Arcopter sailplane were 
important, but certainly not exhaustive. The relative merits of 
the differing degrees of camber available in the adjustable-camber 
wing were not explored. But it was proven that the 7«3~m arc wing 
can carry significant pay loads at low speed. Some evidence was 
obtained also that power requirements for takeoff are low even 
with a fairly short wingspan. Piloted flights, powered by a small 
engine, now more certainly establish the efficiency of the ultra- 
light arc wing as a lifting device. 

The powered Arcopter B-IA was built using the same size wing 
and tail surface but with a simpler structure supporting the pilot 
in a prone position on a steel cable anchored at the wingtips. 
While being lighter and more streamlined, this design has the ad- 
ded advantage of being collapsible for car-top transport. Before 
installation of the power system and landing gear, the new air- 
frame made piloted tethered flights as a foot-launched hang glid- 
er (fig. 24). Empty weight of this configuration, including pilot 
harness and ring tail, is only 39 • 9 kg (88 lbs.). Flights out of 
ground effect were made with a wind of 8.9 m/sec (20 m.p.h.) sus- 
taining a gross weight of 104 kg (230 lbs.), indicating a lift co- 
efficignt of 1.5* Angle of attack at the maximum chord station 
was 13 ^ with respect to the horizon as measured from other photos. 
The attainment of lift coefficients near 2.0 seems a reasonable 
possibility, given a modest increase in angle of attack. 

The Arcopter B-IA (figs. 25-28) has now made successful short 
flights, out of ground effect, under its own power with a pilot 
aboard. Gross weight was 147 kg (325 lbs.). The nucleus of the 

power system is a small 2-cycle engine of 134 cm^ (8.2 in^) dis- 
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placement which drives two opposite-rotating pusher propellers, 
each 1.07 m (42 in.) in diameter (fig. 27). A maximum of 45 kg 
(100 lbs.) static thrust is available at 65 OO engine shaft R.P.M. 
for takeoff and climb. Minimum takeoff airspeed is about 10 m/sec 

( 22.3 m.p.h.). 


Future Experiments 

Flight testing of the Arcopter B-IA has only just begun, and 
experiments will continue. More investigation is warranted be- 
cause the first flights of the B-IA show that the Arcopter wing 
conf iguration offers a maximum of performance from a minimum of 
structure, with possible aeronautical utility ranging from ultra- 
light sport flying to high-speed VTOL transportation . 
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Figure 2.- Arcopter VTOL model, 
flaps fully extended. 



Figure 3-- Arcopter VTOL model, 
with canard control vanes . 
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Figure 4.- Arcopter VTOL model (front view). 
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Figure 5*- Arcopter VTOL model, flaps 
retracted for cruise. 



Figure 6.- Arcopter VTOL model, elliptical 
planform wing (top view). 
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Figure 7.- Arc wing forces, single-element focus. 



Figure 8.- Ultra-light arc wing, 
3-m (10 ft.) wingspan. 
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Figure 9-- 3-in arc wing, 
pitch-stable flight. 
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Figure 10.- 3-ni arc wing, tethered 
flight with payload. 



Figure 11.- 7 . 3-m (24 ft.) arc wing, 
pitch control detail. 




Figure 12.- Moment relation of wing and low C.G. 
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Figure 13 •- Pitching moment of wing with low C.G. 



ANGLE OF ATTACK a DEGREES 
Figure l4.- Arc wing lift and drag (wind tunnel data). 
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Figure 15- - Arc wing pitching moment 
(wind tunnel data). 



Figure l6.- 1.2-m (4 ft.) arc wing 
in wind tunnel. 
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Figure 17 •- Wind tunnel test, 
minimum camber. 
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Figure l8.- Wind tunnel test, 
maximum camber. 



Figure 19- - Arcopter 7 • 3-m sailplane, 
ready for tow. 
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Figure 20.- Arcopter 7.3-m sailplane, 
tethered flight. 
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Figure 21.- Arcopter sailplane on 
mobile force balance. 
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Figure 22.- Arcopter and pilot on 
mobile force balance, detail. 
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Figure 23 ,- Total drag of arcopter sailplane 
(mobile balance data). 



Figure 24.- Arcopter foot-launched 
sailplane, tethered flight. 
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Figure 25.- Arcopter B-IA powered 
ultra-light aircraft. 



Figure 26.- Arcopter B-IA, empty weight 
82.5 kg (l8l lbs. ). 
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Figure 27.- Arcopter B-IA, 


power system detail. 



Figure 28.- Arcopter B-IA In powered flight. 
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SOME NEW AIRFOILS 


Richard Eppler 
Universitat Stuttgart 
Stuttgart, West Germany 


SUMMARY 


A computer approach to the design and analysis of airfoils and some common 
problems concerning laminar separation bubbles at different lift coefficients 
are discussed briefly. Examples of application to ultralight airplanes, canards, 
and sailplanes with flaps are given. 


INTRODUCTION 


In the 1940' s, NACA demonstrated clearly that it is possible to design 
airfoils from pressure distributions in such a way that the boundary layer 
would behave in a desired manner (Refs. 1 and 2). At that time, it was dis- 
covered that the boundary layer would remain laminar longer if the pressure 
minimum occurred further aft on the airfoil. This realization led to the first 
laminar airfoils. Since that time, better methods for designing airfoils from 
pressure distributions have been developed (Ref. 3). Simple methods for com- 
puting the characteristics of laminar and turbulent boundary layers including 
a feasible transition criterion have also been developed (Ref. 4). The 
occurrence of laminar separation bubbles has been detected and studied experi- 
mentally (Ref. 5) and correlated with theory (Ref. 4). Good methods for the 
analysis of the potential flow around a given airfoil have been developed 
(Ref. 6). Thus, it was possible to write computer programs which combine all 
of these methods. These programs allow airfoils to be designed with prescribed 
pressure-distribution properties, the boundary- layer characteristics to be 
determined, and the effects of shape modifications such as plain or variable 
geometry flap deflections to be analyzed. A complete description of such a 
program system will soon be published as a NASA technical memorandum (Eppler and 
Somers). This system is somewhat equivalent to a wind tunnel. Three funda- 
mental differences do exist, however. First, the computer analysis of an air- 
foil is much less expensive than the corresponding wind-tunnel test. Second, 
the total time required to obtain the final results is much shorter. Third, 
much more data, such as development of the boundary-layer shape factor and 
thickness, are available. Moreover, the modification of an airfoil through 
prescribing the pressure distribution, which must be done on the computer, is 
integrated into the program system. This allows a boundary-layer development 
with prescribed properties to be obtained directly. 

Thus, the time has come to use the computer when a new airfoil is to be 
developed. Wind-tunnel and flight tests should be used to obtain a better 
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understanding of fundamental phenomena in support of the theory. Accordingly, 
an appropriate, or even an optimized, airfoil could be developed for each 
application rather than looking for an acceptable airfoil in an airfoil cata- 
log. All such catalogs together could not cover all practical requirements. 

The Reynolds numbers, wing loadings, flaps, takeoff and landing requirements, 
structural constraints, moment restrictions, surface qualities, and many other 
specifications vary over wide ranges. It is not possible to develop catalogs 
for all such requirements. Only for a few applications, such as sailplanes 
with smooth surfaces and model airplanes, have catalogs been used successfully 
(Refs. 7 and 8). Even for these applications, new requirements arise which 
cannot be satisfied by existing airfoils. Other applications (e.g., general 
aviation, remotely piloted vehicles, and hydrofoil boats) are still far from 
having a list of standard requirements. 

So, the tailoring of airfoils to specific applications becomes increasing- 
ly important. This paper presents some general considerations for tailoring 
airfoils and some examples of specific applications. 


GENERAL CONSIDERATIONS 


Airfoil design means to specify an airfoil from its pressure distribution 
in such a way that the boundary layer behaves in a desired manner. This 
approach usually leads to certain problems. Some of these problems are briefly 
discussed in this section. 

The velocity distribution over an airfoil changes with angle of attack. 

An example is given in Figure 1 which shows the velocity distributions of an 
airfoil at seven angles of attack. (Note that all velocity distributions in 
this paper are presented in terms of the ratio (V) of the local potential-flow 
velocity to the free-stream potential-flow velocity.) The differences between 
the different curves are nearly independent of the particular airfoil and are 
approximately proportional to the differences between the corresponding flat- 
plate velocity distributions. Normally the design of an airfoil means the 
specification of the entire velocity distribution at one angle of attack only. 
This is called a one-point design. The design method mentioned previously 
(Ref. 3), however, permits a multipoint design in which the velocities are 
specified along different segments of the airfoil at different angles of attack. 

For Reynolds numbers below about 4 x 10^, one of the most important problems 
concerns laminar separation bubbles which usually occur if transition takes 
place in an adverse pressure gradient. It is well known that this phenomenon 
can cause a substantial increase in the total drag (Ref. 5). This increase 
depends primarily on the Reynolds number R and the degree of adverse pressure 
gradient near transition. At lower Reynolds numbers, less adverse pressure 
gradient is allowed. A so-called "transition ramp" must be introduced ahead of 
the pressure recovery in order to obtain awfully developed, turbulent boundary 
layer. At Reynolds numbers below about 10 , a fully developed, turbulent 
boundary layer is not possible at all and, accordingly, the adverse pressure 
gradient can be only slightly steeper than the one which a laminar boundary 
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layer could overcome without separating. The theory (Ref. 4) as used in the 
program system provides a certain bubble analog. If this analog is prevented, 
the real flow does not normally show an additional bubble drag. 

The problems associated with laminar separation bubbles become more diffi- 
cult as angle of attack changes. As shown in Figure 1, the transition ramp 
introduced on the upper surface at high angles of attack a is reduced and 
even eliminated at lower a. For all multipoint designs, this problem is most 
difficult to solve. Fortunately, another effect helps the situation. For an 
airplane in flight, the Reynolds number changes with angle of attack or lift 
coefficient c^ . Thus, lower c means higher velocity and correspondingly 
higher Reynolds number. This fact can be exploited by requiring a less steep 
transition ramp at lower c^. On the upper surface, it is even possible to 
eliminate the transition ramp required at higher c and, thereby, allow an 
extension of the laminar flow region at lower c^ and higher P,. On the lower 
surface, a laminar separation bubble and even separation of the turbulent 
boundary layer can be permitted at low c and low R. As R increases to 
the free-f light value, the bubble and the turbulent separation should disappear. 
As c^ increases, the adverse pressure gradient should be reduced to an amount 
suitable for a transition ramp. 

All of these features are illustrated in Figure 2 which contains the theo- 
retical section characteristics for the airfoil shown in Figure 1. This air- 
foil was designed for a sailplane. The Reynolds number corresponding to low 
c is approximately R = 3 x 10^. The Reynolds number for high c is about 
R^= 10 . For c < 0.5 and R = 10^, which is not achievable in flight by the 
sailplane, turbulent boundary-layer separation was permitted on the lower sur- 
face. As c^ is decreased from 1.2 to 0.6, the transition point on the upper 
surface moves aft approximately 10% of the chord because the transition ramp 
essentially "disappears." 

Some unpublished wind-tunnel data (Althaus, Universitat Stuttgart, 1975), 
and free-f light data (Ref. 9) are included in Figure 2. The latter data agree 
very well with the theory, while the wind-tunnel results show some discre- 
pancies. The differences in transition point are inconsequential because a 
microphone was used in the wind tunnel to detect transition. This technique 
probably detects only a fully developed, turbulent boundary layer, and there- 
fore, experimental points lay somewhat behind the theoretical ones. Of more 
importance are the differences among the drag polars. The wind-tunnel curve 
for R = 10^ is characteristic of a polar for an airfoil with a small laminar 
separation bubble. That is to say that low drag is achieved at low and high, 
but not medium, lift coefficients. This problem was apparently not experienced 
in flight. Even more important are the drag differences for c^ < 0.2 and 
R = 3 X 10^. Here the free-flight tests indicate that the theoretical results 
are probably more reliable than those measured in the wind tunnel. 

In summary, it is very likely that the "computer wind tunnel" can predict 
at least the differences between different airfoils so reliably that it should 
be used to design an airfoil for a specific application. 
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AIRFOILS FOR ULTRALIGHT AIRPLANES 


Ultralight airplanes usually have only one side of the airfoil covered. 

This means that the airfoil has essentially zero thickness. The structure is 
concentrated primarily near the leading edge and to a lesser extent near the 
trailing edge. The problem, then, is the sharp suction peak which occurs near 
the leading edge at all off-design conditions. A high maximum lift coefficient 
c and a soft stall are desirable for takeoff and landing, whereas because 

^max 

of trie low aspect ratio, the lift coefficient for minimum sinking speed as well 
as for maximum glide ratio is usually somewhat less than c . Good penetra- 

Itld.X 

tion at even lower c^ is also sometimes desired. Thus, the problem is to 
design thin airfoils exhibiting a range of lift coefficient over which the flow 
is not entirely separated. Some thickness is, of course, required near the 
leading edge for structure. The following examples demonstrate what can be 
achieved by carefully shaping the leading-edge region. The first example, 
airfoil 379, is shown in Figure 3 along with its velocity distributions. At 
a = 7° relative to the zero-lift direction, a very high suction peak has 
already occurred on the lower surface near the leading edge. On the upper sur- 
face, a suction peak forms as a increases but the AV /Aa is much less than 

max 

for the lower surface. The pressure recovery is slightly concave, but by no 
means as severe as the recovery typical of the Stratford distribution. This 
shallow, concave pressure recovery together with the rounded, upper-surface 
suction peak results in a soft stall which is most important for the application. 

The section characteristics for this airfoil are shown in Figure 4. A high 
maximum lift and a soft stall are achieved, but below c = 1.0, the lower-sur- 
face flow is separated. The separation is predicted at about x/c = 0.8. This 
is a consequence of the assumption that the flow will reattach in a favorable 
gradient which, in this case, is probably not true. Thus, the flow on the lower 
surface must be considered separated from the leading edge aft. 

An attempt to lower the lower-surface, leading-edge suction peak is shown 
in Figure 5. This airfoil, 378, is much thicker than the previous one (3.88% 
versus 2.10%). As shown in Figure 6, lower-surface separation is now predicted 
below c^ = 0.6, and thus, a much wider range of lift coefficient is available. 

Figures 7 and 8 show airfoil 377, which is similar to 378 except that it is 
shifted to a higher lift coefficient.. Using the design method mentioned above, 
this is easily accomplished. 

The lower surface of this airfoil was then modified so that "zero" thick- 
ness was reached at a more forward x/c. The new shape and its velocity distri- 
bution are shown in Figure 9 and an overlay of Figures 7 and 9 is presented 
in Figure 10. Notice that the lower-surface flow exhibits much more adverse 
pressure gradient after the modification. As a consequence, the flow on the 
lower surface for this case is separated at all lift coefficients. This demon- 
strates the danger involved in arbitrarily modifying an airfoil to a shape which 
only looks appropriate. 
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Airfoil 376 was designed to have the same upper-surface behavior as air- 
foil 377 but to have less thickness and reach zero thickness at about x/c = 0.25 
(Fig. 11). This airfoil has a maximum thickness of 2.21%. It has a certain 
c range over which the flow is not separated, and hence, is much better than 
airfoil 377 modified (Fig. 12). This range is still considerably less than that 
for the original airfoil 377. 

These five examples illustrate the possibilities for thin airfoils. Many 
other constraints probably exist and, therefore, more tailoring would be required 
for this application. 

Another category of ultralight airplanes is becoming more popular, the so- 
called foot-launched sailplane with an empty weight of around 45 kg, full con- 
trols, and an enclosed cockpit. This concept was demonstrated in the 1930 's 
when the "Windspiel" was built. Today's materials allow much more efficient 
structures than were available at that time. 

The airfoil requirements for this application include high maximum lift 
coefficient, soft stall, and low drag down to c ~ 0. Because of the low wing 
loadings involved, penetration always means low c . Airfoil 748 (Fig. 13) was 
tailored for this application which covers a Reynolds number range from 0.6 x 
106 to 3 X 106 (Fig. 14). This airfoil requires a smooth surface for the for- 
ward 45% of the chord. If this can be accomplished, an aircraft with much 
lower wing loading than, say, a Ka-6 or Schweizer 1-26 can achieve the penetra- 
tion of these heavier sailplanes and yet have a minimum speed which would per- 
mit simple takeoff procedures including foot-launch from a ridge with little 
wind. 


AIRFOILS FOR CANARDS 


Because of longitudinal-stability requirements, a canard (forward wing) 
must always operate at a higher c than the main (rear) wing. The maximum 

lift coefficient of the main wing is, therefore, constrained by the c of 

^max 

the canard. Thus, it would be senseless to incorporate lift-increasing devices 
on the main wing if none were included on the canard. Fortunately, the canard 
usually includes an elevator which is deflected down to obtain higher c from 
the main wing. Thus, the elevator acts as a lift-increasing device for the 
canard. This effect, however, does depend on center-of-gravity position. 

The design objectives of airfoils for canards, therefore, include high c-^ 

with small downward flap deflection, low drag at low c^ with no flap defTec- 
tion, and a certain thickness for structural reasons. The Reynolds numbers are 
relatively low because of the small chord lengths. 


Two examples illustrate this application. The velocity distributions for 
the first example, airfoil 1230, are shown in Figure 15. The upper surface is 
designed only for high c^ • This is accomplished by preventing suction 

peaks and by including a certain transition ramp. Even at low c^, only 20% of 
the upper surface can sustain laminar flow. The lower surface can have about 
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50% laminar flow. The theoretical section characteristics are shown in Figure 
16. For positive flap deflection (down), some problems exist at low c^. This 
combination, however, cannot occur in flight. The second example, airfoil 1233 

(Fig. 17), achieves even higher c^ (Fig 18). This airfoil is also thicker, 

nidx 

and therefore, a drag penalty is paid at low c^. The lower surface of this air- 
foil can sustain only 30% laminar flow. An airfoil between these two examples 
has been successfully applied on Burt Rutan's "Defiant" (Ref. 10). 


AIRFOILS FOR SAILPLANES WITH FLAPS 


Sailplanes with normally hinged flaps are a standard application of air- 
foils. The difficulties with this application come from two requirements. 

First, the flap-down case usually corresponds to a Reynolds number of 10^ or 
below. For this case, laminar separation bubbles can be dangerous. This danger 
is increased by the steep adverse pressure gradient immediately downstream of 
the suction peak at the flap hinge. Second, the negative-flap-deflection (up) 
case corresponds to R > 3 x 10^. For this case, transition can occur earlier 
than desired. For a zero pressure gradient at these Reynolds numbers, the 
boundary layer is not stable enough to remain laminar for 60% to 70% of the 
surface and, therefore, a certain favorable pressure gradient is necessary to 
keep the boundary layer laminar. 

Airfoil 662 was designed for this application. The velocity distributions 
for this airfoil with flap deflections (3) of 0 , 10 (down), and -7.5 (up) are 
shown in Figure 19. The pressure recovery on the upper surface for the unde- 
flected-flap case must be less than would be possible for the case where no 
flap deflections were intended. A flap deflection in either direction increases 
the amount of adverse pressure gradient. Severe separation would occur in these 
cases if the pressure recovery for the undeflected case were already approaching 
the separation limit. The flap deflection can, however, be exploited in a 
favorable sense as well. For the flap-down case, a distinct transition ramp 
forms between the original pressure recovery and the suction peak caused by the 
flap. On the lower surface, an additional favorable pressure gradient occurs with 
the flap up which stabilizes the laminar boundary layer at the higher Reynolds 
numbers. Attention to all of these details together with the careful designing 
of the leading-edge region results in the good performance illustrated in Fig- 
ure 20. Notice that, at low c^ and low R, a lower-surface separation was 
again permitted. 

Another application resulted from the practical achievement of the variable- 
geometry concept. A flap which extends the chord 20% while introducing essen- 
tially no disturbances in the flap-retracted configuration was developed by 
F. Mahrer and incorporated into his sailplane, "Delphin" (Ref. 11). This flap 
could only be applied over that portion of the span which required no aileron. 

It was, therefore, desirable to deflect the ailerons down for the high-lift 
case. A negative flap deflection was not allowed. Thus, an airfoil was required 
which would have a laminar bucket that would extend down to around c^ = 0.05 
and which would achieve a high c with a plain and a variable-geometry flap. 

'max 
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The velocity distributions for such an airfoil, 664, are shown in Figure 21. 
The transition ramp between the original pressure recovery and the flap hinge 
is again exploited for the flap-down case. The favorable pressure gradient aft 
of x/c = 0.5, however, had to be introduced for this airfoil because no flap-up 
deflection was possible. The section characteristics for this airfoil are shown 
in Figure 22. 


CONCLUSIONS 


Some new airfoils have been designed for specific applications through the 
use of a computer program. The applications included ultralight airplanes, 
canards, and sailplanes with flaps. The coordinates, moment coefficients, and 
zero-lift angles for all the airfoils presented are given as an appendix. The 
tailoring of airfoils should be encouraged because it is highly unlikely that 
airfoil catalogs will be produced for all possible applications. The relia- 
bility of this theoretical approach increases as more wind-tunnel and flight- 
test data are correlated with the theory. So far, many such theoretically 
developed airfoils have been successfully applied. 
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APPENDIX 


COORDINATES, MOMENT COEFFICIENTS, 

AND ZERO-LIFT ANGLES FOR VARIOUS AIRFOILS 
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a relative to zero-lift direction. 



Figure 2.- Section characteristics for airfoil 603. 
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Figure 3.- Velocity distributions for airfoil 379. 



Figure 4.- Theoretical section characteristics for airfoil 379. 
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Figure 5.- Velocity distributions for airfoil 378. 



Figure 6.- Theoretical section characteristics for airfoil 378. 
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Figure 7.- Velocity distributions for airfoil 377. 



Figure 8.- Theoretical section characteristics for airfoil 377. 
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Figure 9.- Velocity distribution for airfoil 377 modified. 



Figure 10.- Comparison of original and modified airfoil 377. 
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Figure 12.- Theoretical section characteristics for airfoil 376. 
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Figure 13.- Velocity distributions for airfoil 748. 
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Figure 19.- Velocity distributions for airfoil 662. 



Figure 20.- Theoretical section characteristics for airfoil 662. 
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Figure 21.- Velocity distributions for variable geometry airfoil 664. 



Figure 22.- Theoretical section characteristics 
for variable geometry airfoil 664. 
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A COMPARISON OF THE AERODYNAMIC CHARACTERISTICS OF EIGHT 


SAILWING AIRFOIL SECTIONS 

Mark D. Maughmer* 
Princeton University 


SUMMARY 


Because of its light weight, simple construction, and good aerodynamic 
performance, the Princeton sailwing may be a competitive alternative to con- 
ventional wings for many low-speed applications such as ultralight sailplanes, 
man-powered aircraft and high-performance hang gliders. The operational 
characteristics of the sailwing are discussed with some emphasis placed on the 
importance of the trailing-edge cable tension as it controls several aero- 
dynamic properties. The three-dimensional aerodynamic characteristics of 
eight different sailwing profile sections have been obtained from wind tunnel 
tests and the results compared to determine the magnitude of the aerodynamic 
penalties paid for various structural simplifications. For the sectional 
thickness ratios considered in this research, it is concluded that, while the 
basic double-membraned sailwing has exceptional aerodynamic performance, even 
superior for some applications to the conventional hardwing, any notable 
deviation from this configuration results in an unacceptably large performance 
penalty. 


INTRODUCTION 


While there is currently a great deal of interest in the use of flexible 
wings for use on hang gliders, man-powered aircraft, and ultralight sailplanes, 
the design and evaluation of these vehicles is complicated a great deal by the 
fact that very little data is available to aid in analyzing the aerodynamic 
characteristics of such aircraft. Although the data presented herein was 
motivated by a National Science Foundation sponsored research program directed 
toward optimizing a windmill utilizing the Princeton sailwing, reference 1, it 
is hoped that these data will be of some value to designers of flexible-winged 
aircraft. 

The Princeton sailwing, which has been under development since 1948, is a 
unique, semiflexlble wing intended to provide the practical ultimate in a 
light-weight, low-cost lifting surface suitable for a number of low-speed 
applications. Basically, the structural configuration of the sailwing consists 
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of a leading-edge spar with attached ribs which ideally form a rigid framework 
supporting a trailing-edge cable in tension. A non-porous, non-stretchable 
cloth membrane, usually dacron, is then wrapped around the leading-edge and 
attached to the trailing-edge forming the upper and lower sail surfaces. The 
purpose of the pre-tensioned trailing-edge is to impart a chordwise tension in 
the membrane to minimize the deflections caused by the aerodynamic loads. 
Originally, it was not thought that the performance of such a device would 
compare favorably to that of a conventional wing; however, as a result of 
numerous experimental investigations, including those of reference 2, it has 
been found that the aerodynamic efficiency of the sailwing can indeed approach 
that of a hard wing. 

Specifically, the data presented in this report were obtained as a result 
of a wind tunnel program which was undertaken and structured in such a manner 
as to ascertain the relative magnitude of the penalties associated with using 
a readily available circular cross-sectioned leading-edge as opposed to the 
D-section normally used in sailwing construction. In addition, the importance 
of the full double cloth membrane was explored by testing sailwing sections 
which did not utilize the lower membrane as well as several having only a 
partial lower membrane. In total, eight wings, identical in all respects 
except for the airfoil section utilized, were tested and compared. 

In addition to the experimental development the sailwing has undergone, 
it should be noted that it has also received considerable analytical treatment 
such as that discussed in references 3-5. 


SYMBOLS 

b Span 

c Geometric mean chord 

q Dynamic pressure 

t Sectional thickness ratio 

AR Aspect ratio 

Drag coefficient, Drag/qS 

C Lift coefficient, Lift/qS 

J-i 

Moment coefficient about the quarter-chord point of the geometric 
mean chord, M/qcS 

2 

Trailing edge cable tension coefficient, Tension/qb 
L/D Lift-to-drag ratio 

S Wing area 
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MODEL DESCRIPTION 


The tests of the eight different wing profiles, shown in figure 1, were 
performed in the Princeton University 1.2 m by 1.5 m force-balance wind tunnel. 
The wing planform utilized in this study is characterized by a span b of 
.96 m, a geometric mean chord c of .113 m, and a total area S of .108 
square meters. The planform aspect ratio AR = b^/bc is equal to 8.5. 

Relative to the length of the geometric mean chord, the sectional thickness 
ratio t is 11.5 percent. It was possible to adjust the tension of the 
trailing-edge cable in each of the models which, for the results discussed 
herein, was set at 42 N as well as 160 N. These settings yield a trailing-edge 
cable tension coefficient C-j. of .07 and .28 respectively. 

TEST CONDITIONS AND DATA REDUCTION 


Most of the data collected in the series of experiments included in this 
report were obtained with the tunnel speed adjusted to yield a dynamic 
pressure q of 622.4 N/m^. Although the corresponding Reynolds number based 
on the geometric mean chord is calculated to be approximately 230,000, because 
of the turbulence level in the tunnel, the aerodynamic data collected is con- 
sidered to be qualitatively representative of a Reynolds number on the order of 
750,000. In any case, the test results were obtained at a Reynolds number 
above the subcritical value for which laminar separation dominates and are 
therefore indicative of the operating conditions expected for somewhat higher 
Reynolds numbers . 

The mounting of one of the test models in the wind tunnel is shown in 
figure 2. By means of an electrically driven tail-jack, it was possible to 
adjust the wing angle-of-attack to any value between -12 and +24 degrees while 
the tunnel was in operation. Thus, force balance data for lift, drag, and 
pitching moment were obtained at each two-degree angle increment between these 
limits. The test data were then reduced to the standard coefficient form and 
plotted as a function of the wing angle-of-attack as referenced to the unloaded 
(no-wind) orientation of the geometric mean chord, figures 5-8. In addition, 
the performance of each wing is summarized in a plot of lift-to-drag ratio as a 
function of lift coefficient, figures 9-10, and lift coefficient as a function 
of drag coefficient, figures 11-12. 


CHARACTERISTICS OF THE SAILWING 


Because many of the properties of the Princeton sailwing are uniquely 
different from those of a conventional hard wing, it is appropriate to discuss 
its operation to better understand the overall aerodynamic characteristics. 

For example, when the sailwing is at rest (no-wind), the cloth membrane is held 
taut by the trailing-edge cable and is essentially, except for the leading- 
edge, a symmetrical section as the upper and lower surfaces experience the same 
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pressure, figure 2. As the wing experiences airloads in a net lifting 
orientation (wind-on) , the asymmetrical pressure distribution that is estab- 
lished between the upper and lower surfaces causes the membrane (or membranes) 
to displace away from the high pressure regions (underside) and move toward the 
low pressure regions (upperside) . Thus, when a section is at an angle of 
attack greater than that of zero-lift, figure 3, it assumes a positive camber 
distribution that fairs the membrane in smoothly with the shape of the airfoil 
leading edge. It should be noted that the actual shape of the sailwing 
section is a function of the wind velocity, the angle-of-attack, the no-wind 
airfoil shape, and the amount of tension in the trailing-edge cable. Thus, as 
the angle-of-attack is increased, the resulting increased pressure differential 
between the upper and lower surfaces causes the amount of camber in the section 
to increase. This situation not only causes the maximum value of wing 
efficiency, the lift-to-drag ratio, to occur at relatively high lift coeffi- 
cients, but also delays the impending stall. At this point, the importance of 
maintaining the desired trailing-edge cable tension should be noted. As might 
be expected, relaxing the cable tension allows a greater amount of camber to 
be established and therefore a higher maximum lift coefficient is obtainable; 
however, simultaneously, the maximum lift-to-drag ratio is decreased as well 
as the threshold of critical velocity at which detrimental sail luffing occurs. 
Thus, the amount of tension in the sailwing trailing-edge cable controls an 
important trade-off between the maximum lift coefficient and the maximum lift- 
to-drag ratio. It might be considered that as the cable tension becomes 
higher, the sailwing behavior is more like that of a rigid wing while as it 
becomes lower, the behavior approaches that of a high aspect ratio Rogallo 
wing. 


Another interesting characteristic of the lifting sailwing is the upward 
deformation of the trailing-edge in the unsupported mid-span regions of each 
wing-panel, figure 4. The result of this action is a reduced angle of attack 
in these regions and one would expect a local reduction in lift; however, it is 
generally the case that this effect is more than offset by the Increased 
amount of camber that occurs which results in a local increase in the genera- 
tion to lift. In fact, because of this effect, the span lift distributions 
that occur over many of the sailwings that have been tested are often very 
close to that of the elliptical optimum. 

The constant chordwise tension that is a result of the trailing-edge 
cable and the catenary-arc sail cut is responsible for many of the desirable 
features of the sailwing over other flexible designs. One such feature is that 
relatively low drags are present at low angles-of-attack and lift coefficients. 
Furthermore, unlike many flexible wing designs, the sailwing has the ability to 
pass smoothly through the zero-lift condition from that of a positive camber 
and lifting configuration to that of a negative camber and downward loading. 

All in all, through many years of extensive research, the sailwing has 
been found to provide a simple, light-weight, and low-cost alternative to the 
conventional hard wing while not suffering any notable performance penalties 
throughout many low-speed applications. 
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PRESENTATION AND DISCUSSION OF TEST RESULTS 


Representative three-dimensional lift, drag, and moment curves for four of 
the eight sailwing models tested are presented in the plots of figures 5-8. 

For each wing tested, the lift coefficient, drag coefficient, and pitching 
moment coefficient about the quarter-chord point, are plotted as a function of 
the angle-of-attack of the unloaded geometric mean chord. The lift-to-drag 
ratio as a function of lift coefficient for all of the sections tested are 
plotted in figures 9-10 while the drag polars, lift coefficient as a function 
of drag coefficient, are presented in figures 11-13. Finally, the typical 
effect of lowering the trailing-edge cable tension coefficient on the lift-to- 
drag ratio and the drag polar are shown in figures 14-15. 

It is important to note that a direct comparison of these data to those of 
a conventional wing is complicated a great deal by the flexible nature of the 
sailwing. For example, the sailwing data can be likened to that of a rigid 
wing in which an automatic flap is deflected an additional amount for each 
incremental increase in angle-of-attack. This characteristic is responsible 
for the fact that it is generally impossible in the case of a sailwing to 
linearize the drag polar or obtain a meaningful value for the span-efficiency 
factor as is done from wind tunnel test data for a conventional wing. 

Similarly, it should further be noted that at lower angles-of-attack (up to 
approximately five degrees), it is not uncommon for a sailwing to have a lift- 
curve slope which significantly exceeds the theoretical thin airfoil maximum 
for rigid wings of 2ir per radian. This occurs because the section is continu- 
ally varying camber over the angle-of-attack range. At higher angles-of-attack, 
the section is unable to deform proportionally as much as it does when less 
loaded and, therefore, as the angle-of-attack is increased to higher values, 
the lift-curve becomes increasingly more like that of a rigid wing. 

The most notable observation in comparing the compilation of data pre- 
sented in figures 9-12 is the significant performance advantage held by the 
conventional, double-membraned version tested over the more simplified ver- 
sions. This advantage is so great that it is inconceivable of a situation in 
which the potential benefit in weight saving, cost, or more-simplified con- 
struction for any of the modified versions could be justified in relation to 
the performance penalties. 

Examination of these data further indicate an important feature in that 
the sailwing highest L/D values occur in a range of relatively large values of 
lift coefficient. Furthermore, the flexible nature of the sailwing affords it 
low drag coefficient values over a relatively wide range of lift coefficients, 
particularly in the case of the double-membraned sections. 

In order to further emphasize the fact that unlike a conventional wing, 
the shape of the sailwing section is governed by the dynamic pressure (wind 
velocity) , a polar for sailwing model 2 is shown in figure 13 for the case of 
a constant wing loading. This polar was obtained by adjusting the tunnel 
velocity such that the lift force remained constant throughout the angle-of- 
attack range and is representative of the aerodynamic characteristics over the 
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speed range of an aircraft in level flight where the load factor is equal to 
unity. Thus, the high-speed flight conditions correspond to the lower lift 
coefficients while the low-speeds correspond to lift coefficients approaching 
the maximum value at stall. For comparison, the constant velocity polar of 
sailwing model 2 is also shown in figure 13. The difference between these two 
plots, excepting for the small contribution due to the changing Reynolds 
number over the speed range for the constant wing loading case, can be attri- 
buted to the flexible nature of the sailwing. For an equivalent hard wing, one 

would expect these two plots to be nearly identical. From the figure, it 

should be noted that the sailwing maximum lift-to-drag ratio of twenty, 
corresponding to the tangent to the curve drawn from the origin, occurs at a 

fairly high lift coefficient and that the operating range of low drag values is 

fairly wide. It is important to note, however, that in order to maintain a 
suitably high test Reynolds number these data were collected at a wing loading 
of 598 N/m^ which is relatively high for most motorless applications. More 
practical wing loadings would cause the maximum lift-to-drag ratio to corres- 
pond to slightly lower lift coefficient values. 

In attempting to generalize the effect of the different leading edge 
shapes employed, it is apparent that those sections having the smaller radius 
version have a wider region of low drag although, for the most part, the 
actual minimum values of the drag coefficient are slightly lower for the more 
rounded leading edge shape. Furthermore, as expected, the more pointed cross- 
section is accompanied by more abrupt stalling characteristics. Lastly, 
particular attention should be paid to the comparatively poor performance 
delivered by the often employed circular leading edge. 

The effects of lowering the trailing edge cable tension are as expected 
and summarized by the data plotted in figures 14 and 15. Briefly a reduction 
of the trailing edge cable tension leads to a higher value of the lift curve 
slope, gentler stalling characteristics, a higher value of the maximum lift 
coefficient, and a generally lower value of the maximum lift-to-drag ratio; 
however, it should be noted, as observed in figure 14, that the lower trailing 
edge cable tension results in larger lift-to-drag ratios occurring in the 
region near the maximum lift coefficient. 


CONCLUSIONS 


In reviewing the applicability of the data presented herein, it should be 
born in mind that the quantitative information becomes less valid as deviation 
occurs from the equivalent test Reynolds number; however, as the test condition 
is above that of the critical Reynolds number, the trends, relationships, and 
approximate values of the data should remain valid to somewhat higher values. 

In addition, note that the sectional thickness ratios employed for those tests 
are somewhat higher than those used on many current hang gliders and should be 
taken into account when contemplating the use of these test results. In 
particular, it is expected that the importance of the leading edge shape 
should diminish as the thickness ratio decreases and the circular leading edge 
shape should prove adequate, as it often has, for some applications. 
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In considering the findings of these tests in addition to many others, it 
has been demonstrated that the three-dimensional performance of the Princeton 
sailwing is quite competitive with many hard wings of the same aspect ratio. 
Thus, the use of the sailwing should allow the benefits of simpler construction 
and lower costs to be realized without paying any significant price in perform- 
ance. In fact, some consideration should be given to the fact that, unlike 
many of its rigid counterparts, the cambering characteristics found in the 
sailwing cause its three-dimensional lift-to-drag ratio to maximize at a 
fairly high lift coefficient as is desirable for many of the low-speed appli- 
cations for which the sailwing is suitable. Furthermore, relative to many 
conventional wing sections, the sailwing has the favorable characteristic of a 
fairly slow rise in the drag coefficient, and consequently a slow decrease in 
the lift-to-drag ratio, with increasing values of the lift coefficient. 

If one's mind is allowed to freely extrapolate from current trends and 
the results of this test program, one can envision the reality of a ultralight 
sailwing sailplane in which a cockpit lever is used to vary the wing trailing 
edge cable tension. Thus, in operation similar to that of a modern flapped 
sailplane, the pilot would slacken the cable upon entering a thermal to permit 
a slow tight circle with a high value of lift coefficient and upon exiting the 
thermal, pushes the lever to tighten the cable such that a high lift-to-drag 
ratio for inter-thermal cruise is obtained. Perhaps in the not too distant 
future, an aircraft of this type might bridge the gap between the limitations 
of hang gliding and the excessive costs of high performance. 
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Sailvane Model 1 Semi-Sailwing Model 1 Sailwing Model 1 



Sailvane Model 2 Semi-Sait wing Model 2 Sailwing Model 2 



Sailvane Model 3 Semi-Sailwing Model 3 

Figure 1.- Sailvane and sailwing sections tested. 
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Figure 2.- Typical sailwing model mounted in 
wind tunnel with wind off. 
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Figure 3.- Sailwing mounted in wind tunnel 
with wind on. 




Figure 4.- Detail of sailwing trailing edge deformations. 
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Angle of Attack: degrees 


Figure 5.- Aerodynamic characteristics of the sailvane model 2. 
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Figure 6 


Aerodynamic characteristics of the semi-sailwing model 
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Figure 7.- Aerodynamic characteristics of the semi-sailwing model 
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Figure 8.- Aerodynamic characteristics of the sailwing model 2 
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Figure 9.- Composite L/D vs. plot of sailwings with the larger radius leading-edge 
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Figure 10.- Composite L/D vs. plot of sailwings with the smaller radius leading-edge. 
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Figure 12.- 


Composite polar plot of sallwings with the 
smaller radius leading-edge. 
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Sailwing model 2 constant wing loading drag polar compared 
to the constant dynamic pressure drag polar. 
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Figure 13.- 
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Figure lA.- Effect of sailwing trailing-edge cable tension on 
lift-to-drag ratio: sailwing model 2. 
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Figure 15.- Effect of sailwing trailing-edge cable tension on the 
drag polar: sailwing model 2. 
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LENGTH AND BURSTING OF SEPARATION BUBBLES: A PHYSICAL 


INTERPRETATION* 

John M. Russell 

Department of Aeronautics and Astronautics 
Massachusetts Institute of Technology 


SUMMARY 


A physical interpretation of the observed form of the pressure distribu- 
tion beneath a two-dimensional "short" separation bubble (which modifies the 
external inviscid pressure distribution only locally) is given in terms of 
boundary layer concepts (i.e., constancy of static pressure across the layer as 
long as the layer is "thin"). At the mean separation and reattachment points 
(which lie on the same mean streamline) , the local static pressure equals the 
local stagnation pressure, since the velocity is zero at these points. The 
boundary layer hypothesis then implies that reattachment can only occur at a 
point X - Xg downstream of the separation point Xg if the jump in external 


inviscid static pressure 


-[u2(x) 


- tt2 I 


U^(Xg)] is less than or equal to the rise 


in stagnation pressure H(x) - H(Xg) along the separation streamline after 
separation. A simple method for estimating the growth of H(x) along a mean 
streamline entrained into the underside of a growing shear layer through the 
transition region is discussed, and predictions of bubble bursting conditions 
and a lower bound on the bubble length are compared with experiment. 


INTRODUCTION 

Since the early work of Melvill Jones (ref. 1), it has been known experi- 
mentally that a laminar boundary layer on an airfoil, after entry into a region 
of adverse pressure gradient strong enough to cause laminar separation, will 
separate in the laminar state, sometimes achieve transition to turbulent flow 
in the separated shear layer, and (if the adverse pressure gradient is not too 
severe) reattach to the surface to form a closed recirculating flow eddy known 
as a separation bubble. In a typical case (fig. 1), the turbulent boundary 
layer downstream of the reattachment point R either does not separate at all 
over the remaining portion of the upper surface or separates only a short dis- 
tance upstream of the trailing edge. In such a flow, the form of the external 
static pressure distribution outside of the region of separation is approximated 
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reasonably well by the potential flew pressure distribution, though the magni- 
tude of the actual pressure is somewhat smaller. It follows that the lift of 
an airfoil with this type of bubble separation is approximately equal to the 
value predicted by inviscid theory. Following Tani (ref. 2) , we will refer to 
this flow phencxnenon as a separation bubble of the "short" type, to distinguish 
it from another type of closed recirculating flow in which the entire pressure 
distribution differs radically frean the potential flow form. The latter is 
referred to as a "long" bubble. 

Separation bubbles are of great importance in engineering because of the 
role they play in the phenomenon of airfoil stalling, which may sometimes be 
identified with a transition of the flow frean a short type of bubble separation 
to the long type. This abrupt transition, known as "bursting" of the short bub- 
ble, is in seme cases responsible for producing the well-known "critical" Reyn- 
olds number region in the plot of drag coefficient versus Reynolds number for 
circular cylinders and spheres (ref. 2) and for moderately thick airfoils at 
Reynolds numbers lower than about 10^ (refs. 3 and 4). 

Reference 2 contains a thorough background discussion of the phenomenon of 
separaticai bubbles and a review of the experimental work prior to 1964. Refer- 
ence 5 reviews most of the more recent experimental work up to 1976. The reader 
is referred to these sources for general information about the subjects treated 
in this paper. 

In recent years, there have been several attempts to develop semi-empirical 
calculation methods for predicting the length and bursting conditions for sepa- 
ration bubbles of the short type. These range in sophistication from solution 
of the full time-dependent two-dimensional Navier-Stokes equations (with model- 
ing of the turbulence), as in reference 7, to simpler modifications of existing 
boundary layer calculation methods incorporating iterative methods to account 
for viscous-inviscid interaction between the external stream and the flow within 
the bubble (refs. 6 and 8). All of these methods use experimental data on sepa- 
ration bubbles to fix the values of certain numerical coefficients appearing 
in the theory. In each case, the shape of the separation streamline and the 
external pressure distribution between separation and reattachment are unknowns 
in the problem and are not determined until the calculation is finished. Since 
the distributions of these quantities play a crucial role in determining the 
length and bursting of the bubble, there is scane difficulty in using the numerical 
prediction methods (or the theories on which they are based) to deepen one's 
understanding of the fundamental physical processes controlling the flow. 

The objectives of the present study are to explain why reattachment occurs 
at all, to explain why the observed pressure distribution has its characteristic 
form, and to develop a shortcut method for calculating the length of a short 
separation bubble, if one exists, or to determine the conditions under which a 
short bubble cannot exist. 
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PHYSICAL MDDEL 


The Cause of Reattachment 

The principal cause of turbulent reattachment following laminar separation 
has been discussed by Cebici and Bradshaw (ref. 9) . Consider the separation of 
a laminar boundary layer frcm a smooth surface as shown in figure 2. The edge 
of the shear layer has been drawn to show the spreading of the shear layer in 
the downstream direction due initially to laminar and later to turbulent mixing 
of mcmentum. This spreading results in a greater mass flux across the seg- 
ment BB* in figure 2 than across the segment AA' . In two-dimensional flow, 
mass conservation requires that there be a net inflow across the segment AB. 
This entrainment of fluid into the underside of a growing shear layer will be 
greatly increased by the transition from laminar to turbulent flow. If this 
entrainment rate is greater than any reversed flow coming from far downstream, 
then the shear layer simply sucks itself back onto the surface from which it 
separated (ref. 9, p. 366) to form the short separation bubble shown in fig- 
ure 3, in which the entrainment flow is supplied by the splitting of the shear 
layer at reattachment. 


The Form of the Pressure Distribution 

In trying to identify the most important factors controlling the dynamics 
of separaticxi bubbles, it is useful to see how many of the qualitative features 
of the flew can be accounted for adequately by simple boundary layer concepts. 
If, for example, simple boundary layer theory is found to account adequately 
for the observed length and pressure distribution of a bubble, then one can 
infer that the direct effects of streamline curvature and nonzero pressure gra- 
dients in the cross-stream direction (which are neglected in boundary layer 
theory) are indeed unimportant, at least in determining the gross qualitative 
features of the flow. This process of testing hypotheses that certain terms in 
the equations of moticxi are negligible will serve to increase our understanding 
of virtiat a separation bubble is if it allows us to identify the few important 
frean among the many unimportant effects operating at once in the flow. 

In this spirit, we attempt to explain the characteristic pressure distri- 
bution beneath a short separation bubble (cf. fig. 1) by supposing that the 
assumptions of classical boundary layer theory hold. In particular, the static 
pressure across the separation bubble in figure 3 is assumed to be constant, and 
the velocity component in the direction normal to the wall is assumed to be an 
order of magnitude smaller than the component parallel to the wall. With these 
assumptions, the pressure distribution along the wall between points S and R 
in figure 3 can be inferred. The pressure gradient between points S and K 
must be small since the velocities and shear stresses in the slowly moving 
interior of the bubble are driven only by entrainment and are, therefore, small. 
In the region downstream of the point K, the streamwise velocity along each 
closed recirculating flow streamline will be zero at the downstream extremity 
of that streamline. Since the vertical velocity at such a downstream extremity 
is small, the static pressure is very nearly equal to the local stagnation pres- 
sure at that point. It follows that the distribution of pressure along the line 
connecting all such downstream extremities of streamlines (represented as the 
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dotted line in fig. 3) will approximately equal the distribution of stagnation 
pressure between the lower edge of the shear layer (at point K) and the sepa- 
ration streamline (at point R) . This difference in stagnation pressures may 
be substantial because the velocity on the separation streamline (which lies 
close to the center of the shear layer at the streamwise station K) will typ- 
ically be much larger than the velocity on the lower edge of the shear layer 
at any given streamwise station. 

In general, we expect the pressure distribution beneath the separation bub- 
ble in figure 3 to be as shown in figure 4, which is in good qualitative agree- 
ment with the examples of short separaticn bubbles in reference 2. The rather 
sharp corner in the pressure distribution in figure 4 is interpreted here as 
the point where the lower edge of the shear layer impinges on the wall. Accord- 
ing to our model, this point also coincides with the center of the innermost 
closed streamline in figure 3. Note that the point K is not identified with 
the transition from laminar to turbulent flow (which we assume takes place some- 
where upstream of the point K) . 

It has been assumed in all experimental studies on separation bubbles known 
to the author that the sharp corner in the pressure distribution coincides with 
the transition point. Gault (ref. 10) proposes this as a means of experiment- 
ally determining transition, while Van Ingen (ref. 11) uses the same assumption 
to fix the numerical values of the adjustable constants in a form of his "e to 
the ninth" method for transition prediction, developed especially for use in 
separation bubbles. In view of the fact that transition from laminar to turbu- 
lent flew occurs quite readily in a free shear layer in zero pressure gradient 
without modifying the static pressure, one is at a loss to see why transition 
as such should have any effect on the pressure distribution at all, much less 
account for the abrupt, nearly discontinuous pressure rise that occurs at the 
rear of a short separaticxi bubble. The rather arbitrary and ad hoc character 
of the conventional interpretation of the pressure distribution renders it less 
favorable from a fundamental viewpoint than the interpretation in terms of stag- 
nation pressures given above. 


Length and Bursting of Short Bubbles 

We have included among the defining features of short separation bubbles 
the condition that the bubble modify only locally the static pressure distribu- 
tion about the body on which the bubble occurs; that is, the pressure distribu- 
tion everywhere ahead of and behind the bubble is approximately equal to the 
value predicted by potential theory. In particular, the pressures at the sepa- 
raticn and reattachment points (which are the inner limits of the portion of the 
pressure distribution lying outside the bubble) must both be points lying very 
near the inviscid pressure distr ibuticxi curve. These endpoint conditions, 
together with our model of the pressure distribution within the bubble in terms 
of stagnation pressures on the recirculating flew streamlines, provide a clue 
as to what determines the length of a short separation bubble, when it is pos- 
sible for one to exist, and why it is sometimes impossible for a short bubble 
to exist. 
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Now the difference in static pressure between the points S and R in 
figure 3 equals the net rise in stagnation pressure along the separation stream- 
line, since the velocity is zero at both endpoints. It follows that if a plot 
of the rise in stagnation pressure along the separation streamline downstream 
of separation is drawn on the same graph as the rise in inviscid static pressure 
downstream of the same point, then an estimate of the bubble length is the 
streamwise distance from separation to the first point where the curves cross. 

One difficulty with such a graphical procedure for the calculation of the 
bubble length is that the curves may cross at more than one point, giving a 
choice of several possible bubble lengths rather than a single one. Another 
difficulty is that the stagnation pressure along the whole length of the sepa- 
ration streamline is not known under general conditions. It can, however, be 
estimated by assuming that the growth of the shear layer downstream of separa- 
tion is nearly equal to that of a free shear layer in zero pressure gradient, 
for which the mean velocity and, hence, stagnation pressure and stream function 
profiles are known from experiments on transition in a laminar mixing layer 
(ref. 12). 

A calculation based on similar physical concepts was made by Tani (ref. 2) 
to estimate the maximum attainable value of the coefficient of pressure recovery 
within a bubble cr defined by the relation 

Pr ■ PS 

a = (1) 

P 2 



Owen and Klanfer (ref. 13) had proposed the value Cfj^ax “ 0.35 from experi- 
mental evidence. By noting that the velocity along any streamline entrained 
into the lew velocity side of a growing shear layer cannot exceed the velocity 
on the dividing streamline (u/U)u_g^ = 0.5873 ..., Tani argued that a lower 
bound on the nondimensional stagnation pressure rise existed and was given by 

\ = [(u/U)d.sJ ^ = 0.3449 ... (2) 

max 



which is in excellent agreement with the figure 0.35 given by Owen and Klanfer. 
Tani did not suggest the extension of the stagnation pressure idea to the calcu- 
lation of the bubble length, however. 

The method of estimating the reattachment position as the streamwise coord- 
inate of the point of intersection of the curves of inviscid pressure rise and 
stagnation pressure rise on the separation streamline after separation provides 
a necessary condition for the existence of a short bubble, namely, that the two 
curves indeed intersect. A possible condition for the impossibility of a short 
separation bubble (i.e., a "bursting" condition) might then be that the inviscid 
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pressure gradient downstream of separation be so steep that the two curves never 
intersect. The borderline case in which a short bubble is just barely possible, 
but would become impossible with an infinitesimal increase in the inviscid pres- 
sure gradient, is when the two curves are tangent to each other at a single 
point. Later in this paper we will use this tangency condition to derive a 
close approximation to Caster's "bursting line" (ref. 14), which has been found 
by several experimentalists (see ref. 5) to be a reliable empirical formula for 
predicting the boundary between the short and long type of bubble separations. 


DEVELOPMENT OF THE SHEAR LAYER 

For reasons indicated in the last section, it is desirable to have an 
explicit formula for the development of the stagnation pressure H(s) along 
the separation streamline in a short bubble. 

For simplicity, we will assume that the actual distribution of stagnation 
pressure along the separation streamline can be approximated by the stagnation 
pressure distribution along a particular streamline entrained into the underside 
of a self-similar shear layer in zero pressure gradient, for which the velocity 
profile has been calculated from the Blasius equation in the laminar case 
(ref. 15), and has been measured experimentally in the transitional and fully 
turbulent cases (refs. 12 and 16, respectively). It has been found experimen- 
tally (cf. ref. 12, fig. 12) that the functional form of the mean velocity pro- 
file is experimentally indistinguishable for these three cases, the effect of 
transition being a large and abrupt increase of the rate of spreading of the 
shear layer measured by, say, the growth of the momentum thickness 6 (s) . 

Before treating the problems of the location of the transition region, the 
growth of the resulting turbulent layer, and the determination of the value of 
the stream function (relative to the reference streamline in a self-similar 
shear layer) appropriate for representing the separation streamline in a sepa- 
ration bubble, it is necessary to recall certain properties of self-similar 
shear layers, and, in particular, to obtain explicit approximate formulas for 
the velocity and stream function profiles. 


The Laminar Shear Layer 

The ncMnenclature for self-similar shear layers is as shown in figure 5, 
which illustrates the flow downstream of a hypothetical splitter plate, the 
upper surface of which moves with the velocity U-| (precluding the formation 
of a boundary layer there) , while the lower surface is stationary. The fluid 
above the plate (x < 0, y > 0) moves with the uniform velocity U-| while the 
fluid below the plate (x < 0, y < 0) is at rest. The effect of these boundary 
conditions is the elimination of any length scale other than the streamwise 
coordinate x, so that self-similarity in laminar flow is to be expected. 

The steady two-dimensional Prandtl boundary layer equations for zero pres- 
sure gradient. 
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The boundary conditions appropriate for a free shear layer are 
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The exact solution of equation (5) satisfying the boundary conditions (9) , 
(10), and (11) has been tabulated by Lock (ref. 15) to four significant figures. 
The tabulated solution provides a useful check on the accuracy of approximate 
representations of the velocity profile in terms of elementary functions. 
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Approximation of the Velocity Profile 


For the present purpose, it is sufficient to obtain an approximate solu- 
tion to equation (5) in terms of elementary functions whose forms are chosen 
to reproduce as many of the properties of the exact solution as is convenient. 
For example, the condition f(0) = 0 implies that f''(0) =0 (frcm equa- 
ticxi (5)), so that the velocity f has an inflection point at the origin, or, 
equivalently, the shear stress is maximum there. In addition, the conditions 
f(r|) * r| as n -*■ +oo and f(ri) « Constant as ri imply that f 

behaves like a complementary error function for large positive r| and like a 
simple exponential for large negative n. These conditions, together with the 
requirement that the velocity profile be a smooth monotonic function of ri, 
suggest the aj^roximate forms 


n > 0 : f(n) 



n < 0: f’ (n) 


f ' (0) - sech^ 
2 


bri + tanh~^ 



( 12 ) 


(13) 


each of which has an inflection point at r| = 0 and the correct asymptotic 
behavior as |t)| °°. Equations (12) and (13) are constructed so that the 

velocity f (T)) is continuous at the juncture n = 0. The requirement that 
the slope of the velocity be continuous at n = 0 provides one equation relat- 
ing the constants a and b: 





(14) 


Equation (14) provides an expression for the ratio a/b in terms of f'(0), 
where either a or b remains free to be determined by a suitable 
normalization . 

The only important unknown that remains is now f'(0), which we will deter- 
mine by a momentum integral method. Integrating equation (5) fran -®° to 0, we 
obtain 


2 f 


+ ff 



(f')2 dn 


0 


while the integral frcm 0 to °° gives 
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2 f " + f (f ' 

- 0 


"] CO ^00 

1) - f (-1 + f ') dn = 0 

- 0 ^0 


Applying the boundary conditions for a shear layer, these equations become 

2 f"(0) = r (f ) 

«^_oo 

which are forms of the Von Karman monentum integral equation (e.g., ref. 17, 
p. 146) for each half of the flow, and which have direct counterparts when the 
flow is turbulent (e.g., ref. 18, p. 227). 

If equations (12), (13), and (14) are substituted into equation (15), then 
a cubic equation for the quantity 1 - f ' (0) is obtained (see appendix) , the 
relevant solution of which gives f'(0) = 0.58923, which agrees well with the 
tabulated value f ' (0) = 0.5873 in reference 15. 

For our normalization condition (needed to determine either a or b from 
equation (14)), we note that equation (15) implies 

e 

- E 9 

z 

1 ° 

= f - f"(0) + f"(0) = -f(-°°) (16) 

J — oo 

which according to reference 16, has the numerical value -f (-°°) = 1.2386. 

It follows from equations (12), (13), (14), and (16) (see appendix) that 
a = 0.88544 and b = -0.30160, which completely determines all of the constants 
in our approximate velocity profile. For r) < 0 (which is the region of 
interest in the separation bubble problem) , the approximate velocity profile 
given by equation (13) is uniformly accurate to within 1 percent of the tabu- 
lated values in reference 16. 
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f ' (1 - f ') dn + 
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f ' (1 - f) dn 
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( 15 ) 


Development of the Momentum Thickness 

The nondimens ional velocity profile f ' (n) , together with the definition 
of the transverse length scale (x) in equation (6) and the free stream veloc- 
ity U-| completely determine the flow in a laminar shear layer. It is conve- 
nient to express S- (x) in terms of information available at a station x = x-j , 
rather than in terms of the distance x downstream of the splitter plate in 
figure 5. Frcxn equation (6) , 
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where A-] is defined to be S,(x^). Rewriting in terms of the momentum thick- 
ness 9(xi) and using the fact that 0(x)/8,(x) = 0* = 1.2386, we have 


0(x) 

0 (xi ) 


{QA) 


Ui 0(xi) 


■L V 


X - x-| 
0 (Xl ) 



(Laminar flow) (17) 


As mentioned before, we will assume that the functional form of the mean 
velocity profile in a self-similar turbulent shear layer is the same as in 
laminar flow. Thus, we continue to use equations (12) and (13) for the velocity 
profile, though the definition of the transverse length scale 2- (x) (which 
equals (Vx/U-|)^/2 in laminar flow) must be altered. Since the equation 
0=0* 5, (x) holds in general, this may serve as the definition of 5- (x) in 
turbulent flow, provided that an adequate empirical formula for the development 
of 0 (x) can be found. 


From the data of Sato (ref. 12), such a formula is not hard to develop. 
We may take, for example, in the transitional and fully turbulent regions of 
the flow 


0 (x) 


0 (x)‘ 
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and where [j9 /©si lam given by equation (17) with x-j = Xg. A plot 
of equation (19) with the numerical values X-i/6g = 34, ~ ^35, and 

Nj. = 56 is shown in figure 6. Also shown are the measured values of 
0(x)/0s - [0(x)/0g] ]_am for three different flows in reference 12. The agree- 
ment is seen to be quite reasonable. 

The condition (s - Xg)/0g = represented a transition criterion and 

appears to hold quite generally for two-dimensional laminar separated flows of 
the type under discussion. For example, figure 12 of Freymuth (ref. 19), which 
shows the growth of the streamwise fluctuation velocity component |u'|jj,ax 
a free shear layer downstream of separation for a variety of separation Reynolds 
numbers, indicates that the range of (x - Xg)0g in which the disturbances 
grow exponentially is almost completely independent of Reynolds number for 
61 i R(5^ i 334. This range was found to be 0 ^ (x - Xg)/6g S 60. Strong non- 
linearity (indicating the onset of transition) appears for (x - Xg)/0g ^ 60, 
which agrees well with the value = 56 used here. A similar independence 
of the onset of transition with Reynolds number can be found in separation 
bubbles. For example. Caster's results (ref. 14, fig. 11) of intermittency 
measurements in seven short separation bubbles in the range 136 S R^^ i 394 

show that the onset of transition occurs for 50 ^ (x - Xg)/0g S 80 for all but 
one of the bubbles. Unlike the case of fully attached flow, it therefore 
appears that transition in fully separated flew is relatively easy to predict 
in terms of the critical value of a single parameter, at least in the case of 
nominally two-dimensional incompressible flow in the Reynolds number range 
under discussion under conditions of small background turbulence and acoustic 
noise. 


VELOCITY DISTRIBUTION ALONG A STREAMLINE 


The stream function can be obtained from our approximate velocity pro- 
file (13) by integration. Applying the boundary condition = 0 on T) = 0, 
we have, from equations (7) and (13), 


vlJ = U-, 8-(x) f(B) 


Ui S<(x) 
b 


3 

2 


f’ (0) 



bq + tanh“^ 




( 20 ) 


Now the value of the stream function at q = -«>, which represents the mass 
flow across the entire lower half of the shear layer in figure 5, is a function 
of X and may be used to define the (constant) stream function along the 
streamline that just enters the underside of the shear layer at any station x. 
In our model of the separation bubble, we will suppose that the separation 
point x = Xg corresponds to the point where the separation streamline enters 
the underside of the separated shear layer in a bubble in just this way. That 
is, we will take 

Vljg = Ug Jl(Xg) f(-») = -Ug0g (21) 
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as the value of the stream function (relative to the shear layer) that repre- 
sents the separation streamline in a bubble. 


The distribution of velocity along the streamline can now be found 

easily. Fran equations (20) and (21) (with U] = Ug) , we have 


tanh 


bn 


+ tanh“^ 



^s 2b /e\ 1 

0 3 f (0)\£y 


( 22 ) 


It follows from equation (13) and the identity sech^ (x) 


tanh^ (x) that 




-2b6* 
3 f ' (0) 



(23) 


Figure 7 is a plot of equation (23) for various separation Reynolds numbers. 
Equations (17) through (19) have been used to evaluate the quantity 0/0g. 


STAGNATION PRESSURE ALONG Ip = il^g 

Equation (23) , which was calculated from the known behavior of a free shear 
layer in zero pressure gradient, provides a measure of the stagnation pressure 
distribution along the streamline il) = ipg, since in zero pressure gradient, the 
stagnation pressure H(x,i1j), defined by 


H(x,\l)) 



(24) 


is simply equal to the dynamic pressure (p/2)u^. 

We recall that in a general shear flow, the stagnation pressure varies 
along a mean streamline according to the approximate equation 


3h 

9s 


9t 

9n 


(25) 


where T is the combined viscous and turbulent shear stress 


9u 

T = p — - pu'v' 
9n 


(26) 
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and where s and n are orthogonal curvilinear coordinates in the plane of 
the flew oriented in the streamwise and cross-stream directions, respectively. 
Equation (23) thus represents the cumulative effect of such cross-stream stress 
gradients in the downstream direction. 

The use of equation (23) to model the actual stagnation pressure distribu- 
tion along i]; = in a separation bubble involves the assumption that the 
modification of the shear stress distribution due to the abrupt rise in static 
pressure at the rear of the bubble is small, or, equivalently, that the turbu- 
lent shear stresses are dominated by "memory" effects. Such an approximation 
of memory-dominated shear stresses is familiar in many analyses of turbulent 
flows following entry into a region of severe adverse pressure gradient (see 
especially refs. 20 and 18). The validity of this approximation when applied 
to the bubble problem has yet to be established, however. 


With due regard for this uncertainty, we may take the nondimens ional form 
of equation (24) as our working equation for the distribution of stagnation 
pressure along the streamline 4^ = ipg: 


H 


P 

2 




^ . ( 1 .]" 

P 2 \Us/ P 2 \Us /plane 

2 ® 2 ^ layer 


(27) 


where the last term on the right-hand side is given by equation (23) . 


CASTER'S BURSTING CRITERION 


M. Caster (ref. 14) 
parameters R0^ and the 


has obtained an empirical correlation in terms of the 
pressure gradient parameter P, defined by 


Ur - Ug 6g (Ur - Ug)/Ug Ug9 g 

P = = (28) 

^r — ^s ^ ^^r ~ ^ 


that distinguishes between flow conditions permitting a short separation bubble 
and the conditions when a short bubble is not possible (i.e., a "bursting" con- 
dition). Figure 8 shows Caster's so-called "bursting line," which we will 
attempt to derive from equation (27). 

If we define a pressure coefficient c^ relative to the free stream con- 
ditions at the separation point in a bubble, then 
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where Cp = 1 - (U/U^) ^ is the (X>nventional pressure cxsefficient. In terms 
of Cp, equation (27) becomes, on = ipg. 



layer 


At the reattachment point, the streamwise velocity u on 'V = 'I's 
so that 


(30) 


equals zero. 
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,U, 


plane 

layer 


(31) I 


where the pressure recovery coefficient a is the same as in equation (1). 

Equation (31) together with equation (28), which at reattachment can be 
written as 


(32) 


P = 


P - cpr - 1 

(X — Xg ) /0 g 


R0< 


define the largest possible P that can be achieved by a short bubble at a 
given ^g* That is, if c^^ from equation (31) is substituted into equa- 
tion (32) and the resulting P tabulated as a function of (x - Xg)/dg for a 
particular *^gf then an extremum of P is found. For the Reynolds numbers 

R0g = 100, 200, 300, and 400, we find that the extreme values of P are 

-0.0996, -0.1868, -0.2732, and -0.3593, occurring at (x - Xg)/0g = 103, 111, 
115, and 117, respectively. A faired curve through these points is shown in 
figure 8 and exhibits suprisingly good qualitative agreement with Caster's 
enpirical curve. 

The agreement between theory and experiment is particularly pleasing in 
view of the fact that all of the empirical input to the present calculation 
method was obtained from measurements of transition in free shear layers in 
zero pressure gradient at a remote distance from any wall. Our calculations 
therefore lend support to the view that the processes of bubble bursting and 
shear layer reattachment are linked to the processes of laminar-turbulent 
transition and shear layer growth in a fundamental way. 
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COMPARISON WITH EXPERIMENTAL BUBBLE LENGTHS 


As a caution against overly optimistic conclusions drawn frcan the present 
analysis, we present in figure 9 the results of an attempt to verify the 
reattachment condition (31) directly. In the figure, we have plotted the mea- 
sured pressure recovery coefficient a and length (x^ - Xg)/0g for a number 
of separation bubbles fran various sources. On the same scale, we have drawn a 
copy of figure 7, which represents the nondimens ional stagnation pressure along 
the separation streamline for a plane layer. If equation (31) held exactly, 
then all of the experimental points would lie on the dashed curves. 

In fact, the experimental points are scattered quite widely about the 
analytic curves. More serious than the scatter in CT, however, is the discrep- 
ancy between the measured length of the bubbles, which are typically in the 
range 150 = (x^ - Xg)/6s = 300, and the values calculated from our analytic 
derivation of Caster's bursting line. We have indicated the latter points in 
figure 9 by "plus" symbols, which lie, to within graphical accuracy, on the 
family of straight lines drawn through the origin and the points of tangency 

2 

with the various curves of (u/U)^=t|;^ for the plane layer. As can be seen, 

the calculated bubble lengths are typically a factor two or three too small. 

This error is significant, and indicates that the present theory is not suf- 
ficiently refined for practical use in the quantitative prediction of bubble 
lengths . 

It is not clear which of the many effects neglected in the theory are 
responsible for this discrepancy. We may speculate, however, that the neglect 
of changes in the shear stress gradient 9T/9n on the separation streamline 
during reattachment is important. Some such change must occur, since in the 
real flow the shear stress on 4^ = v|;g vanishes at reattachment, which places 
the reattachment point below the point of maximum T in the layer. The crude 
theory developed above assumes that as the separation streamline gets swallowed 
up into the underside of the shear layer, the shear stresses on the streamline 
tend asymptotically toward the value of the center of the layer (which is the 
maximum value). It follows that 9T/9n is underestimated by the crude theory. 
This implies (from equation (25)) an unrealistically low prediction of the stag- 
nation pressure H. It is not obvious how the present theory could be refined 
to take such effects into account. 


CONCLUSIONS 

From our investigation of the j^ysics of separation bubbles, we may state 
the following tentative conclusions: 

1 . The form of the static pressure distribution about a short separation 
bubble can be explained qualitatively in terms of boundary layer concepts. In 
particular, the sharp corner in the pressure distribution at the rear of a short 
bubble is more properly interpreted as the point where the lower edge of the 
reattaching shear layer grazes the surface than as the point of transition from 
laminar to turbulent flow, as has traditionally been assumed. 
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2. The observation that the static pressure distribution about the body 
on which a short separation bubble occurs very nearly coincides with the 
inviscid pressure distribution at all points other than those between separa- 
tion and reattachment leads to a plausible hypothesis regarding the factors 
controlling bubble bursting and a simple but crude calculation method for esti- 
mating the length of a short bubble when it is possible for one to exist. Spe- 
cifically, the reattachment point Xj- coincides with the earliest downstream 
station x at which the stagnation pressure on the separation streamline is 
large enough to support an abrupt rise to the local inviscid static pressure 
at X. 


3. The attempt to approximate the rise in stagnation pressure on the sepa- 
ration streamline in a bubble by the corresponding rise along a streamline 
entrained into the underside of a laminar mixing layer undergoing transition 
leads to a reasonably accurate prediction of Caster's bursting line, but signif- 
icantly underestimates the length of short bubbles in general. This inaccuracy 
is probably due to the failure of the calculation method to take account of 
changes in the viscous and turbulent shear stress profiles occurring prior to 
reattachment that produce (according to equation (25)) an increase in the actual 
stagnation pressure along the separation streamline beyond what exists in a free 
shear layer. 

4. The present results are in every way consistent with and support the 
view of Cebeci and Bradshaw and others that reattachment is a direct consequence 
of the increased entrainment of fluid into the underside of the growing sepa- 
rated shear layer downstream of transition, the entrainment causing the shear 
layer to suck itself back down onto the surface. 
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APPENDIX 


We rewrite equation (15) in the equivalent form 


/^OO 

I f'- 


f ) - (1 - f’) 


2 ] dn = r (f )2 dn 

'-'-^00 


Substituting equations (12) and (13) into equation (A1 ) gives 

|[i - f (0)] ii - [1 - f (0)J^^1 = 13 
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and 


I 3 = rV- sech2 
J_ooV 


bn + tanh"^ f — 

[f3j 


b dn = 


4J3 - 9 


(assume b < 0 throughout) 

The right-hand side of equation (A2) can be rewritten by means of equation 
as 


I 3 f -f ' (0)3 2fi 1 

= I3 - 

b 1 - f'(0) a 


Substituting the identity 


-f (0)3 = [i - f ( 0 )] ^ - 3[l - f ( 0 )] ^ + 3[l - f ( 0 )] - 1 


(Al) 


(A2) 


(A3) 


(A4) 


(A5) 


(14) 


(A 6 ) 
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into equation (A6) and combining with equation (A2) gives the cubic equation 
[l - f(0)] 3^-12 - I3 - f'(0)]2|li + 313 


P / 2i/tt\ 2Jtt 

+ [1 - f • (0)] 1^-313 — = 0 


(A7) 


which may be solved numerically by Newton's method to give the solution 
1 - f ' (0) = 0.41077, or f ' (0) = 0.58923. 

The constants a and b may be calculated from the normalization condi- 
tion (16) in the following way. The first integral on the right-hand side of 
equation (16) may be written 



(f')2j dn = 


f (0) 


f (0)2 
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I3 


-f (0) 

1 - f'(0) 



^|f'(0) I4 - [f'(0)J^I^ 


(A8) 


where we have used equation (14), and where 


I4 



bn + tanh“^ 



b dn = 


s[3 - 3 
2 


(A9) 


The second integral on the right-hand side of equation (16) is just the left- 
hand side of equation (A2) . According to equation (16), the sum of these two 
integrals equals -f(-“) = 1.2386, which immediately gives a = 0.88544. Equa- 
tion (14) then gives b = -0.30160. 
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Figure 1.- Characteristic pressure distribution beneath a 
separatic«i bubble of the "short" type. 


EDGE OF SHEAR LAYER 


MEAN STREAMLINES 



Figure 2.- Two-dimensional steady separation frcm a 
solid wall showing entrainment of fluid into the 
underside of the shear layer due to spreading in 
the downstream direction. 
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- EDGE OF SHEAR LAYER 


MEAN STREAMLINES 



Figure 3.- Recirculating flew pattern within a bubble 
showing splitting of the shear layer at reattach- 
ment to supply the entrainment (vertical scale 
greatly exaggerated) . 



Figure 4.- Pressure distribution inferred from figure 3. 
Small entrainment velocities beneath shear layer imply 
small pressure gradients between S and K. Steep 
gradient between K and R due to splitting of shear 
layer at reattachment and gradient of stagnation pres- 
sure across lower half of shear layer. 
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- - EDGE OF SHEAR LAYER 



Figure 5.- Nomenclature for free shear layers in 
zero pressure gradient. 



Figure 6.- Verification of empirical formula for growth of momentum 
thickness downstream of transition. Data frcan reference 12. 
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Figure 8.- Caster's bursting line derived from present theory 
which is based entirely on data frcxn free shear layers. 



Figure 9.- Attempt to verify reattachment condition (31). 
Measured lengths and pressure recovery coefficients 
frcan various sources: □, Gaster (ref. 14) "series I"; 

A, Gaster "series II"; O, Van Ingen (ref. 11); V, Tani 
(ref. 2, fig. 17, based on data from Gault, ref. 10). 

If equation (31) held exactly, then all measured points 
should lie on dashed curves. Note: each set of data 

represents a series of runs at various tunnel speeds, 
which, to avoid cluttering the figure, have not been 
shown. 
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WING SHAPE OPTIMIZATION FOR MAXIMUM CROSS-COUNTRY SPEED, 


WITH MATHEMATICAL PROGRAMMING 
Gunter Helwig 

Rensselaer Polytechnic Institute 
SUMMARY 

A computer program was developed to calculate numerically the speed and 
circling polars of an aircraft when the lift and drag characteristics of the 
wing airfoils are known. The planform of the wing is described by variables 
which are optimized so that the cross-country speed of the glider is maxi- 
mum for the particular type of thermal model. Two thermal models will be 
compared and it can be shown that with a greater wing area than now normally 
used the performance can be increased. 
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SYMBOLS 

half span 

section drag coefficient 

section lift coefficient 

total drag coefficient of the aircraft 

total lift coefficient of the aircraft 

induced drag coefficient 

parasite drag coefficient 

acceleration of gravity 

height 

radius of thermal 

area 

time 

net climb rate of the glider 
sink rate for the speed polar 
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V 

sc 

V 

thermal 

V 



W 

Z 


X . 
1 

a. 

X 

P 


sink rate for the circling polar 
vertical velocity of thermal 
airspeed 

average cross-country speed 

cross-country speed for the i-th thermal 

total weight 

object function 

variables 

weighting factors 

air density 


(p sweep angle 

A aspect ratio 

bank angle 


INTRODUCTION 

The most important part of an airplane is a well designed wing. Only a 
few variables are necessary to describe the planform (fig. 1). The halfspan 
b and the sweep-angle are fixed. The airfoils with their lift and drag 
characteristics must be prescribed. The total lift of the airplane is 


r = r + C jj- 

wing tail 


( 1 ) 


The C is calculated either by Multhopp's method (ref. 1) or, for swept wings. 


by Truckenbrodt 's method (ref. 2). The tail lift is in general so small that 
it can be neglected. From the lift calculation the induced drag (ref. 1) 
is obtained, too. The total drag is i 




+ 


profile 


s. ^ 
1 


Do 


( 2 ) 
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(3) 


Do 




Sp (p = 1,2. .P) 
P 


is the drag coefficient of the p-th part, with surface area S^. The 

profile drag coefficient C can be calculated from the measured 

profile 

Reynolds-number-dependent c^-versus-c^ plots (fig. 2) . 


profile n 


N) 


(4) 


From fig. 2 (measured data), fig. 3 can be determined, and by linear interpo- 
lation, the local c, at the strip n with surface S is obtained. The air- 
’ d n 

n ^ 

craft equivalent parasite area y Sp is assumed to be constant for a 

P ^ 

given aircraft. 


SPEED AND CIRCLING POLARS 

With the weight as the fifth variable (x^ = weight), the principal per- 
formance characteristic of the glider, which is its cross-country speed, can 
be calculated. The speed polar is given by the equations 
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CROSS-COUNTRY SPEED 


The net rate of climb of the glider in the thermal (fig. 4) is 


V (r) = v^, 1 (r) - V (r) 

c thermal sc 


( 9 ) 


The maximum climb rate must be calculated as a function of the diameter of the 
thermal. With these two polars (speed and circling) the cross-country 
speed can easily be determined (fig. 5) , 


V 


R 



( 10 ) 


(D = distance; t = time from point A to C) . The loss in height is 

H = t,v 
1 s 


The gain in height is 


H = 


• V 

c 


Therefore, 



( 11 ) 


( 12 ) 

(13) 


and with the time t^^ to go from point B to C, and the total time t, re- 
spectively, 


^ D 

^1 V ’ 


(14) 


t = ti -1 t2 


(15) 


the cross-country speed by using the equations (13) to (15) and (10) is 
finally 

V • V 

V = 

R V -t- V 
c s 


(16) 


THE OBJECT FUNCTION 

The program can be used with different types of thermal models. Two 
examples are given in fig. 6 and fig. 7 (ref. 3,4). As object function for 
the optimization, an average cross-country speed is defined for each thermal 
model 
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(i = 1,2,3 for Carmichael) 
(i = 1,2,.. 4 for Horstmann) 


(17) 


Z 


1 1 


Z a . = 1 and 1 > a, >0 (18) 

1 — 1 — 

1 

The 's are the cross-country speeds from the i-th thermal and is its 
i 

weighting function. The a^’s must be chosen by the designer and then the 
optimization is achieved for the particular distribution of thermals assumed. 


THE OPTIMIZATION PROCEDURE 

The maximum cross-country speed is calculated by a penalty optimization 
program (ref. 5,6), so that constraints like ~ 1»2»«*5) (a^ = 

lower bound) or others can be observed. Using the performance polars 
from equations (1) and (2), two sub-optimization problems are 
solved: (a) determination of max v^ by using C^^ and r as variables in 

equations (7) to (9), and (b) determination of max V (eq. 16) by using 

C as a variable through equations (5) and (6), with max v fixed. V 

i-» C K, . 

has to be calculated for each thermal (index i) . Then Z (eq. 17) can 
be calculated and optimized. 


RESULTS 

Only gliders with a span of 15m = 2b and (J> = 0° were optimized. In 
diagrams 1 and 2 an optimization is shown with the profile FX-61-184 (ref. 7) 
for two thermal models with different sets of a^'s. For a convenient repre- 
sentation (this is only one possibility), = 02 was used for the 

Carmichael thermals and a, = a., and a„ = a, for the Horstmann thermals. For 

12 3 4 

comparison, the performance of the D-38, an almost optimized competition 
glider in the 15m class from the Akaflieg Darmstadt (ref. 8), is shown in 
both diagrams. The curves were calculated by changing only the weight of 
the glider. The difference between the two curves gives the gain in per- 
formance over the D-38. On the average, the gain in performance is more 
than 2% for both thermal models. The Carmichael thermals deliver a higher 
cross-country speed because it is narrow and strong. This leads to a higher 
optimal wing area for Carmichael thermals compared to the Horstmann thermals 
(see table 1 and 2). In general, 2 the optimized wing area is significantly 
larger than the normally used 11m area for rigid 15m wings. 
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In diagram 3 the wing loading and the aspect ratio are plotted versus 
the variable sets of for both thermal models, using the numbers from 

tables 1 and 2. The difference in the aspect ratio, AX , is about 1 and 
therefore the wing area difference, AS, is about Im^. The average wing 
area is 14. 4m^ for the Carmichael and 13. 5m^ for the Horstmann thermals. 
Clearly the influence of the shape of the wing planform is small. The 
aspect ratio is almost constant from weak to strong thermals. This is true 
for both thermal models. The wing loading is the dominant factor influencing 
the cross-country speed. The dotted line (diagram 2) which was obtained by 
changing the weight, using optimal design for all a^'s equal, shows the im- 
portance of the weight (wing loading) but a relatively small influence of 
the wing shape. The proper wing loading, again, is more important for the 
Horstmann model (greater gradient) than for the Carmichael model for 
maximum cross-country speed. 

Because of the dominating rule of the weight for the gliders , the 
curves in diagram 4 were calculated for various masses by using the 
optimal design with all a^'s equal. The optimal variables are given 

in table 3. The comparison between the 7 profiles shows that the best rigid 
profile is the FX-61-184, and the best profile with flaps is the FX-K-170 


CONCLUSION 

An increase of performance of more than 2% is possible with an increase 
in wing area. It is not necessary to design a glider for extreme (strong 
and weak) thermal conditions. The use of average weather (all a^'s equal) 

for the optimization of a glider results in an almost optimal design for 
all weather conditions. But the glider must be built as light as possible 
and should be able to carry up to 150 kg water ballast. The thermal model 
(Carmichael, Horstmann, or perhaps others) plays a minor part in the design 
but is extremely important for the proper choice of the water ballast to 
maintain maximum performance. Here it matters whether to believe in the 
Carmichael, Horstmann, or other thermal models. The program can then be 
used for a single variable optimization (the weight) to calculate the 
optimal water ballast curves for any type of weather condition. 
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Table 1: 


Optimal Variables for 15m Gliders with the Profile FX-61-184 
(Carmichael thermals, = 0.04m^) 


OPTIMAL 



OPTIMAL 15m 

GLIDERS 


D-38 

VARIABLES 

a^a.^0 

03=1 

a^=a 2 = 0.2 
03 = 0. 6 

01=02=03 

Oj^=O2=0. 45 
03 = 0.1 

o^=O 2 = 0 . 5 
03=0 


(m) 

1.31 

1.19 

1.23 

1.14 

1.19 

0.94 

y ^2 (“) 

1.11 

1.07 

1.09 

1.04 

1.02 

0.753 

x^ (m) 

0.49 

0.47 

0.51 

0.45 

0.45 

0.376 

x^ (m) 

3.78 

4.17 

3.87 

3.62 

4.25 

4.5 

(kg) 

249 

289 

327 

360 

381 

300 

DERIVED 

VALUES 







S (m^) 

14.89 

14.6 

14.76 

13.65 

14.2 

11 

\ 

15.11 

15.4 

15.22 

16.48 

15.85 

20.45 

max 

34.3 

34.8 

34.8 

35.7 

35.5 

37 
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Table 2: Optimal Variables for 15m Gliders^with the Profile FX-61-184 

(Horstmann thermals, = 0.04m ) 


OPTIMAL 

VARIABLES 


OPTIMAL 

15m GLIDERS 


D-38 

ct^=a2=0 

a_=a,=0. 5 
3 4 

oij^=a2=0 . 25 

a~=a,=0. 25 
3 4 

a^=a2=0- 4 

a„=a =0.1 
3 4 

a^=a2=0*5 

a„=a =0 
3 4 


(m) 

1.32 

1.19 

1.149 

0.999 

0.94 

X 2 (m) 

0.99 

0.884 

1.012 

0.995 

0.753 

(m) 

0.507 

0.384 

0.394 

0.389 

0.376 

(m) 

4.05 

4.48 

4.28 

3.679 

4.5 

(kg) 

254 

328 

390 

432 


DERIVED 






VALUES 






S (m^) 

14.53 

13.2 

13.5 

12.64 

11 

A 

15.48 

17.04 

16.67 

17.8 

20.45 

max C^/Cjj 

34.6 

35.9 

36.0 

36.7 

37 
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Table 3: Optimal 15m Gliders with Different Profiles 

(Horstmann thermals, = 0.04m^) 





“l “2 

“3 “4 







Profile 





v£) 

<N 

CO 
1— 1 

r-H 

> 

vO 

rH 

1 

CO 

00 

rH 

rH 

CO 

0 

m 

0 

OPTIMAL 

VARIABLES 

1 

o 

VD 

1 

X 

fn 

I 

rH 

vO 

1 

g 

1 

VO 

VO 

1 

X 

1 

rH 

VO 

1 

X 

pH 

1 

1 

Xi 

1 

X 

1 

X! 

1 

x; 

1 

X 

(m) 

0.94 

1.07 

0.98 

1.19 

1.18 

1.03 

1.077 

(m) 

0.76 

0.927 

0.919 

0.884 

0.99 

0.866 

0.911 

(m) 

0.379 

0.395 

0.37 

0.384 

0.38 

0.39 

0.37 

x^ (m) 

4.37 

4.08 

3.37 

4.48 

4.9 

4.26 

4.14 

x^ (kg) 

290 

311 

316 

328 

348 

344 

367 

DERIVED 

VALUES 








S (m^) 

10.76 

12.8 

11.52 

13.2 

14.08 

12.1 

12.74 

A 

20.9 

17.58 

19.53 

17.04 

15.98 

18.6 

17.66 

max 

36.4 

35.7 

34.3 

35.9 

38.5 

37.6 

37.1 
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X, = chord at the root 
Xg = chord at the break 
Xj = chord at the tip 
X 4 = span to the break 




Fig. 2‘. Reynolds-number-dependent 
lift and drag characteristics 


Fig. 3'. Reynolds-number-dependent 
Cj for C|=constant 

^profile 
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V 




Fig. 5; Typical Flight Section 
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V 


thermal 


[m/s] 



Fig. 6= Carmichael Thermals 
T| = strong 

T^= wide 

T^= weak 


V,. ,[m/s] 

thermal 


6 - 



Fig. 7: Horstmann Thermals 
T| = narrow, strong 

T^ = wide, strong 

T^= narrow, weak 

T^= wide, weak 



Diagram 1: Optimal performance of a 15m glider with the profile FX-61-184 

for the Carmichael thermals. The numbers on the curves are total 
weights of the gliders in kg. = O.OAm^ = const.) 



(28.1) (31.3) (34.4) (37.5) (mph) 

Diagram 2: Optimal performance of a 15m glider with the profile FX-61-18A 

for the Horstmann thermals. The numbers on the curves are total 
weights for the gliders in kg.' (C^ S = O.OAm^ = const.) 
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Wing Loading, 
Kg/m^ 


Carmichaei 
Thermais Oj 



Aspect Ratio X 


A Horstmann — 
o Carmichael — 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 a.+a^ 

1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 Thermair 


Diagram 3: Comparison Between Carmichael and 

Horstmann Thermals. 



Diagram 4: Performance of 15m gliders using different profiles for the Horstmann 

thermals. Curves are calculated by using the optimal design ( ^ design 
point) with all o. 's equal. The change of the total mass from point 
to point is 50 Kg (Cjj^S = 0.04m^) 
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FURTHER DEVELOPMENTS IN SIMPLE TOTAL ENERGY SENSORS 


Oran W. Nicks 
Langley Research Center 


SUMMARY 


In 1976, research results were published on a simple total energy probe 
concept using principles of laminar flow around a cylinder. A number of 
probes employing these principles have been built. Additional tests have 
been conducted to further support earlier findings and options for probes 
made of a single bent-up tube. 

Total energy pressure relationships are reviewed and flow fields around 
cylinders normal to and inclined to the flow are described. A variety of 
bent-up probe configurations were tested to explore variations in geometry. 
Test results are presented on the effects of sensor length, hole location, 
and angle of sweep. Comparisons are made with other probe tests reported 
in the literature. 

A brief summary of damping restrictors and their use in filtering gusts 
is presented. Flow field effects, indicating the variables involved for 
different mounting locations on aircraft, are discussed. 


INTRODUCTION 


Since the 1976 publication of research results on a simple total 
energy sensor using principles of laminar flow around a cylinder, reference 1, 
a number of developments in their application have occurred. One objective 
of the research was to help sailplane owners improve their soaring instru- 
mentation with a simple "do it yourself" design for a total energy sensor, 
references 2 and 3; the principles and broader applications are outlined in 
the Patent description, reference 4. 

Many sensor probes have been made embodying the principles advanced; 
some are in use and reported to be performing quite satisfactorily, a number 
of modifications have been reported to suit the individual application, and 
some difficulties have been encountered with home-built sensors and 
applications for various reasons. In order to expand the general reference 
knowledge of the principles and sensitivities involved, this report provides 
information and data on further analyses and tests of simple total energy 
sensors using principles of laminar flow around a cylinder. 
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TOTAL ENERGY - PRESSURE RELATIONSHIPS 


First, it is appropriate to briefly review the physical relationships 
of gliding flight which make a total energy sensor useful. Figure 1 shows 
how the useful total energy of a sailplane can be visualized in terms of 
altitude and velocity. The potential energy is directly proportional to 
altitude, the kinetic energy is directly proportional to the square of the 
velocity, and the sum of potential and kinetic energies is the total energy. 
At a particular instant, the best indications of the total energy state are 
provided to the pilot by the altimeter and airspeed indicator. For real-time 
energy management, however, the rate of change in this total energy state 
is most important. The basic requirement for a total energy sensor is to 
help provide the pilot such an indication. 

When a sailplane dives, it gains kinetic energy at the expense of 
potential energy; the opposite occurs in a zoom; however, the only change in 
total energy during such maneuvers is caused by the drag of the sailplane 
and atmospheric energy variations, if any. The steady state drag effects 
are proportional to the drag times the velocity squared as indicated by a 
sailplane "polar" (where "polar" is defined as sink rate vs. airspeed), and 
the most significant atmospheric effects are rising and sinking air currents. 
Secondary effects on changes in total energy are caused by drag increases 
during rapid accelerations such as sharp pullups and tight turns, and 
horizontal wind gradients or shears can be significant near the ground. 
However, these are generally ignored as secondary during non-aerobatic 
flight at soaring altitudes. Thus, in simple terms, the rate of change in 
useful total energy may be indicated by a simple variometer instrument, if 
the total energy sensor connected to it provides proper compensation for 
exchanges in velocity and altitude. 

We have said that the total energy TE of a sailplane of mass M 
gliding at a velocity V and at altitude H is: 

TE = MgH + 1/2 MV^ 

The energy per unit mass of the sailplane can be written: 

TE = gH + 1/2 (1) 

Differentiating to obtain the rate of change gives: 

d(TE) = gdH + 1/2 d(V^) (2) 


Assuming constant altitude, making use of the relation dP^ = 


pgdH, 


where P is the ambient static pressure, and dq = 1/2 pd(V ), where 
° 2 

q = 1/2 pV we have: , 

_ -dP 1/2 d(V‘^) 

d(TE) = 5. + 

P P (3) 
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Assuming the total energy of the sailplane remains constant, 
d(TE) = 0 and eouation (3) becomes: 


dP^ - dq = 0 

Integration of equation (4) yields the following: 


(4) 


- q = p. 

0 ^ s 


(5) 


where is the constant of integration. 

P^ obviously is a pressure and must be positive, since for 

q = 0, Pg “ Pq • Changes in the pressure P^ therefore provide an 

indication of changes in total energy of the sailplane. Putting equation (5) 
in coefficient form: 


As indicated in reference 1, a number of sources pointed to the nearly 
correct pressure relationship on the downstream side of a right circular 
two-dimensional cylinder aligned normal to the flow, when the size of the 
cylinder and flow characteristics produce Reynolds numbers based on diameter 
from about 5,000 to perhaps 350,000. Within this Reynolds number range, the 
flows are described as subcritical, well established laminar flows before 
separation occurs. 

Figure 2 shows the nature of the streamlines calculated for the flow 
around a cylinder in this regime, reference 5. The separated flow region 
results in relatively constant base pressures over the aft 110° to 160° of 
the cylinder. This large, relatively constant pressure region is the reason 
that a sensor made from a cylinder is insensitive to angles of yaw or 
circumferential hole position accuracy within this region. The separated 
flow does tend to fluctuate, however, and a high frequency vorticity can be 
sensed with dynamic instrumentation. Fortunately, these rapid pressure 
fluctuations can be damped and need not compromise a total energy sensor 
output for practical application. 

A sample plot of a typical pressure distribution, figure 3, is shown 
to illustrate the nature of the pressure distribution in the Reynolds number 
range of greatest interest. This plot shows the relative sameness of the 
pressures on the aft side of the cylinder, corresponding with the separation 
region downstream in the streamline diagram. Such data provided the 
inspirations to use a small cylinder as a means of achieving the desired 
pressure relationships for a total energy sensor. 



( 6 ) 


THE FLOW OF AIR AROUND A RIGHT CIRCULAR CYLINDER 
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THE INCLINED THREE-DIMENSIONAL CYLINDER AS A TOTAL ENERGY SENSOR 


Reference 1 described experimental studies which included the discovery 
that pressures on the aft side of a cylinder could be modified near the flat 
end of a three-dimensional cylinder to the value desired for total energy 
relationships. Further testing showed that sweeping the cylinder forward 
into the airstream about 20°, combined with a specific hole location relative 
to the end of the cylinder, produced the proper pressure relationship with 
an insensitivity to flow inclinations of at least - 10°. This is more than 
adequate for use during soaring flight, as pitch variations for sailplanes 
of less than 6° to 8° are normal. 

The combination of insensitivities to sideslip and angle of attack is 
especially desirable features of the simple probe. In addition, data were 
provided which indicated that the sensor holes could be of various diameters, 
saw slots, or multiple orifices as long as the average dimensions with 
respect to the cylinder end were maintained. Two specific total energy 
probe configurations were described which had been based on laboratory 
results and flight tests. 

The most scientifically significant findings of reference 1 are 
summarized below for reference during the discussions to follow: 

1. Flow normal to a right circular cylinder at subcritical Reynolds numbers 

produces pressure coefficients very close to the value needed for 
useful total energy rate of change indications. 

2. For the velocity and altitude operating range of sailplanes, practical 

cylinder diameters of about 4.76 mm (3/16-inch) to 6.35 mm (1/4-inch) 
provide Reynolds numbers within a range from about 8,000 to 30,000, 
where sensor pressure coefficient and drag coefficient remain 
practically constant. 

3. Orifices on the downstream side of a cylinder provide pressures that 

are relatively insensitive to sideslip angles. 

4. Whereas pressures on the downstream side of a two-dimensional cylinder 

produce coefficients that tend to be too negative, it is possible to 
obtain predictably biased pressures with three-dimensional flow effects 
on a practical sensor by locating rearward facing orifices a given 
distance from the end of the cylinder; furthermore, the variation in 
pressure coefficient with hole distance from the end of the cylinder 
tends to be linear in the region of interest for coefficients near 

Cp = -1.0. 

5. Pressure coefficients on the aft side of such a cylinder remain 

relatively constant over a range of forward sweep angles from about 
10° to 30°; thus, a nominal cylinder orientation of 20° forward sweep 
provides a sensor with i 10° insensitivity to pitch changes. 
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BENT-UP PROBES - PURPOSE OF TESTS 


The fin mounted probe configuration described in references 1, 2, 3, 
and 4 was made by joining two straight sections of tubing. Many home- 
builders have made total energy probes from a single piece of tubing bent to 
provide the 20° forward swept portion, thus eliminating the two-piece 
manufacturing difficulties and the chance for leaks at the intersection of 
the two tubes. It is understood that many of the probes made in this way 
have been made of 6.35 mm (1/4-inch) steel tubing commonly used for hydraulic 
brakelines or fuel lines. This size tubing offers the strength and stiffness 
to support an extension of about 40 cm (15 inches) ahead of the fin leading 
edge. 


Homebuilders who had difficulties with such probes report that pressure 
coefficients too low in absolute value were achieved, thus causing under- 
compensation in flight. One of the suspected causes of the low pressure 
coefficients was the likelihood that the bent-up sections of tubing were 
not long enough to sustain the two-dimensional flow field below the orifice. 
Further wind tunnel testing has been conducted on simple bent-up probes to 
determine the effects of various probe lengths on pressure coefficients and 
to better define the geometry of suitable probes made in this manner. 


WIND TUNNEL TEST SETUP AND PROBE CONFIGURATIONS 


The small wind tunnel used for the tests outlined in reference 1 was 
used for these tests. An atmospheric tunnel with velocities on the order of 
20 meters (60 feet) per second, it provided Reynolds numbers based on 
diameter with the 6.35 mm (1/4-inch) tubing of 8,000 to 10,000. 

The probe test setup was made to allow a common set of sensors to be 
used. This insured that no differences in results were caused by 
manufacturing differences on the sensors or hole location geometry. 

A simple mounting arrangement allowed changes in the probe angle of 
attack so that sensitivity to sweep angle or angle of incidence could be 
determined. The angle of attack could be varied over a range of 35°. Since 
the probe had been designed with a forward sweep angle of 20° as the nominal 
mounting position, this meant that the probes were tested with forward 
sweep into the airstream over a range of +5° to +40°; for simplicity, all 
data are presented on that basis. The mounting and angle of attack changing 
system allowed the sensor to remain in the core of the wind tunnel flow where 
flow was uniform and velocities constant. As each probe configuration change 
was made, leak tests were performed to insure sealed joints. 

Three types of configurations were tested with a combination of four 
sensors and two extenders, giving a total series of 13 configurations. These 
allowed a range of geometric parameters to be tested including sensor length, 
hole position from the probe end, hole position circumferentially, and two 
holes at a fixed orientation. These are shown in figure 4. 
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WIND TUNNEL RESULTS 


As indicated in figure 5, variations in the length of the straight 
section of the sensor above the bend did have an effect on the pressure 
coefficient for sensors with the end geometry described in reference 1. 

Sensor lengths of 7 and 9 diameters produced pressure coefficients lower 
than the desired Cp = -1.0, whereas Y/D's of 11 and 13 both produced 

Cp = -1.0. These data support the findings from earlier tests where 

Y/D's of 12 or greater gave satisfactory results. Since the drop-off in 
coefficient apparently began to occur for lengths between 9 and 11 diameters, 
sensor lengths of 11 diameters or greater from the bend appear necessary to 
insure the proper flow effects at an orifice located two diameters from the 
probe end. 

These data also confirm that a nominal forward sweep angle of 20° is 
a good choice to allow variations in flow direction which may result from 
downwash, slight mounting misalignments, or attitude changes during flight, 
without effects on the sensor pressure. 

After it was determined that lengths shorter than 11 diameters did not 
produce the desired pressure coefficient with the hole at X/D = 2.0, an 
experiment was performed to determine whether locating the sensor hole 
nearer the end of the probe might counteract this effect by capitalizing on 
the varying effects of hole position discovered earlier. Using bent-up 
probe configurations having a fixed length/diameter ratio of Y/D = 7, 
results were obtained with three hole positions as indicated in figure 6. 

These data show that a C = -1.0 should be achievable for a hole position 

P 

between 1.5 and 1.75 diameters from the sensor end; however, it is also seen 
that the range of insensitivity to sweep has lessened somewhat when compared 
to probes with a hole location 2.0 diameters from the end. 

In a personal communication, Frank Irving of Imperial College referred 
to data on pressure distributions around cylinders normal to the flow showing 
less variation in pressure coefficient at circumferential hole locations 
other than 180°.^ He reported that he had tested probes with two holes 

located at 0 = - 130° which gave good results. 

His comments led to a review of earlier data from references 1, 6 and 7 
for cylinders normal to the flow. A slight trend toward greater dispersion 
of coefficient at 0 = 180° as a function of Reynolds number was evident, 
but perhaps more interesting are trends from reference 8 shown in figure 7 
for two-dimensional cylinders swept at various angles to the flow. The trend 
toward more negative pressure coefficients at higher 0 values is seen, 
along with the interesting fact that differences in coefficient for various 

sweep angles are less at lower 0 values. This led to the conclusion that 

a broader range of insensitivity to incidence changes might be achieved if 
circumferential hole positions less than 0 = 180° were used. 
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Figure 7 shows the variations in pressure coefficient with sweep to be 
very slight at 6 = 140°; for symmetry in a probe, it would be desirable to 
have two holes at a ±0 ^alue. Since it is much easier to position two 
holes 90° apart at 6 = - 135°, this was chosen as a practical compromise to 
keep the two-hole probe as simple as possible to construct. 

Tests with the 135° azimuth hole position were conducted using the same 
probe sensor sections tested earlier with the holes reoriented to 9 = 135°, 
thus insuring consistent and comparable results with the same sensor sections 
that had been tested at 9 = 180°. The tests were conducted at various 
X/D's for only two Y/D values of 7 and 11. These two probe sensor lengths 
were used because the results from earlier tests, shown in figures 5 and 6, 
indicated that coefficients of the desired values could be achieved for these 
lengths. 

The data from tests of the shortest bent-up probe, Y/D = 7 
configuration, are shown in figure 8. As indicated, the holes at the 
9 = 135° position do increase the range of insensitivity to forward sweep 
over the range of interest. These data show the best X/D to be 1.75, 
producing a coefficient Cp = -1.0 over a 20° range of sweep angles. 

However, the mid-range appears biased such that 25° of forward sweep might 
be better than 20° as a nominal. 

In figure 9, similar data are shown for the sensor length, Y/D = 11, 
found earlier to give consistent results with previous tests of probes 
employing the same and greater lengths. By comparing these data with 
figure 5, it is seen that the configuration with the 0 = 135° hole position 
gave results very comparable to the 9 = 180° hole position, except for the 
slightly lesser sensitivity to forward sweep. 

Finally, as a confirmation check of an actual two-hole configuration, 
a new two-hole probe was made for Y/D = 7.0, X/D =1.5, 9 = -135°. The 

results of this test showed consistent results with the single hole tests 
for 9 = 135° as shown in figure 10. This was expected since no variation 
in yaw occurred during tests of the single hole probe, but it provided 
positive assurance. The insensitivity to sweep was significantly extended 
for the short probe having two holes, although the particular hole location, 
X/D = 1.50, resulted in slightly over-compensating coefficients. However, 
as indicated in figure 8, the desired coefficient should be obtainable by 
locating two holes at X/D = 1.75. Thus, it has been shown that for a short 
sensor section with a straight length of 7 diameters, a two-hole configuration 
will produce good compensation. In general, it has been found that the closer 
the hole is to the end of the probe, the more sensitive the end effects and 
more rapid changes in coefficients are likely to occur for small variations 
in the bevel or chamfer. This tends to make tolerances more important; for 
this reason, it appears that more consistent results can be expected if 
bent-up probes are made with sensor lengths of 11 diameters or greater with 
hole positions 2 diameters from the end. 
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Perhaps a further word should be said about the shape of the end of the 
sensor and its effects. References 1, 2, 3, and 4 pointed out that the three- 
dimensional effects of the probe end were affected by the amount of rounding 
or chamfer of the end. To obtain consistent results, a squared-off end 
with a very slight chamfer to "break" the sharp edge was recommended. 

In reference 9, Wells discussed the matter of beveling edges and described 
a method of making a chamfering tool for this purpose. He indicated values 
of about 0.066-0.018 mm (0.004-0.007 inch) are typical. Such precision in 
chamfering has not been found essential, but experience has shown that 
rounding off the edges too much tends to produce over-compensation, and care 
is recommended in beveling the edge. While it is difficult to specify and 
measure dimensions for the beveled edge, it is believed that chamfered 
surfaces of about 0.013 mm (0.005-inch) are suitable. 


COMPARISON OF VARIOUS CONFIGURATIONS 


In addition to a number of informal reports of experiences with total 
energy probes based on the principles outlined in reference 1, several 
results have been presented in periodicals which are worthy of mention. 

In reference 9, Wells discussed a method used for calibrating total energy 
probes and some results of tests. Performance data were not shown; however, 
the article indicated satisfactory results with a bent tube approach, 
although the benefits of forward sweep were questioned. The importance of 
care in beveling the edge was specifically discussed. 

Diplom-Physiker Westerboer (ref. 10) described probes made and tested 
in West Germany of both fin and fuselage mounted types. He indicated that 
good compensation was achieved and also mentioned other results obtained 
in Europe. He specifically mentioned that an experimenter from the 
Braunschweig Akaflieg had also found 60° to 80° forward sweep (equivalent 
of 10° to 30° measured from vertical) to be optimum. Dimensions given 
for a bent-up tube version included a straight section for the sensor 
portion of about 10 diameters. The test results reported in figure 5 show 
that greater than 9 diameters are required to achieve the desired sensor 
pressure coefficient; thus, it would seem that Diplom-Physiker Westerboer 
also confirmed the suitability of this geometry for a bent-up probe con- 
figuration. 

An article (ref. 11) by Charles W. Shaw described a probe made for a 
nose mounting instal lation. The report stated that excellent compensation 
was achieved, along with improved response rate, and no effects due to 
sideslip. The probe was mounted on the nose cap of the fuselage and 
projected eleven inches above the surface at about the 20° forward sweep 
angle. 

A probe using these principles was discussed in reference 12 by 
Frank Irving. Although a complete description did not accompany the data 
presented in the article, a photograph of a probe mounted on a sailplane 
fin indicated it to be a bent tube design. In the description, it was 
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stated that the design incorporated the 70° bend inclination to the 
airstream (20° forward sweep), the same as suggested in reference 1. The 
data show the probe as being relatively insensitive to incidence or angle 
of attack changes, although variations of 4 percent were indicated at about 
-10° of incidence. 

In the same article, data were presented for a so-called "Modified 
Nicks" probe which undercompensated; however, no details of the modification 
were given, so it is not possible to assess the reasons for the under- 
compensating pressure coefficients. Data from this report on the best 
probe tested by Irving are shown in figure 11 for comparison with the 
configurations reported in reference 1 and the bent-up probes reported here. 
In a personal communication, Irving described his probe as a bent probe 
configuration having two holes at 0 = i 130° at an X/D value of 1.5, 
with a Y/D of about 7. Based on these data, it appears that the probe 
tested by Irving has a greater variation during pitch changes than the 
bent probe versions having sensor lengths of 11 diameters or more as 
tested and discussed in this report. 

As a matter of general interest, it should be mentioned that probes 
made like configuration A of reference 1 have been tested in the NASA 
8-foot Transonic Pressure Tunnel over a Mach number range from M = 0.15 
to M = 0.75. Although test results have not been published, a variety 
of X/D values were tested confirming the X/D = 2 as most suitable for 
obtaining Cp values of -1.0. The preliminary data show coefficients 

within 5 percent of the desired value up to M = 0.3, with only slightly 
more variation to M = 0.75. The probes are being installed on a transport 
aircraft for experiments in wind shear detection and total energy management. 


DAMPING RESTRICTORS 


The separated flow region behind a cylinder operating at subcritical 
Reynolds numbers produces a fluctuating pressure which may couple with the 
dynamic characteristics of a sensitive variometer and cause needle 
oscillations or "vibrata" effects on an audio signal. Enlarging the orifice 
diameter or changing the tubing volume connecting the sensor to the 
variometer can affect these natural frequencies without modifying the 
average signal pressure, but a recommended solution to this effect involves 
the use of a damping restrictor-volume combination. Even if the natural 
frequencies of the sensor-variometer system do not cause oscillations, 
gustiness will produce fluctuations which tend to compromise the usefulness 
of total energy readings. 
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The simplest form of damping or gust filter can be made with a simple 
capillary restrictor and volume placed in the line between the total energy 
sensor and the variometer. Most mechanical variometers have time constants 
of about 2-7 seconds, although some electric variometers and a few 
mechanical models are capable of response times of less than a second 
(ref. 13). It is doubtful that response times of much less than two 
seconds are useful when gusty conditions exist; however, a gust filter 
will become the limiting factor if it has a slower response than the 
instrument, and the selection of a filter response rate must take this 
into account (ref. 14). 

A good discussion of gust filters and ways to make them appeared in 
reference 15. The time constant for a restrictor-capacitor pneumatic gust 
filter is given by the equation: 


T = (1.06 X 10"^) ^ 

PD^ 


(7) 


where: 

T = gust filter time constant, sec 

C = capacity of the gust filter flask, cubic inches 

L = length of capillary tube restrictor, inches 

P = absolute pressure of the atmosphere at expected flight 
levels, pounds per square inch 

D = internal diameter of capillary restrictor, inches 

Theoretically, it is desirable that the filter capacity volume be larger 
than the variometer flask when the variometer system uses such; however, 
a capacity equal that of the vario capacity seems to suffice. For electric 
varios having a built-in volume, a small capacitor may be suitable. For 
example, good results have been achieved with a fin mounted sensor having 
a built-in restrictor made of 0.508 mm (0.020 inch) inside diameter 
capillary 25.4 mm (1 inch) long, when used with an electric vario having 
a small internal volume. In this case, the length of tubing from the 
fin to the vario served as a capacitor. 

Using the equation above, restrictor lengths have been calculated 
and presented in figure 12 for two restrictor inside diameters as a 
function of capacitor volume, for time constants of 1.0 and 2.0 seconds, 
at an altitude of 1524 meters (5000 feet). From these relationships, it 
can be seen that restrictors made of 0.508 mm (0.020 in.) or 
0.794 mm (1/32-in.) inside diameter tubing are of practical lengths for 
a total energy system. Gust filtering is very important, and every total 
energy probe of the types described herein should have a restrictor added 
into the probe or installed in the line nearby. For fast response varios. 
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about 2.5 cm (1 inch) of 0.508 mm (0.020 inch), or about 8 cm (3 inch) of 

0.794 mm (1/32-inch) tubing should be helpful for gust filtering without 
introducing undesirable delays, even if the capacitor volume is small. 


FLOW FIELD EFFECTS ON TOTAL ENERGY SENSORS 


Obtaining the pressure coefficient for good compensation is best 
achieved if the sensor can be located in the freestream, unaffected by 
attitude changes of the aircraft. It is not necessary that the local 
pressure be the same absolute value as the freestream; it is only necessary 
that the static pressure, relative to freestream, not vary as the aircraft 
attitude changes. Because of this, a desirable sensor location must take 
into consideration local flow field changes during maneuvers. There are 
several aspects of flow fields which may be important: 

1. Boundary layer growth along the body 

2. Flow angularities caused by the windshield and the wing body 

intersection 

3. Downwash caused by lifting surfaces deflecting the flow 

4. Movable control surfaces which may propagate pressure influences 

upstream 

5. Induced velocities above wing or fuselage 

The boundary layer consideration is largely relevant if probes are located 
on the aft portion of the fuselage. Flows tend to parallel fuselage 
surfaces aft of the wing, so that a location roughly mid-way between the 
wing trailing edge and the tail offers relatively constant flow conditions 
for total energy sensors, provided that the sensor is located far enough 
from the body to avoid the boundary layer at all angles of attack or yaw. 
For aft fuselage mounting on the upper side, the sensor element should 
be located about 7 inches above the surface to insure avoidance of boundary 
layer fluctuations as attitude changes. 

Sensors have been located successfully on the noses of sailplanes; 
however, for this location there often are significant flow angularities 
as the flow streamlines are diverted around the body. Canopy bumps may 
cause local effects which would be undesirable, for example, and when 
positioning at the proper sweep angle, it should be recognized that stream- 
lines parallel the surface at the surface. 

High performance sailplanes usually achieve some laminar flow on the 
nose portion of the fuselage; a performance penalty may result with a probe 
in the laminar region which triggers an early transition from laminar to 
turbulent flow. This is not a problem for training sailplanes or others 
which do not depend on laminar flow for performance. Judgement must be 
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used in determining the proper sensor sweep angle on a curving surface, 
and experimenting with flight tests may be necessary. 

The vertical fin location usually offers near freestream conditions, 
provided the probe is positioned so that the fin, rudder, stabilizer and 
elevator (especially for Tee Tail configurations) are taken into account. 

The principal downwash in the flow field at the fin location is caused 
by the wing deflecting the air to produce lift. A series of calculations 
have been made to cover the range of effects for typical sailplanes during 
cruise and climb conditions. The downwash flow angle is a function of the 
lift coefficient being achieved at a given time. Reference 16 provides a 
thorough discussion of the mechanisms affecting the downwash as well as 
analytical methods for use in calculations. Based on these techniques, 
and the dimensions for short coupled sailplanes like the 1-26, the down- 
wash angle at the fin tip in degrees is about three times the lift 
coefficient, . In the cruise condition, the lift coefficient for the 

1-26 is about 0.5, making the downwash angle only 1.5 degrees. In the 
climb condition, the lift coefficient is about 1, making the downwash about 
3°. For high performance sailplanes having longer wings and fuselages, 
the downwash values decrease to about half those for a 1-26; that is, the 
range of downwash angles at the fin may be about 1.2 to 1.8 times 

degrees. The range of lift coefficients may be somewhat greater due to 
flaps; however, the total downwash variation for high performance sailplanes 
may still be less than 3°. 

For a fin installation, the sensor should be positioned at least 5 to 
10 times the maximum fin thickness ahead of the leading edge (ref. 17). 
Severe rudder deflections may cause significant lateral flow inclinations; 
however, the insensitivities of the simple probes described herein are a 
real advantage. Horizontal tail movements affect the downwash flow field 
to some extent. When attitude changes are being made, transients may be 
noticed; however, the effects can be minimized by smooth movements of the 
control surfaces. Sailplanes that are well balanced will not have very 
large tail lift coefficients, and therefore small downwash effects. 

In summary, a sensor location insensitive to changes in attitude is 
necessary for operation over a broad range of locations. Aft fuselage and 
vertical fin locations can be suitable for the probes discussed. Nose 
installations may be acceptable for low performance sailplanes; however, 
they must be positioned carefully. 
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CONCLUDING REMARKS 


Experimental pressure coefficients suitable for total energy 
compensation have been obtained using principles of laminar flow around 
a small inclined cylinder. To obtain the correct flow relationships, the 
sensor orifice should be located carefully with respect to the end of a 
3-dimensional cylinder; several options for providing the proper relation- 
ships have been extended by the current study of probes made of bent-up 
tubing. Total energy-pressure relationships have been reviewed to explain 
the principles involved and further explanations of 3-dimensional effects 
have been presented. 

In general, it has been shown that probe sensors with lengths as 
short as 7 times the outside diameter of the tubing used can be made to 
work with certain orifice locations. On the other hand, data have shown 
that sensitivities to manufacturing tolerance and flow incidence angles 
are reduced when sensor lengths of 11 diameters or greater are used. 

Comparative results from a number of experimenters have verified the 
principles and findings previously presented. The most significant of 
these probe dimensions are the sensor hole location geometry and the best 
angle of sweep for compensation that is insensitive to range of angles of 
incidence. 

Damping restrictors are useful to filter gusts and may be simply made 
by installing a small section of capillary tubing in or near the total energy 
probe, in series with an appropriate capacitor volume. 

Flow field effects around aircraft can affect the compensation of 
total energy sensors and must be considered. Among the effects are the 
boundary layer growth, flow angularities, downwash caused by lifting 
surfaces and movable control surfaces which may propagate pressure 
influences. The significance of these effects and ways of accounting for 
them are discussed. 
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PHYSICAL RELATIONSHIPS 
GLIDING FLIGHT 



POTENTIAL ENERGY. 

PE = HMg ; ALTITUDE xMASS x g 
KINETIC ENERGY. 

KE= ^ ; MASS X M LOC I JY)^ 

TOTAL ENERGY . 

TE = PE + KE = HMg + 

FIGURE 1 
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CALCULATED STREAMLINES 
PAST A CIRCULAR CYLINDER 

(REF. 5) 



FIGURE 2 
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TYPICAL PRESSURE DISTRIBUTIONS 

TWO-DIMENSIONAL CYLINDERS 
(REFS. 6 AND?) 



FIGURE 3 


237 


BENT-UP PROBE CONFIGURATIONS 


CONFIG. 

NO. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Y/D 

7 

9 

11 

13 

7 

7 

7 

7 

7 

11 

11 

11 

7 

X/D 

2.0 

2.0 

2.0 

2.0 

1.5 

1.75 

1.5 

1.75 

2.0 

1.5 

1.75 

2.0 

1.5 

HOLE 

POSITION 

e = 180° 

e = 135° 

9 = 
±135° 


X 



CONFIG. CONFIG. CONFIG. 
NOS. 1-6 NOS. 7-12 NO. 13 
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EFFECT OF PROBE LENGTH 
BENT-UP PROBE, X/D = 2.0 
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□ 9 

O 11 
O 13 


CONFIG. NO. 
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FIGURE 5 
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EFFECT OF HOLE POSITION 
SHORT BENT PROBE, Y/D = 7 



PRESSURE 

COEFFICIENT 

-C 

P 



X/D 
□ 1.5 
O 1.75 
A 2.0 


CONFIG. NO. 
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FIGURE 6 


PRESSURE DISTRIBUTIONS AROUND 
A TWO-DIMENSIONAL CYLINDER AT 
VARIOUS SWEEP ANGLES 

(REF. 8) 



FIGURE 7 


POSITION EFFECT OF ONE HOLE 
LOCATED AT 0 = 135°; Y/D = 7 
BENT-UP PROBE 




X/D 
O 1.5 
□ 1.75 
O2.0 


CONFIG. NO. 
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8 
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to 


FIGURE 8 
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PRESSURE 

COEFFICIENT 


POSITION EFFECT OF ONE HOLE 
LOCATED AT5 = I35“, Y/D = II 


BENT-UP PROBE 


FLOW 

► 




135 ° 



r , FORWARD SWEEP, deg 


CONFIG. NO. 
10 
11 
12 


FIGURE 9 
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COMPARISON OF ONE HOLE VERSUS TWO HOLES, 
0= 135® Y/D = 7, X/D = 1.5 

BENT-UP PROBE 


135° ±135° 



PRESSURE 

COEFFICIENT 



FIGURE 10 
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COMPARISON OF CONFIGURATIONS 
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FIGURE 11 
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RESTRICTOR-CAPACITOR RELATIONSHIPS 
FOR TWO TIME CONSTANTS 
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VOLUME, 

liters 


PROBE 



0 2 4 6 8 10 12 14 16 18 20 


RESTRICTOR LENGTH, cm 
FIGURE 12 
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HOW ACCURATE IS NETTO? 


Stephen du Pont 
Soaring Society of America 

SUMMARY 


The historical origin and general history of the MacCready vertical 
current total energy variometer (now termed "netto") , including its optimum 
airspeed selector ring are reviewed, and some later developments of it are 
discussed. Polars of three sailplanes of different spans are charted for 
straight and circling flight, then plotted to reveal their parabolic anomaly 
and the effect of circling flight sink rate. These effects are further 
analyzed for their influence on the transient compensation of netto variometers 
as well as the speed ring. Some other disturbances due to the quality of 
sailplane preparation and flight dynamics are listed. Conclusions are drawn 
about the problems to pilots from imperfect netto variometer compensation 
and its effect on the maximization of ground speed from the speed ring. 

A modification for improvements to the speed ring and computer is 
suggested. 


DISCUSSION 


Ideally a variometer would be "compensated" to eliminate needle deflec- 
tion resulting from speed changes so as to show only the vertical motion due 
to vertical air currents surrounding the glider. But the sailplane is an 
isolated energy system in which changing speed requires an exchange of its 
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potential energy for kinetic energy, causing changes of altitude. Its drag 
forces also increase generally with airspeed squared, producing a parabolic 
variation of sink rate with airspeed (ref. 1). 

The pure total energy variometers are compensated through the change of 
dynamic pressure with airspeed which is converted by either a venturi (ref. 2), 
Braunschweig tube, etc., or an elastic bellows (ref. 3) to cancel the unwanted 
"stick thermal" coming from the zoom and dive maneuvers necessary to change 
speed in a glider. But they do not produce a fixed variometer needle with 
airspeed variations because of the drag force variation. Thus they leave some- 
thing to be desired. 

In 1949, Paul MacCready Jr. had disclosed his invention of the "speed 

ring" in a paper (ref. 4) read at the lAS-SSA meeting at Elmira while its 

author was busy winning the National Soaring Championship there. In 1954 

MacCready first disclosed at the lAS-SSA meeting in New York his new "vertical 

current" variometer (ref. 1) today known as "netto" which more effectively 

than anything else known even today improved the "compensation" of sailplane 

variometers. At the same time it greatly simplified the use of the speed 

ring. The theory of the device was to leak a small calibrated flow of air 
2 

proportional to V outwards through the variometer, causing it to add the 
sink rate for still air to the variometer indication. The result was to 
indicate the vertical air current rather than the climb or sink rate of the 
sailplane. MacCready pointed out that with this arrangement the speed ring 
now indicated directly the speed to fly, dispensing with "iteration", that 
is the need to chase the needle while bringing the sailplane speed to 
a number never quite stabilized on the ring by the variometer needle. Today 
it seems incomprehensible that the soaring world took twenty years to appreciate 
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the device which twenty five years later remains the best. It is perhaps 

natural that a meteorologist would be preoccupied with the measurements 

of vertical currents, but only one who was also a sailplane expert would have 

combined the need with the tool and come up with such an ultimate solution. 

There are several ways of plumbing the laminar leak (refs. 1, 3, and 5) 

used for the vertical current compensator to various total energy systems so 

2 

as to produce the specified flow proportional to V . MacCready had placed 

the leak between the pitot and the variometer capacity line since he had used 

a venturi type total energy compensator (ref. 2) which is connected outboard 

of the variometer static vent. This produces a pressure differential across 

2 

the leak of twice the dynamic pressure, which we call 2q, where q = 1/2 pV , p 

being density. Thus, + q - (q + P^) = 2q, where P^ is the static pressure 

2 

of the altitude of the sailplane. Because q contains V the pressure across 

2 

the leak is proportional to V and the flow through a laminar leak is pro- 
portional to the pressure across it. 

In "Soaring" for 1975 (ref. 5) Don Ott's arrangement of the leak with a 
total energy venturi was described. He ran the netto leak from the capacity 
line to static pressure, and showed that cockpit pressure was close enough, 
so he left the leak simply open to the cockpit. Here the pressure across 
the leak is: 

p, - (-q + P,) - q 

2 

still giving the specified V pressure variation of MacCready. 

Where a bellows or diaphragm total energy compensator (ref, 3) is used 
which is driven by the difference of pitot pressure and the capacity (Burton, 
PZL, Schuemann etc.), the leak parallels the bellows and the pressure 
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differential is: 


and again we have MacCready's V pressure variation. 

With the arrangement shown by MacCready, the capillary must be twice as 
long as with the other systems mentioned here. Reference to the 1975 Soaring 
article (ref. 5) will yield the formula for calculating the capillary. 

In 1954 the Ka-6 sailplane was three years in the future. But many 
modern sailplanes today incorporate camber changing flaps, broad drag bucket 
airfoils, and very rigid and smooth wing skins, and realize broad areas of 
laminar flow in the boundary layer. Still the shape of the polar is much as 
MacCready had described it: "approximately" and "fairly exactly" and "within 

a few inches per second" of parabolic, which means of course that the sink rate 
is nearly proportional to airspeed squared. It is nearly proportional, but 
not quite, as figures 4 to 6 herein show. Those are generally sharper curved 
at the low end than the parabola, that is the sink rate in that area is 
greater. The PIK does not show any droop at the fast end while 1-26 and the 
AS-W 17 do show it. 


SCOPE AND LIMITS OF STUDY 


We investigate here specifically a) those netto errors, transient and 
steady, that come from the parabolic anomalies of some typical real world 
polars, b) those netto errors that are due to the variations from these 
straight flight polars that occur in banked circling flight, and c) errors 
from the netto speed ring. We do not go into some other errors due to 
1) Flight off design altitude of variometers. 
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2) Transient variations in L/D from maneuvers causing accelerations that 
change the wing loading and from L/D changes caused by uncoordinated skidding or 
slipping. 

3) Changes in L/D from wing loading variations due to weight of crew, 
equipment, ballast, center-of-gravity effect, etc. 

4) L/D variations from flap position, aileron rig, aileron flap couplers, 

etc. 

5) L/D variations coming from air leaks, ventilators, loose fairings 
and gear doors, etc. 

6) Pitot static position and airspeed calibration errors. 

7) Total energy probe and vent position errors from wing pressure field 
and wake, probe yaw and angle of attack errors, and uncoordinated skidding 
and slipping probe errors. 

8) Plumbing hose length and capacity resistance in fittings, and hose 
pinch, restrictors, filters, electric damping, etc. as they affect variometer 
indication. 

9) Dirty wings and skins, skin stress wrinkles, assymetrical ballast, 

etc. 

Such errors may or may not be transient, will be difficult to predict, 
detect, and measure, may be cumulative, and may have comparative values that 
are significant compared to the errors analyzed. 

In the light of the foregoing a knowledgeable pilot has a right to 
wonder about the reliability of his netto for finding better air as well 
as for the maximization of ground speed from the speed ring or computer. 

He may ask 
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1) How important is the effect of the parabolic anomaly of polars versus ' 

2 

the true parabolic characteristics of the V driven leak? 

2) If the variometer that is read in circling flight, and the speed ring 

set to it, is compensated according to a polar derived in straight flight I 

(since this indicates too low a climb and sets the ring too low (slow), will 
this significantly affect the maximization of ground speed? 

3) If the speed ring that is derived for circling climb (no progress 
being made along course while climbing) is set during straightaway climb 
along course, this sets the ring too low (slow) (refs. 6 and 7). Does this 
affect the maximization of ground speed? 

4) What are the effects of a mismatched leak to accomodate parabolic 
anomaly (it would be too short, too much flow) or one for the wrong sailplane 
or incorrect polar? (Capillary somewhat too short or too long.) 

5) Is netto well enough compensated both for climb circling or straightawa] 
so that it can be relied upon during small speed changes? 

To obtain answers to these questions we have considered polars of the 12- 
meter 1-26 sailplane, the 15-meter PIK-20 and the 20-meter AS-W 17 (ref. 7). 

The polars have then been reworked by recalculating the slow ends for 40 and 50 

I 

degrees of banked turning (30 degrees has an insignificant effect) by the method 
of the Appendix (taken from ref. 8). (See figs. 1 to 3.) 

The results, including the noncircling polar, have been replotted with the 
sink rate against airspeed squared. When any parabola is so plotted, it is a 
straight line passing through the origin. This makes it easy to inspect 
the parabolic anomalies of the curves. True parabolic characteristics of 
netto leaks can be compared to polars. If the polar were parabolic and the leak 
matched it they would plot the same, as one straight line. If the leak 
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is enlarged to compensate for some droop from the parabolism of the polar, 
the leak plot will be lower, cutting across the tighter curve of the polar 
(figs. 4 to 6). The leak sink rate in the midspeeds will then be lower and 
higher at the end portions. The difference, or the extent of the anomaly, 
can thus be directly measured on the plot along the sink direction under 
the square of the airspeed of interest, from any straight line that passes 
through the origin. 

The effects on netto speed rings of the parabolic anomalies and of 
circling sink rates, as well as climbing along course with a ring designed for 
circling climb, have been computed for climb rates of 200 and 800 ft/min. 

The classic graphic analysis is used before and after modifying the polar 
values to reflect the above mentioned deviations. The work was done directly 
on the preplotted polar analysis sheets from reference 7 and is not reproduced 
here. Instead, the results are shown in table 1. There they can be seen to 
have an insignificant effect on the maximization of ground speed by the speed 
ring. 


TRANSIENT ERRORS OF NETTO COMPENSATORS 


Even small parabolic anomalies of the polars cause annoying transient 
netto compensation errors whose trends and speed ranges can be seen in figures 
4 to 6, and are discussed below and further detailed in table 2. 

The time lags of variometers (they vary widely between models and installa- 
tions) (ref. 9) will modify these transient errors. In the slow portion of the 
polars (see figs. 4 to 6), slowing down causes an erroneous indication of 
worsening air (decreasing rise, increasing sink). Variometer time lag tends 
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to reduce this error. The effect begins as slight at about 55 knots and 
becomes larger with decreasing speed. This is an important regime of soaring 
flight where after the pilot has had a signal of better air he may be slowing 
down in an effort to locate the lift and use it, and possibly to circle and 
center it. If circling has begun during the slowing, the compensation errors 
are increased. Here the pilot should stay alert to his sense of vertical 
acceleration, heed the roll motion of his wing tips, and pay less attention to 
the netto variometer. As the charts show, the effects become significant 
with the 1-26 below 50 knots, the PIK -20 below 55 knots, and the AS-W 17 
below 55 knots. 

In the slow zones of the polars, speeding up will indicate erroneously 
better air than actual and the errors are increased by variometer time lag. 

But in this flight regime, the netto and its compensation become less important 
in the search for lift, because a pilot will probably be accellerating into 
his speed ring glide towards the next gaggle or cloud. 

At the fast zones of the polar, slowing down will erroneously indicate 
improving air and variometer time lag will worsen the error. This is an 
important regime of soaring. The effect lessens as the speeds are lowered, 
vanishing with the 1-26 above 75 knots and with the AS-W 17 at above 85 knots. 
The PIK-20 is not affected. 

In the fast portion of the polars, speeding up shows better air than actual 
and variometer lag will decrease the error. Searching for lift is less 
important in this regime, making the error less disruptive. The error increases 
with speed above the speeds noted in the previous example. Table 2 shows some 
values for smaller speed changes and steepened banks. 
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Pilots may check the calibration of their netto by flying in still air 


at three or preferably five speeds in the slow, medium, and fast portions 
of the polar. At the steady speed in still air the netto should indicate 
zero. The non netto total energy variometer should indicate the sink rate 
for each speed tested taken from reliable polars such as those derived by 
Bikle and Johnson in various issues of Soaring and those found in reference 7. 
Individual airspeed calibration errors may be troublesome here. It is 
convenient to tabulate the speeds and sink rates of interest for use in the 
cockpit. Indicated airspeeds are used. 

Errors of transient compensation of netto can be noted in the slow and 
faster ends of polars by slowing down and speeding up while noting the behavior 
of the needle as it deviates from zero. Again the air must be still for this. 
What the pilot observes here in still air, he can apply to his actual soaring. 


ERRORS OF THE NETTO SPEED RING 


Calculations of speed rings and related computers (as well as of netto 
leak sxze) have usually been based upon the straight flight polars (refs. 1, 
4, 9, and 10), but the ring is often set during banked circling climb where 
the actual vertical sink rate is greater than the netto variometer is indi- 
cating. This causes uhe ring or computer to be set too low (slow). Rings 
based on the pure MacCready mode (fig. 7) consider that no progress is being 
made along course while in circling climb (ref. 4). The so called street 
speed ring (refs. 6 and 7) of figure 8 acknowledges that some course distance 
is used up during straightaway climb, thus shortening the distance to go to 
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the next lift and so allowing a steeper glide at faster speed. The pure MacCready 
mode ring would be set too low (slow) while in climb along course. A street 
speed ring would seem better when set for both circling and along-course climb 
netto. 


CONCLUSIONS 


1) A Pilot who flys glued to his variometer may well be confused by the 
compensation errors of netto in slow straightaway climb and in steeply banked 
climb if his speed is not steady. 

2) In slow speed climbing by netto, holding speed steady or use of a 
standby non netto total energy variometer may be helpful. 

3) While shortening (increasing the flow of) the netto leak may give 
better steady speed matching to the anomalies of the polar, this will not help 
transient compensation errors. 

4) The^ time lag of variometers has an effect on the transient compensation 
of netto. Where this is favorable speeding up, it will be unfavorable slowing 
down, and vice versa. 

5) The netto errors studied have a negligible effect on the maximization 
of ground speed from the speed ring. 

6) The match of the netto speed ring or computer might be very slightly 
better through use of the street ring construction. 

7) The study hints that proponents of steady speeds instead of a miriad 
of little speed changes may have a point, due, if for no other reason, to the 
fact that even the best variometer system is plagued with compensation errors. 

A slow response variometer is favored by some to delude themselves to 
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thinking the variometer is compensated correctly, while fast repsonse is 
the goal of certain variometer makers who see better compensation giving 
more accurate information. 

8) One total energy variometer with netto and a total energy vario- 
meter without netto might be the best solution for most of us, using the 
one that fits the soaring situation of the moment. 
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APPENDIX 


Determination of sink rate due to circling flight for figures 4 to 6 
is done by the method treated in "The New Soaring Pilot" by Welch and Irving. 
There it is assumed that the sailplane has the same L/D so long as the circling 
flight angle of attack is the same as in straight flight. The formulas are 


= V 

4> o 


/ 


cos , 


( 1 ) 


and 


S, = S 


cos , 


/ 


cos 


( 2 ) 


where the subscripts 

(j) and 0 are bank angles, V is airspeed and S is sink rate. 
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TABLE 1 

MAXIMIZATION OF GROUND SPEED AS AFFECTED BY PARABOLIC ANOMALY AND CIRCLING SINK RATES WITH NETTO SPEED 
RING AS COMPARED TO SOLUTION FROM CLASSIC MacCREADY GRAPH. 


Glider 

types 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


Actual 
glider 
climb , 
fpm 

Still air 
sink 

best sink 
speed, 

bank A0° 

Polar sink 
deviation 
from para- 
bolic, fpm 

bank 40° 

Actual 
airmass 
rise, fpm 
(l)+(2) 

MacCready 
speed to 
fly, knots 

Netto 
indicated 
airmass 
rise at 
zero bank 
(l)+(2)+(3) 

Netto 

ring 

indicated 
speed to 
fly, knots 

MacCready 
ground 
speed , 
knots 

Netto 

ring 

ground 

speed 

knots 

1-26 


40° at 40 
knots 

40° at 40 
knots 








200 

275 

-80 

475 

53 

555 

55 

22.5 

22.5 


800 

275 

-80 

1075 

71 

1440 

73.5 

43.5 

44 

PIK-20 


40° at 51.5 

40° at 51.5 








200 

195 

-20 

495 

70 

415 

70 

32 

32 


800 

195 

-20 

995 

92 

1015 

87 

59 

59 

AS-W 17 


40° at 52 

40° at 52 








200 

160 

-50 

360 

64 

410 

75 

38.5 

36 


800 

160 

-50 

960 

105 

1010 

97 

66 

65 


TABLE 2 


SOME TRANSIENT COMPENSATION ERRORS OF NETTO VARIOMETER FOR THREE 
SAILPLANES, NOT CONSIDERING VARIOMETER LAG 



Slowing 
5 knots, 
bank angle 

50° 

Slowing 
5 knots, 
bank angle 
zero 

Slowing 
5 knots, 
bank angle 
zero 

1-26 




A IAS 

40 to 35 

40 to 35 

43 to 33 

As ink 

- 15 fpm 

- 30 fpm 

- 50 fpm 

at sink 

- 350 fpm 

- 180 fpm 

- 210 fpm 

PIK-20 




AIAS 

60 to 55 

40 to 35 

55 to 45 

Asink 

-30 

-10 

-20 

at sink 

-250 

-130 

-130 

AS-W 17 




AIAS 

60 to 55 

50 to 45 

55 to 45 

Asink 

-25 

-20 

-20 

at sink 

-200 

-120 

-110 
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IAS , EAS , knots 


Still air 
sink, knots x 100 
and ft/min 



Figure 1. - Performance polars ; 1-26, 12 meter. 


IAS , EAS , knots 
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Figure 2. - Performance polars: PIK-20, 15 meter. 


IAS , EAS , knots 


Still 

sink. 


air 

ft/min 



Figure 3. - Performance polars : AS-W 17 , 20 meter. 


IAS, EAS, knots 


Still air 
sink, ft/min 


UO 50 60 TO 80 90 



2 

Figure U. - Sink versus V : 1-26, 12 meter. 
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IAS , EAS , knots 


Sink, 
ft /min 



2 

Figure 5- - Sink versus V (see fig. U): PIK-20, 15 meter. 



2 

Figure 6. - Sink versus V (see fig. U): AS-W 17, 20 meter. 



A Circling climb at zero speed 
along course . 

B Average groundspeed, for one 
climb glide cycle, begin and 
end at same altitude. 

C Tangent at polar. 

D MacCready optimum airspeed 

between climb zones . 

EC Ring layout arc, ft/min between 
index and speed to fly. 

FG Tolerance of speed to fly for 
insignificant effect to average 
ground speed. 


Figure 7- - Pure MacCready speed ring 

derivation, all climbs circling. 



A Rate of climb along course at 
speed H 

B Average ground speed for one 
climb glide cycle, begin and 
end at same altitude. 

C Tangent at polar (see above FG) 

D Ring optimum airspeed between 
climb zones. 

EC Ring layout arc, ft/min between 
index and speed to fly. 

H (I) flight speed during climb 
along course . 

I Climb speed 


Figure 8. - Street speed ring derivation, 
climbs are along course . 
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THE APPLICATION OF MICROPROCESSOR TECHNOLOGY TO IN-FLIGHT COMPUTATIONS 


Patricia L. Sawyer and Dan M. Somers 
Hampton, Virginia 


INTRODUCTION 


During the past decade, electronic microcomponent technology has 
developed to the point where microprocessors are appearing in an expanding 
number of diverse applications. Improved manufacturing techniques have 
allowed the computing power of the microprocessor to increase at the same 
time its physical size has decreased. The mass production of these tiny 
high-powered components has made them readily available for nominal costs to 
industry, research, government, and private individuals. 

The application of this new technology proposed by the authors is a 
microprocessor-based system to essentially replace current panel instrumenta- 
tion in sailplanes. A microprocessor-based system is particularly attractive 
for in-flight computations because of its accuracy and dependability, but 
more importantly because of its potential for flexibility of output. Given 
the two basic input values of dynamic and static pressures, the system could 
perform the functions of a conventional panel by providing the pilot with 
airspeed, altitude, and variometer readings, but could also be programmed to 
provide him with additional data at his request. Since all microprocessor- 
based systems contain an internal clock, time history data and moving averages 
could be computed. The microprocessor could easily handle computations in- 
volving integrations or requiring the manipulation of many data points or 
solutions of complex equations. 

This paper will present a modular design of a general purpose micro- 
processor-based computer to perform in-flight computations for cross-country 
soaring pilots. The basic requirements for the system will be discussed. 
Several specialized applications of the computer will be presented, including 
real-time pilot feedback and flight-test data acquisition and reduction. 


CURRENT INSTRUMENTATION 


A brief examination of the instrumentation in current, competition sail- 
planes will reveal part of the motivation behind the microprocessor-based 
panel. An excellent overview of present-day instrumentation is contained in 
reference 1 . 

Today's airspeed indicators and altimeters, although simple and reliable, 
suffer from position errors caused by incorrect static-pressure sources. In 
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the case of the airspeed indicator, the relocation of static ports in the 
fuselage between the wing and the tail has eliminated the major portion of 
this error. The remaining correction can still be significant, however, 
especially when performing flight tests. The static-pressure source for the 
altimeter is usually cockpit pressure, for reasons which will be discussed 
later. In addition to this position error, altimeters have added problems 
caused by high relative friction forces within the instrument. This friction 
results in the sticking problems known to all sailplane pilots. 

Present-day variometers are of two types — those which measure pressure 
differences and those which measure mass-flows. The pressure variometers have 
no altitude errors, whereas the mass-flow variometers exhibit a rather large 
altitude error, approximately 10%/3000 ft. 

Of the four available pressure variometers, two are mechanical and two 
are electrical. The first, and presently most common, mechanical type is the 
classic vane variometer such as the PZL or Winter. This type is quite re- 
liable, but has several disadvantages which include relatively slow response 
due to friction, close-tolerance construction, delicate mechanisms, and sus- 
ceptibility to dirt. The second mechanical type, the taut-band variometer 
produced by Schuemann or Bohl i , has a fast response. It also has some inherent 
problems, however, associated with tuning the required capillary and producing 
the delicate metal -band assembly. 

The two electric types have the common virtues of fast response and audio 
capability. They also share the disadvantages of possible unreliability 
because of the required power source and wiring, and a response which is "too 
fast" and thus requires damping. The most common of this type, the capillary- 
leak variometer made by Ball, has further difficulties arising from the 
capillary which must be tuned and the close tolerances on the linear variable 
differential transformer (LVDT). The second type is the altitude derivative 
variometer (ref. 2). This variometer is expensive to produce because it 
requires a pressure transducer. 

Mass-flow variometers, such as those made by Cambridge or Moore, incorpor- 
ate thermistors which require some degree of matching. 

All of the variometers mentioned above require total -energy and/or netto 
compensation. Total -energy compensation is usually provided by either an 
Althaus-type venturi or a diaphragm. The venturi, as well as similar devices 
such as Braunschwieg, Nicks, and Irving tubes, provides accurate compensation 
independent of altitude. They are, however, relatively sensitive to horizontal 
gusts, require careful positioning on the sailplane, and add drag. They also 
do not work well, if at all, in rain. A diaphragm compensator, such as the 
PZL, Winter, or Burton, does not increase the drag of the sailplane. It does, 
on the other hand, compensate exactly at only one altitude. Although less 
sensitive to gusts than a venturi, it is difficult, if not impossible, to 
correctly adjust. Recently, electronic compensation has been introduced, 
which incorporates thermistors like those in the electric mass-flow variometer. 
These thermistors require some matching. 
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In addition to total -energy compensation, netto compensation is also 
desirable. This compensation requires a capillary which can be tuned to 
approximate, but not necessarily exactly "fit," the given sailplane's polar. 
With this system, it is difficult to account for weight changes such as the 
addition of water ballast. 

Regardless of type, all current instrumentation suffers from several 
problems. None of today's instruments are truly temperature-compensated. 
Because the entire system is pneumatic, the various pressure sources and 
instruments must be connected by tubing. This tubing adds volume to the system 
and introduces the possibility of leaks. Refinements to the variometer sub- 
system require additional restrictors and/or gust filters. While solving some 
problems, these devices have problems of their own, such as time-constant 
matching (ref. 3). Also, the interaction between the various instruments can 
cause difficult problems. These problems arise from the interaction of several 
diaphragms (airspeed, total -energy , and altimeter) through the variometer on 
the common static-pressure line. Disconnecting the altimeter from the static- 
pressure line alleviates part of this problem. 

In summary, the instruments currently available are certainly adequate 
for safety, but are somewhat inadequate for efficient cross-country flight 
for the following reasons: 

(1) The information they provide is not accurate for all altitudes, 
attitudes, and environments. 

(2) To correct some of the errors inherent in the instrumentation, the 
devices for measurement and display have become inordinately delicate and 
compl icated. 

(3) The current equipment offers no facility for flexibility or expansion 
of the type of information provided. 


EARLY ATTEMPTS AT IN-FLIGHT COMPUTERS 


During the past few years, several attempts have been made to expand 
the amount and vary the type of information provided to the pilot by his 
instruments. In 1972, A. Raouf Ismail, of Cambridge Instruments, suggested 
that an analog air-data-computer (ADC) could provide the pilot with average 
rate-of-cl imb and speed-to-fly information during flight, but cautioned that 
this approach could be bulky, power consuming, complicated, and expensive 
(ref, 4). That same year, John Firth, of Ottawa, Canada, reported on his work, 
begun in 1969, involving the coupling of an electric variometer to a simple 
analog computer (ref. 5), The system performed normal variometer functions 
with fully electronic total -energy compensation. In addition, the system 
computed an average rate-of-cl imb over variable time periods and displayed the 
correct speed-to-fly using a nose-up/nose-down pitch indicator. Variable 
damping and sensitivity, variable time constants, and a dry/wet switch allowed 
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more variables to be adjusted during the computations, and hence provided the 
input necessary for a more accurate output. 

During an average cross-country flight, the analog system was found to 
be quite reliable, and its flexibility enhanced the usefulness of the informa- 
tion provided. During the World Championships in Yugoslavia in 1972, however, 
the instrument was used very little because the soaring conditions did not 
match those for which the instrument was developed, the battery power on board 
had to be conserved for the radio and gyros, and the pilot's constant attention 
was required for navigation. 

The analog air-data -computer was a definite improvement over standard 
instrumentation, but was still greatly limited in its flexibility, and, 
therefore, in its accuracy. In 1974, William Foley presented results of a 
study involving the use of mini-calculators in soaring flight (ref. 6). 

Foley's use of a programmable digital calculator partially solved the problem 
of inadequate flexibility encountered earlier, but introduced other problems 
which seemed insurmountable given the current technology. 

The calculator was not connected to the instruments; thus, all data 
necessary for program execution had to be entered manually by the pilot. This 
visual -manual interface was not only slow, but held potential for error. The 
calculator itself was an obstacle in the cockpit. In order to enter data, a 
cumbersome kneeboard was required, and to read the LED output, the calculator 
had to be carefully shaded from sunlight. Depletion of battery power was a 
concern, and limited memory within the calculator demanded that the pilot load 
programs from magnetic cards during flight. 

Because of these inconveniences, Foley described his mini -calculator as 
interesting, but impractical. He recognized that the first requirement of any 
instrument system must be that it provide better or more timely information 
than the pilot can obtain using his own resources. 


REQUIREMENTS FOR AN IN-FLIGHT COMPUTER 


The following requirements can be used to define an in-flight computer 
with sufficient capability to insure that this important condition of faster 
and better information is met. 


Output Capability 

1. The in-flight computer must display its output during flight in a 
format easily readable by the pilot. For example, if the pilot is accustomed 
to a needle-and-scale altitude instrument, then the in-flight computer must 
not force the pilot to read a digital display. 

2. The system must provide, as a minimum, all information currently 
provided by standard panel instrumentation. For example, altimeter, airspeed 
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indicator, and variometer functions must be performed with all types of 
compensation, filtering, and signal conditioning. 

3. The system must provide audio capability. 

4. The system must be capable of digital, alphanumeric, and graphic 
display via an easily readable, glare-free output device. 

5. The system must be compatible with other input/output devices, such 
as flexible disks or printers, for pre-flight system development or post- 
flight data reduction. 


Computational Capability 

1. The in-flight computer must be driven by a system clock of sufficient 
speed to provide real-time information to the pilot. Real-time, in this 
context, means that the pilot must detect no hesitancy or flicker in the 
output displays. 

2. The system must be more accurate in its computations than standard 
current panel instrumentation. 

3. The system must provide hardware and software signal conditioning 
capabil ity. 

4. The system must be reliable and must provide a self-test capability 
to verify its integrity during pre-flight checkout. 

5. The system must include adequate memory so that time history data 
may be gathered and stored for post-flight analysis. 

6. The system must be built around a microprocessor of sufficient 
sophistication so that the flexibility that programming allows is not over- 
shadowed by the complexity of the actual programming process. For example, 
the system must support high-level languages and feature an operating system 
that is human -oriented. 

7. The system must accept standard programmable read-only memory (PROM) 
so that standard software can be transported and incorporated easily. 


Input Capability 

1. The in-flight computer must contain an interface to a standard 
alphanumeric terminal to be used for program development. 

2. The system must allow the pilot to initialize certain predetermined 
parameters during pre-flight checkout and also during the flight. For 
example, the pilot must be able to zero his altitude gage before takeoff if 
he desires. 
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3. The system must feature an interrupt system that will allow the 
pilot to stop the processing of one program and begin the execution of 
another in flight with a minimum of effort. 

4. Through the use of default parameters, the system must require a 
minimum of data entry during flight. 


System Configuration 

1. The processing and memory modules of the system must be portable 
and easily removed from the sailplane for post-flight analysis or program 
devel opment. 

2. The entire system must be small and lightweight, and must not be an 
obstacle in the cockpit. 

3. The system must require a minimum of power, and this power supply 
should be separate from that required by the radio. 

4. The system must be inexpensive. 


A MICROPROCESSOR-BASED IN-FLIGHT COMPUTER 


The proposed microprocessor-based in-flight computer system consists 
of three main subsystems for pressure sensing, computations, and input/output 
of information. A functional block diagram showing the major tasks of each 
subsystem is presented in figure 1. 


Pressure Sensing Subsystem 

The major components of the pressure sensing subsystem are two pressure 
transducers — one to measure dynamic pressure and the other to measure static 
pressure. Transducers with adequate accuracy and reliability are very small 
and can easily be located at the point of measurement, thus eliminating the 
need for tubing and its related problems. The analog signals representing 
dynamic and static pressures must be conditioned either by circuitry within 
the pressure transducer, or by the hardware or software of the computations 
subsystem. A photograph of a typical pressure transducer with adequate 
sensitivity for sailplane applications is shown in figure 2. 

Using the values of static and dynamic pressures, altitude and airspeed 
can easily be computed. Uncompensated and compensated rates-of-cl imb are 
merely functions of these pressures with respect to time. Because the system 
contains a memory and has access to stored and real-time information, the 
values it computes can be corrected for changing parameters, such as altitude 
or location on the polar curve. This flexibility will result in an output 
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that is far more accurate than that available from current instrumentation. 


Computations Subsystem 

The computations subsystem contains several important components which 
enable it to handle all of the processing functions of the system. The 
components are of two main types -- hardware, the actual circuitry that 
defines the capabilities of the system; and software, the programs that 
define how the hardware will be exercised to meet the needs of the applica- 
tion. The "brain" of the system is the microprocessor, the selection of 
which must be guided by several requirements. The microprocessor must 
operate at a speed sufficient to provide real-time information to the pilot, 
yet, at the same time, require a minimum of power. It must possess a rapid 
interrupt-handling feature to enable the pilot to redirect computations 
instantaneously. The processor word size should be of sufficient length to 
allow for adequate accuracy during the computations and the addressing of a 
sizeable memory block. Finally, the processor should be a standard, easily 
obtainable part that is well -supported by development and checkout software, 
such as high-level languages and debug packages. 

A microprocessor, such as the one shown with two memory units in 
figure 3, would satisfy these requirements. A typical 16-bit microprocessor 
will allow the direct addressing of over 65,000 words of memory and the in- 
direct addressing of many more words. The instruction set is adequate for 
this application and contains special instructions that allow the rapid 
processing of interrupts. Several 16-bit microprocessors already boast 
large user communities, are readily available, and are relatively inexpensive. 
At this writing, a typical unit price is under $25.00. 

I ' t 

The microprocessor cannot function as a computer without supporting 
components. The memory requirements of the system will depend on the type 
and number of programs to be executed. The amount of memory can be expanded 
at any time up to a point depending upon the addressing scheme of the 
processor and the space allotted for memory on the processor board. Memory 
in the system will be of two types -- programmable read-only memory (PROM) 
containing the library of predefined programs, and random access memory (RAM) 
which may contain instructions or data. The cost of memory represents a 
large portion of the cost of the computations subsystem, and the memory will 
also occupy a large portion of the processor board. The cost of memory is 
rapidly decreasing, however, even though the capacity of each memory unit is 
increasing. At this writing, 1000 words of memory, packaged as one unit, 
cost about $10.00. 

The processor board will also contain anal og-to-digital hardware and 
some signal conditioning circuitry, depending upon how much conditioning is 
performed by the pressure transducers themselves, or in the software of the 
system. 

The board may also contain special purpose hardware to perform 
repeatedly required functions more quickly than software can perform them. 
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For example, if the applications program involves the manipulation of values of 
very large or small magnitudes, then these values must be stored in memory as 
floating point numbers represented by an exponent and a mantissa. Most micro- 
processors available today perform computations involving floating point num- 
bers much more slowly than those involving fixed point numbers. For this 
reason, hardware has been developed to perform these floating point computa- 
tions very rapidly outside the microprocessor. The computations subsystem may 
operate more efficiently if specialized hardware to meet specific needs is 
included. 

The system will contain a synchronizing master clock and whatever 
circuitry is required to drive output gages, display screens, or audio 
devices, and to allow input from switches or keyboards. The entire system 
will be driven by battery power, which must be regulated by special 
circuitry on the board. Figure 4 shows a processor board built for another 
application which illustrates the arrangement of a 16-bit processor, memory, 
clock, special purpose floating point hardware, and input/output circuitry 
on a single board. 

All processing will be incorporated into a modular computations sub- 
system, possibly consisting of several boards that may be mounted near the 
panel area of the cockpit or elsewhere in the sailplane. Because the 
communication between each of the three subsystems is electrical, not 
mechanical, the location of the computations subsystem is arbitrary. 


Input/Output Subsystem 

System input and output will be of two types -- real-time and pre- or 
post-flight. The format of real-time output must meet the pilot's first 
requirement that information be timely and easily understood. Since 
standard needl e-and-scal e gages inform the pilot of absolute values, as well 
as the rate of change in these values, then this type of output display 
should continue to be used at least for airspeed and rate-of-cl imb informa- 
tion. Additional output desired by the pilot, such as speed-to-fly or 
altitude, may be displayed using a digital output device or a screen display. 
The choice of a proper screen display is important. A standard raster scan 
device, such as a small television screen, while easy to interface to the 
system, is heavy, cumbersome, requires a great deal of power, and does not 
provide adequate contrast in the cockpit environment. A liquid crystal 
display, however, uses ambient light to increase contrast and, therefore, is 
well-suited to cockpit applications. The digital display shown in figure 5(a) 
and the top view shown in figure 5(b) illustrate the small size of the liquid 
crystal display. These devices are also capable of graphic output, as shown 
in figure 6. A typical panel configuration, including an LCD display, is 
shown in figure 7. 

Input to the system during flight must be minimal and uncomplicated. 
Select switches and a small numeric/command keyboard must be located in the 
cockpit where they may be accessed easily, but will not inhibit flying. 


274 


Pre- and post-flight input and output can be more complicated and time- 
consuming than real-time, but should be streamlined by a carefully designed 
human-oriented system. A functional grouping of the input and output devices 
would spare the pilot the task of selecting one of an array of identical 
pushbuttons while flying. 

The entire computations subsystem is modular and portable so that it can 
be removed from the cockpit and connected to an alphanumeric terminal for 
system development. Therefore, the interface required to drive such a terminal 
must be included in the system. 


APPLICATIONS OF THE IN-FLIGHT COMPUTER 


Because of the flexibility of the system, there are as many applications 
of the in-flight computer as there are pilots with special interests. Two 
examples of possible applications are presented here. 

The first example is a final glide analysis. At the time that the pilot 
begins his final glide, he would enter his estimated distance to destination 
into the system. Since the sailplane's altitude is a parameter known to the 
system, a display such as the one shown in figure 8 would be output to the 
graphic display device in the cockpit. As the pilot progresses on his final 
glide and recognizes a landmark, he would again enter his estimated distance 
to destination. The system would then update the display with this new 
information, and also inform the pilot of a new speed-to-fly based on his 
actual altitude relative to his projected altitude at this point. An updated 
display is shown in figure 9. 

Some pilots may dispute the value of additional real-time information. 

The second example is one involving flight-test data acquisition. Before 
takeoff, the data to be gathered would be defined by the pilot and entered 
into the system. The data acquisition program would begin execution upon 
command and cause the storage of this data during the flight. After the pilot 
has landed, a data reduction program could be executed to display the output 
to the pilot in a useful format. 


REMARKS 


The state of electronic technology has changed drastically during recent 
years. Vast improvements in capability have caused electronic components to 
be considered for diverse applications, such as the in-flight computer dis- 
cussed in this paper. In the opening address at the 1978 Conference on 
Man-Computer Communication in Amsterdam, W. J. Doherty of IBM stated that, 

"it is not facetious to say that this industry is currently rapidly 
approaching its infancy." (ref. 7). 
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The impact of this dynamic technology will continue to be felt in the 
years to come. At this time, it is possible to construct a microprocessor- 
based in-flight computer to replace current panel instrumentation and greatly 
improve the quality and diversity of information available to the pilot. The 
prototype system, however, will be much more complicated, consume more power, 
occupy more space, provide less capability, and be more expensive than a system 
built several years from now to meet the same requirements. Hardware costs and 
power requirements are falling, even while capability per unit is rising. 
Figures 10, 11, and 12 show trends in the cost of processors and the cost and 
density of memory that substantiate this projection. 

While the capability/cost trend of the in-flight computer will probably 
not rise as sharply as that of more marketable items, such as hand-held 
calculators, the system should be within economic reach of the average soaring 
pilot within the next few years. 
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I/O SUBSYSTEM COMPUTATIONS SUBSYSTEM 



Figure 1.- Functional diagram of an in-flight computer system. 
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Figure 2.- Typical pressure transducer. 



Figure 3.- Typical 16-bit microprocessor (left) 
and memory units (right). 



Figure 4.- Typical prototype microcomputer. 
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(a) Front. 



(b) Top. 

Figure 5.- Typical digital liquid crystal display. 
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Figure 6.- Liquid crystal module capable of graphics output. 
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Figure 8.- Application: output display at beginning of final glide. 
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Figure 9.- Application: updated output display at landmarks 

during final glide. 
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Figure 11.- Cost trend of RAM memory. 


100 rxio3 



Figure 12.- Trend of amount of memory per package. 
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DESIGN OF PROPELLERS FOR MOTORSOARERS 


E, Eugene Larrabee 

Massachusetts Institute of Technology 


SUMMARY 


An efficient method has been developed for the design of propellers of 
minimum induced loss matched to an arbitrary operating point characterized by 
disc loading (thrust or power), air density, shaft speed, flight speed, and 
number of blades. A consistent procedure is outlined to predict the perfor- 
mance of these propellers under off-design conditions, or to predict the per- 
formance of propellers of general geometry. These procedures are particularly 
helpful for motorsoarer constructors and propeller builders constrained to 
use unusual powerplants under unusual conditions. The examples discussed 
include a man powered airplane, a hang glider with a 7.5 kW (10 hp) 8,000 rpm 
engine, and an a i rp I ane- I i ke motorsoarer. 


I NTRODUCT I ON 


For wings and propellers alike, there are spanwise or radial circulation 
distributions which minimize the kinetic energy loss associated with the gener- 
ation of lift or thrust. These circulation distributions give rise to simple 
induced velocity distributions which help determine efficient wing or blade 
geometry. Everyone is familiar with the elliptic span loading and the corres- 
ponding uniform downwash velocity of the vortex sheet shed by a wing of mini- 
mum induced drag, but not so many are familiar with the Betz-Prandtl (ref. I) 
or Goldstein (ref. 2) radial circulation functions and the correspond i ng uni- 
form "displacement velocity" of the hellcoidal vortex sheets shed by a propel- 
ler of minimum induced loss. The purpose of this paper is to demonstrate the 
application of these ideas to geometry determination and performance prediction 
for propellers of motorsoarers and other unusual aircraft. I have discussed 
some of these ideas in another paper (ref. 3). See Table I for notation used 
in this paper. 


THE DISPLACEMENT VELOCITY 


Consider an elementary helical vortex filament lying in an hellcoidal 
vortex sheet which forms part of the slipstream of a propeller, as shown in 
figure I. The vortex filament is constrained to move everywhere perpend icu I ar 
to itself with a velocity Wg, which is the same as the local slipstream velo- 
city, If the filament helix angle is tf)^, the axial velocity of the filament 
is Wg cos (f)g and its angular velocity is Wg sin where rg is the helix 

radius. If we were unaware of the helix angular velocity, however, we would 


suppose that it had only a displacement velocity v' = Wg/cos ({> 5 , in the same 
way that a rotating barber pole has a displacement velocity even though it 
has no axial velocity. 


Betz, in reference I, was the first to show that the condition for mini- 
mum induced loss operation of a propeller (or a windmill, for that matter) 
corresponds to radially constant displacement velocity. The axial and swirl 
components of the vortex sheet motion are then given by 


w 


axial 


v' 


2 

cos 


( I ) 


w 


swirl 


v' cos (p sin d) 
s ^s 


( 2 ) 


MOTION OF THE ENTIRE SLIPSTREAM 


Prandtl, in an appendix to reference I, pointed out that the slipstream 
fluid between the vortex sheets moves at a fraction F of the sheet velocity, 
which he evaluated by analogy with the known solution for the flow about an 
infinite array of semi - i nf i n i te plates moving perpendicular to themselves with 
velocity v, as shown in figure 2. The plate solution spacing parameter, f, is 
recalculated according to the helicoidal vortex sheet spacing and the radial 
distance from the outer edges of the sheets: 
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B /A^ + I , , r 
2 X R 


(3) 


Here X is the advance ratio 


X 5 V/fiR = (V/nD)/ir (4) 

and B is the number of blades; slipstream distortion is neglected. The corres- 
pond i ng average axial and swirl velocities at a certa in radius in the slip- 
stream are then 


w . , = F v' cos^ 6 
axial s 

(la) 

w . , = F v' cos d) sin (b 
swirl s ^s 

( 2 a) 

F = (2/tt) cos * (e 

(5) 
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THE RADIAL CIRCULATION DISTRIBUTION 


The radial circulation function correspond i ng to this minimum induced loss 
slipstream motion is found by setting the circulation about a slipstream tube 
equal to the total vorticity trailed by the blades at the corresponding radius, 
and Introducing "light loading" approximations: 


BF = 2tt r F V ' cos (j) s i n 6 
s s ^s 


r = r 
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(Os s 4. 


( 6 ) 

(7) 

( 8 ) 


(f) = tan (V/J^r) (9) 

The resulting circulation function is conveniently. written in a normalized form 
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(G for Goldstein or Glauert) where 

X = Qr/V = (r/R)/X 


( 10 ) 


(II) 


Equation 10 seems too simple to be true. Goldstein, in his doctor's 
thesis (ref. 2), verified its essential correctness, however, for propellers 
operating at low advance ratios or with many blades, where the vortex sheets 
are nearly flat, parallel, and closely spaced. The PrandtI-Betz and Goldstein 
circulation functions are compared in figure 3. It should be noted that a 
radial plot of G is identical with a radial plot of the ratio of the average 
axial slipstream velocity (Increment) to the displacement velocity. 


DETERMINATION OF THE DISPLACEMENT VELOCITY 


Following Goldstein (ref. 2) we relate the displacement velocity to the 
disc loading (thrust or torque) by resolving Joukowsky's law into two orthogo- 
nal components: 
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= pV ( I + a) FB 


(12b) 


Here the subscript "L" means that only the lift forces are being considered. 
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V 


and a' and a are the swirl and axial components of the induced velocity at the 
lifting lines. Retaining the light loading assumptions, and taking the induced 
velocities to be half the vortex sheet velocities in the developed slipstream, 
we obtain: 


I v' 

= 2 SF ♦s 
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The circulation is given by the Betz-Prandtl approximation: 


r = 


'2ttVv'' 

Bn 

V J 


(10 restated) 


The blade profile drag contributions to thrust and torque are given by 
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The radial 
cient, 


gradients of thrust 
= 2 P/pV\r 2, final I y 


coefficient, T^ = 2 T/pV^itr 2, and power coeffi- 
may be written as 


dT dl, 

c I 


dl 


dC dC 


dC dK 


dC 

dJ 




dP dJ, 


d? 


where ^ = r/R, C = v'/V, and 


dl, 


d? 


4 C G 


D/Ll 


dl2 

d? 


2 e G 


f. d/l^ 

f ' 

X 

1 

1x2 + |j 


dJ 


dK 


= 4 C G 




dJ, 


dC 


2 5 G 




'■x^ + 


(15a) 

(15b) 

( 16) 

( 17 ) 

(18) 
(19) 


288 


Equations 16, 17, 18, and 19 can be numerically integrated radially to give 
four integrals l|, I 2 , J|> and J 2 which depend only on X and B and the pre- 
sumed radial distribution of profile D/L ratio. The displacement velocity 
ratio is then easily found with these integrals and the propeller disc loading: 





/ 


4T 


c 2 
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Equations 20a and 20b are the propeller counterpart of the induced angle of 
attack of an elliptically loaded wing, Ci^/Trfb^/S) . 


If the propeller is to absorb a given amount of power, one calculates the 
power coefficient, P^, and the displacement velocity ratio, C, from equation 2 
20b; the thrust coefficient and the efficiency are then given by 
and n = Tj;/P^, respectively. The alternate procedure, when the thrust is spec- 
ified, is obvious. 


For moderately loaded propellers operating at low advance ratios, 
equations 20a or 20b may give values of C which are large compared to X. 

In this case a second approximation of the radial gradients of thrust and power 
coefficient is given by 
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Equations 21 and 22 can then be integrated radially to find better values of T^ 
and P(, appropriate to the value of C obtained from equation 20a or 20b. 
Following Theodorsen (ref. 4), one might consider a third level approximation 
in which G(X,B) is recalculated with a "vortex advance ratio", = X( 1+5/2), 
to account for slipstream distortion. 
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DETERMINATION OF THE PROPELLER GEOMETRY 


The propeller chord distribution is controlled by the choice of lift 
coefficient for the required circulation: 


ip C = p w r = P 


(G) 


This can be written as p 
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Lift coefficients must be chosen with regard to structural constraints on 
th ickness-to-chord ratios at inner radii and local Mach numbers at outer radii; 
also they must be consistent with the D/L ratios that have been used to find C- 
Some consideration must be given to off-design conditions as well; for example, 
a propeller designed for cruise can be expected to develop larger lift coeffi- 
cient increases at inner radii than at outer radii when it is operated at lower 
advance ratios, as in climbing flight. 


Traditionally propellers have been built with flat bottom airfoil sections 
such as the Clark Y. Considering the large th ickness-to-chord ratios needed 
structurally at the inner radii, and the inherent variation of lift coefficient 
with camber (proportiona I to th ickness-to-chord ratio), one can design the pro- 
peller to operate with radially constant zero angle of attack. In this case 
the propeller will have constant "true geometric pitch", given by: 


geometric 
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Modern computational airfoil theory (ref. 5) shows that the lift coefficient 
for Clark Y airfoils of varying th ickness-to-chord ratio is given by 


c^ = 0.062 t 4.21 (t/c) t 0.0971 a°; 0.07 < ^< 0.19 


0 

when they are operated at a Reynolds number of I x 10° and a Mach number of 0.2, 


The theory presented so far has assumed uniform flow at flight velocity V 
through the propeller disc at vanishingly small values of This is not a 
realistic assumption for propel lers turned by direct drive piston engines which 
are often quite large compared to the propel ler radius. If the axial velocity 
distribution, averaged around the propel ler ^i sc at radius r, is given by uV, 
it is customary to "depitch the propeller" (u(^)<l) so that the blade angle is 
given by 


3 = tan 




+ a 


(28) 


This has the effect of preserving the prescribed circulation function. 

The performance consequences of propeller-fuselage interaction are considered 
in the next section. 
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PERFORMANCE OF ARBITRARY PROPELLERS 


Unlike an untwisted eiliptical planform wing, which has elliptic loading 
over a range of angles of attack, a mi n imum- i nduced loss propeller has minimum 
induced loss loading only at its design advance ratio. An arbitrary propeller 
theory is needed to calculate its off-design point performance, or the perfor- 
mance of any general propeller. The theory given here is a radiaily graded 
momentum theory like Glauert's (refs. 6 and 7), but it will return the design 
performance of a minimum induced loss propeller when applied to the design 
conditions and geometry calcuiated by the methods described before. 


The axial and swirl components of the induced velocity at the blade ele- 
ments are found by setting the changes of axial and swirl momentum within a 
given annulus of the slipstream equal to the axial and torque loading of the 
corresponding blade elements as shown in figure 4: 
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where (see figure 4) 
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C = c„ sin (J) + c . cos (j) (33) 

X £ d 

In the absence of the propeller, the velocity in the flow field about the 
fuselage or nacelle is assumed to be given by an average axiai component u and 
an average radial component v at a distance r from the propeller shaft. 

We account for only the axial component 

u = u/V (34) 


Equations 29 and 30 can be solved for the induced 
terms of the dimensionless thrust and torque loading: 

. aC 


a 


u + a 


• 2 . 
s I n (J) 



velocity components 

Be 

27Tr 


i n 


(35) 


, I aC , 

a ' _ J X J_ 

I - a’ 4 sin<j) cos(p F 


291 


Equations similiar to these appear in Glauert's article in Durand's "Aerodyna- 
mic Theory" (ref, 7), with the vortex spacing factor F in the numerator in- 
stead of the denominator, just as his widow and R. McKinnon Wood left them. 


The induced velocity components are evaluated at each radial station by an 
iterative process outlined below: 

At each value of 5 the following are known: 


X, F, Of 3, u, c^ c^(ct), c^ c^(Cj^) 
Choose a I 

Calculate = 6 - Uj 
Calculate c^, c^ 

Calculate C , C (eqs. 32 and 33) 
y X ^ 

Calculate a and a' (eqs. 35 and 36) 


Calculate (})_ 


tan 


-I 


X (u + a ) 

C (I - a') 


Calculate cb - d) 

“l ^1 

if - 4,^^ > 0,a^< a, 

"V, ■ ♦a, ' “2 * “l 

Iterate until |0 - 4> | i s less than some small quantity, 

ot a 


n 


n 


Reta in(l),C ,C ,a,a' 
n y X n n 
' n n 


The wing theory analog of this computation is to suppose that the induced angle 
of attack at any spanwise station y of a non-e I I i pt ica I I y loaded wing of span b 
is given by 

(c/b)c. 


a. 


I 


'£ 


i nduced 


/T 


(37) 


(2y/b)' 


where c and c^ are the chord and section lift coefficient at the same station. 

The quantity /l - (2y/b)^ vanishes at y = b/2 in the same way that F vanishes 
at r = R, and it may be shown that equation 37 yields otj = C[_/7r(b2/S) for 
an elliptically loaded wing of elliptic planform. 


The values of C , C^, and a' are then integrated radially to find the 
thrust load and the power absorption of the propeller in the fuselage (or 
nacelle) flow field. These may be conveniently written in terms of coeffi- 
cients based on the shaft speed n (revolutions/sec) 

C.^ = ■ j (D = 2R) (38) 

pn D 
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T 3^5 
pn D 


dC, 3 

T _ Tr 

dC 4 

dCo ^4 

P _ 7T 

"d^ r 


2 ^3 


cos (j> 


I - a ' 
cos (j) 


roc 




(39) 

(40) 


(41 ) 


T_he a', cos ((> choice is preferred for numerical precision over the a, sin (p, 
u choice. 


The thrusting propeller is surrounded by a static pressure field with an 
appreciable axial variation, both upstream and downstream. Koning (ref. 8) 
has estimated its value: 


AP _ ^c 

Py2 “ 2 

2 


x/R 

/(x/R)^ + I 


downstream; 
tractor prope I I er 


(42) 


AP 

2 


T 

c 

2 


x/R 

/(x/R)^ + I 


upstream; 
pusher prope I I er 


(43) 


Here x 
causes 
by 


is distance downstream from the propeller. This axial pressure gradient 
the prope II er-bear i ng-f use I age or nacelle to have a buoyancy drag given 


ACp = ^ 

buoyancy ref 0 


AP 


dx, 


dx, 


(44) 


where is the reference area for drag coefficients, £ is the body length, 

and Sfj is the body cross section area at the distance X[^ behind the body nose. 
The net thrust of the propeller-body combination is then given by, 

^ref '^^Dbuoyancy 

ttR^ 

while the Installed efficiency of the propeller becomes 



C (7tA) C (V/nD) 



Figure 5 shows an application of the arbitrary propeller theory just des- 
cribed to the prediction of the performance of a scale model of a light air- 
plane propeller when tested as an "isolated" propeller, and when run at the 
nose of a representati ve fuselage. This theory Is computationally more demand- 
ing than the design theory presented in the previous section since it requires 
extensive estimates of the propeller airfoil section properties at several 
radii, a good estimate of the three dimensional flow field surrounding the 
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fuselage (or nacelle) at the propeller location, an iteration procedure to de- 
termine the induced velocities, and numerical integrations to determine blade 
loading and body buoyancy drag in the propeller pressure field. Limited ex- 
perience with it at M.l.T. shows that it gives reasonable results, and these 
are being experimentally confirmed (1979). In common with other radially 
graded momentum theories it fails to take account of the effect of circulation 
at every radial station on the downwash (or "inflow") at a particular station, 
but it is made to be consistent with the induced velocity pattern for a mini- 
mum induced loss propeller through PrandtI's analytic vortex spacing velocity 
fraction F rather than through tabulated values of Goldstein's circulation 
function. The next step up to a "prescribed" or "free" discrete vortex model 
of the "rotor" and Its "wake" is much more diffcult. 


APPLICATIONS 


(1) Man powered airplane. Here we redesign the "Gossamer Condor" propeller. 
The design conditions, summarized in figure 6, correspond to climbing flight in 
ground effect at an angle of 1°; approximately 30^ of the 53.3 N (11.8 lbs) of 
thrust is required to overcome the component of airplane weight along the 
flight path. The figure shows the radial variation of profile D/L ratio and 
the radial gradients of the integrals l|, I 2 , J|> and J 2 - The design thrust 
coefficient, T^- = 0.3175, requires a displacement velocity ratio, 5 = 0.2671, 
which corresponds to a power coefficient, P^, = 0.3914, and an efficiency, 

n = 0.8II3. The powerplant output required is 328 watts (0.44 hp). 

Since the displacement velocity ratio is moderately large, it is worth- 
while to recalculate the thrust coefficient and the power coefficient using 
eqs. 21-24. The results are summarized in figure 7, which also compares the 
propeller geometry determined by the methods of this paper with the geometry 
actually employed. The agreement of blade angles is very good, especially 
when one takes into account the difference between the zero lift angles of the 
Clark Y airfoils assumed in the design calculations and the Stratford pressure 
recovery airfoils used on the "Gossamer Condor". In my opinion the propeller 
calculated here would be more efficient than the one actually flown. 

(2) Powered hang glider. Soarmaster, Inc. supplies a powerpack consisting of 
a West Bend (Chrysler) two stroke, single cylinder engine developing 7.46 kW 
(10 hp) at 10,000 rpm, a centrifugal clutch, a chain and sprocket reduction 
gear, and an extension shaft turning a pusher propeller. This is a suitable 
powerplant for hang gliders of 12 m (40 ft) span; figure 8 presents the options 
available for propellers intended to absorb the engine power at a flight speed 
of 13 m/sec (30 mph). The diameter of the direct drive propeller is limited 

to 690 mm (27 in) by a tip Mach number of 0.85; its efficiency is very poor 
because of the excessive disc loading. Gear reductions and larger propellers 
lead to progressive improvements in performance. Figure 9 gives the geometry 
of the largest propeller considered, a 1372 mm (54 in) diameter propeller 
turned at 1946 rpm by a 9:37 sprocket pair driven at 8000 rpm. It has 617 mm 
(24 in) nominal pitch, and the typical wide root chord - narrow tip chord 
geometry of a propeller matched to a low advance ratio; this is in spite of a 
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design lift coefficient of I near the hub and 0.5 at the tip. Soarmaster 
supplies two propeller options: a 1067 x 483 mm (42 x 19 in) or an 1118 x 356 
mm (44 x 14 in) "laminar" airscrew, both of fiber reinforced plastic. 

Table 2 summarizes the propeller parameters covered in this study: when 
two values are given for Tc, Pc. and n, the second set corresponds to the 
improved velocity polygon geometry corresponding to eqs. 21-24; note the rela- 
tively good agreement, even for ^ values of more than I. Examination of 
figure 9 and Table 2 suggests that still larger propel lers and larger reduction 
ratios would improve climbing performance; this has to be balanced against the 
weight penalty and the reduction of ground clearance at the tail. 

(3) Motorsoarer. The Ryson ST-100 is a 17.58 m (57.67 ft) span two seated 
aircraft, with a flying mass of 748.4 kg (1650 lbs), fitted with a Hoffmann 
HO-V-62 propeller of 1.7 m (67 in) diameter. This propel ler has a low pitch 
setting, a high pitch setting, and can be feathered for glider mode operation. 
Figure 10 shows three design points which might be considered In the selection 
of such a propeller: sea level climbing performance at 40 m/sec (90 mph); 
sea level top speed at about 68 m/sec (152 mph); and cruise at 75% power at 
full throttle at 1981 m (6500 ft) altitude and 65 m/sec (145 mph). The circled 
points show the performance that may reasonably be expected from minimum 
induced loss propel lers designed for each of these flight conditions by the 
methods of the paper. 

Figure II shows how a compromise propeller may be designed which will 
give nearly this performance at two of these points. The displacement veloci- 
ties are calculated assuming minimum induced loss loading and a somewhat pessi- 
mistic radial distribution of D/L ratio. Blade lift coefficients are assigned 
at 5 = 0.3 and ^ = 0.7 so that the blade chord to radius ratio, c/R, as given 
by equation 26, is the same for both flight conditions. The c/R ratio Is then 
calculated at other radii, assuming a linear radial variation of cj_. 

Reasonable assumptions are then made about the radial variation of thickness 
to chord ratio, t/c, to give the radial variation of blade angle (eq. 28). 

The compromise c/R ratio and blade twist, AS, are then chosen to minimize 
differences between the two conditions. In general, highly loaded, low 
advance ratio flight conditions demand high lift coefficients near the hub 
Betz (ref. 9) was of the opinion that Coriolis forces within the rotating 
blade boundary layer favored such a distribution. 
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TABLE I 

SYMBOLS AND NOTATION (follows Glauert; ref. 8) 

aV axial component of induced velocity (m/sec) 

'J^r rotational (swirl) component of induced velocity (m/sec) 
number of blades 
b wing span (m) 

c blade chord (m) 

C(j section (profile) drag coefficient 

C£ section (profile) lift coefficient 

C[_ wing lift coefficient 

Cp power coefficient (Cp = P/pn^D^) 

Cj thrust coefficient (Cj s T/pn^D^) 

Cx blade element torque load coefficient 

Cy blade element thrust load coefficient 

D^ drag; also propeller diameter (m) 

F slipstream velocity fraction (eg. 5) 

f vortex sheet spacing parameter (eg. 3) 

G circulation radial distribution function (eqs. 10,11) 

I|,l 2 thrust loading integrals (eqs. 16,17) 

J I ,J2 power loading integrals (eqs. 18,19) 

L lift 

n revolutions per second 

P shaft power (kW) 
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Pc power coefficient (Pc = 2 P/pV3ttR2) 

I Q propeller shaft torque (Nm) 

R propeller tip radius (m) 

r propeller general radius (m) 

S wing plan or fuselage cross section area (m^) 

1 T thrust (N) 

thrust coefficient (Tc = 2 T/pV2itr2) 

£ axial velocity of fuselage flow field at r (m/sec) 

uV average axial velocity at r (m/sec) 

V flight velocity (m/sec) 

V radial velocity of fuselage flow field at r (m/sec) 
v' displacement velocity (m/sec); see fig. I_^ 

W resultant velocity at blade element^CW + w)(m/sec) 

w induced velocity at blade element (w = aV t a'fir)(m/sec) 

Wg slipstream velocity ( i ncrementa I ) (m/sec) 

X velocity ratio (x = flr/V) 

y spanwise location (m) 

a section angle of attack (rad); a° (degrees) 

3 section blade setting angle (rad); 3° (degrees) 

r circulation (m^/sec) 

Z, displacement velocity ratio (? = v'/V) 

n efficiency (q = Tc/Pc = (V/nD)Cj/Cp) 

X advance ratio (X = V/fiR) 

C radius ratio (^ = r/R) 

p air density (kg/m^) 

a blade solidity (a = Bc/27rr) 

(p helix angle (rad); cj) = 3 - ot 

n shaft speed (rad/sec) 


TABLE 2 


POWERED HANG GLIDER PROPELLERS 
V = 13.41 m/sec (30 mph) 

7.457 kW (10 hp) @ 8000 engine rpm 
p = 1.225 kg/m^ (760 mm Hg, I5°C) 


Gear 

Ratio 

2R 

m 

Pc 


Tc 

n 

p i tch 
d i ameter 

7TX 

1 : 1 

0.690 

13.500 

2.544 

4.650 

0.344 

0.356 

0.146 

9:27 

1 .000 

6.426 

1.778 

2.888 

0.449 

0.507 

0.288 

9:27 

1 .219 

4.323 

3.881 

1.283 

2.179 
2. 103 

0.504 

0.541 

0.361 

0.220 

9:37 

1.372 

3.348 
3. 145 

1 . 151 

1 .840 
1.821 

0.550 

0.579 

0.450 

0.301 
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no. / VORTEX SHEET MOn ON 

M/MMOM //VDUCBD/OSS PROPEILER 



E/G. 2 PRANDTIS VELOCITY FACTOR , F 
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f/0,-5 RADIAL CIRCULATION D/5TRIBUr/OH5 

A1/A/fMUM /A/DUCrO lOS5 P/POPSJJLSA^S 
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FKS.5 radially GRADDD MOMfMim THDOnr 



FL6. PROPDLLDR D£5K^P FU/VCT/ORS 
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no. 9 PROP£Ll£R OEOPJPTRY 
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